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We provide a theoretical perspective on the glass transition in molecular liquids at thermal equi- 
librium, on the spatially heterogeneous and aging dynamics of disordered materials, and on the 
rheology of soft glassy materials. We start with a broad introduction to the field and emphasize 
its connections with other subjects and its relevance. The important role played by computer 
simulations to study and understand the dynamics of systems close to the glass transition at 
the molecular level is spelled out. We review the recent progress on the subject of the spatially 
heterogeneous dynamics that characterizes structural relaxation in materials with slow dynamics. 
We then present the main theoretical approaches describing the glass transition in supercooled 
liquids, focusing on theories that have a microscopic, statistical mechanics basis. We describe both 
successes and failures, and critically assess the current status of each of these approaches. The 
physics of aging dynamics in disordered materials and the rheology of soft glassy materials are 
then discussed, and recent theoretical progress is described. For each section, we give an extensive 
overview of the most recent advances, but we also describe in some detail the important open 
problems that, we believe, will occupy a central place in this field in the coming years. 
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I. INCIPIT 

Gl asses belon 
ter dTabor. 199ll 



to a well-known state of mat- 
we easily design glasses with de- 
sired mechanical or optical properties on an industrial 
scale, they are widely present in our daily life. Yet, a 
deep microscopic understanding of the glassy state of 
matt er remains a challenge for condensed matter ph ysi- 
cists (|Angell. 19951 iDebenedetti and Stillinger. 200l . 

Glasses share similarities with crystalline solids since 
they are both mechanically rigid, but also with liq- 
uids because they both have similar disordered struc- 
tures at the molecular level. It is mainly this mixed 
character that makes t hem fascinating even to non- 
scientists ( Zanotto. 19981 ) . Given that glasses are neither 
normal liquids nor standard solids, they are quite often 
not described in any detail by standard textbooks. For 
instance, glasses are not described in textboo ks on con- 
densed matter (IChaikin and Lubenskv. 2 0001). or solid 
state physics ( Ashcroft and Mermin. 19761) . they only 
made it in the latest edition of the refe rence textbook 
on liquids ( Hansen and McDonald. 20061 ). while statisti- 
cal mechanics textbooks usually culminate with a presen- 
tation of our current understanding of phase transitions 
i n pure material s using renorm alization group concepts 
(jChandler. 19871 ISethna. 2006T) . leaving out disordered 
systems. 

As we shall describe in detail in this review, modern 
statistical mechanics approaches to the glass transi- 
tion involve good knowledge of advanced liquid state 
theory, field theory, renormalization group, solution 
of lattice models, percolation, replica calculations, 
and concepts de veloped for f ar-from-equilibrium , 
driven systems (Young. 19981 : iBarrat et al. . 2003 
iBinder and Kob. 20051 ). These development are all 
posterior to the mid-70's important breakthroughs on 
phase transitions: the c anonical spin glass Hamiltonia n 
was introduced in 1975 ( Edwards and Anderson. 19751 ). 
to be solved in infinite dimension on ly several years 
later (jParisi. 19801: iMezard et al. 1988(1 . mode-couplin g 
theory was developed in the mid-80's ( Gotze. 19991 ). 
just be fore kinetic lattice glass model s were intro- 
duced (jFredrickson and Andersen. 1984al) . The aging 
and rheology of disordered systems such as spin glasses 
or soft materials emerged as broad research fields during 
the 90's. In this paper, we shall review the fruits that 
have grown out of these important seeds. Given none of 
these advances has allowed the derivation of a complete, 
well-accepted theory of amorphous media, we present a 
large number of different approaches. We try to discuss 
both successes and failures, we explain similarities and 
differences between them, and we present the current 
status of each approach. Thus, the article takes at times 
a somewhat subjective tone. 

A glass can be obtained by cooling the temperature of a 
liquid below its glass temperature, T g . The quench must 
be fast enough that the more standard first order phase 
transition towards the crystalline phase is avoided. The 




FIG. 1 Glassy phases occur at low temperature or 
large density in many different systems spanning a broad 
range of lengthscales, such as atomic (top left, atomic 
force spectroscopy image of an alloy linear size 4.3 
nm (jSugimoto et al.. 20071 ) 1. colloidal (top right) systems, 
foams (bottom left, a beer foam with bubbles of submillimeter 
size) and granular materials (bottom right, a fertilizer made 
of millimeter size grains). 



glass 'transition' is not a thermodynamic transition at all, 
since T g is only empirically defined as the temperature 
below which the material has become too viscous to flow 
on a 'reasonable' timescale (and it is hard to define the 
word 'reasonable' in any reasonable manner). Therefore, 
T g does not play a fundamental role, as a phase tran- 
sition temperature would. It is simply the temperature 
below which the material looks solid. When quenched in 
the glass phase below T g , liquids slowly evolve towards 
an equilibrium state they cannot reach on experimental 
timescales. Physical properties are then found to evolve 
slowly with time i n far from equ ilibrium states, a process 
known as 'aging' ( Struik. 19781 ). 



The subject of the glass transition has quite broad im- 
plications. A material is said to be 'glassy' when its typ- 
ical relaxation timescale becomes of the order of, and 
often much larger than, the typical duration of an ex- 
periment or a numerical simulation. With this generic 
definition, a large nu mber of syste ms can be considered 
as glassy materials ( Young. 19981 ). One can be inter- 
ested in the physics of liquids (window glasses are then 
the archetype), in 'hard' condensed matter (for instance 
type II superconductors in the presence of disorder such 
as high-T c superconducting materials), in charge density 
waves or spin glasses, in 'soft' condensed matter with nu- 
merous complex fluids such as colloidal assemblies, emul- 
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sions, foams, but also granular materials, proteins, etc. 
Glass physics thus covers a remarkably broad range of 
timcscales and lcngthscales, as illustrated in Fig. [T] All 
these materials exhibit, in some part of their phase dia- 
grams, some sort of glassy dynamics characterized by a 
very rich phenomenology with effects such as aging, hys- 
teresis, creep, memory, effective temperatures, rejuvena- 
tion, dynamic heterogeneity, non-linear response, etc. 

These long enumerations explain why this research 
held has received increasing attention from physicists in 
the last two decades. 'Glassy' topics now go much beyond 
the physics of simple liquids (glass transition physics) 
and models and concepts developed for one system of- 
ten find application s elsewhere in physic s, from algorith- 
mics to biophysics ( Mczard et at, 2007t) . Motivations to 
study glassy materials are numerous. Glassy materials 
are everywhere around us and therefore obviously at- 
tract interest beyond academic research. At the same 
time, the glass conundrum provides theoretical physicists 
with deep fundamental questions since standard statisti- 
cal mechanics tools are sometimes not sufficient to prop- 
erly account for the glass state. Additionally, simulating 
in the computer the dynamics of microscopically realistic 
material on timescales that are experimentally relevant is 
not an easy task, even with modern computers. Finally, 
the field is constantly stimulated by new, and sometimes 
quite beautiful, experimental developments to produce 
new types of disordered materials, or to obtain more mi- 
croscopic information on the structure and dynamics of 
glassy systems. 

The outline of the article is as follows. Sec. [XT] pro- 
vides a broad introduction to glassy materials. The is- 
sue of dynamic heterogeneity is tackled in Sec. IIIIl while 
the main theoretical perspectives, characterized by a mi- 
croscopic, statistical mechanics basis are summarized in 
Sec. IIVI Aging and noncquilibrium phenomena occupy 
Sec. El Finally, we present a set of general and conclud- 
ing remarks in Sec. I VII 



II. A BROAD INTRODUCTION ABOUT GLASSES 



A. Some phenomenology 



1. The basic phenomenon 



A vast majority of liquids (molecular liquids, polymeric 
liquids, etc.) form a glass if cooled fa st enough in o rder to 
avoid the crystallization transition (jAngell. 19951 ). Typ- 
ical values of cooling rate in laboratory experiments are 
0.1 — 100 K/min. The metastable phase reached in this 
way is called 'supercooled phase'. In this regime the typ- 
ical timescales increase in a dramatic way and they end 
up being many orders of magnitudes larger than micro- 
scopic timescales at T g , the glass transition temperature. 

For example, around the melting temperature T m , the 
typical timcscalc r a on which density fluctuations relax, 
is of the order of y 'ma 2 /fcsT ', which corresponds to few 



picoseconds (m is the molecular mass, T the tempera- 
ture, ks the Boltzmann constant [which will often be set 
to unity in the later theoretical sections] , and a a typical 
distance between molecules). At T g the typical timescale 
has become of the order of 100 s, i.e. 14 orders of mag- 
nitude larger! This increase is even more remarkable be- 
cause the corresponding temperature decrease is, as a 
rule of thumb, about \T m — not a large value if one con- 
siders fcgT m as the typical energy scale. The increasing of 
T a is accompanied by a concomitant increase of the shear 
viscosity r\. This can be understood by a simple Maxwell 
model in which rj and T a are related by rj = GooT a , where 
Goo is the instantaneous (elastic) shear modulus which 
does not vary considerably in the supercooled regime. In 
fact, viscosities at the glass transition temperature are 
of the order of 10 12 Pa.s. In order to grasp how viscous 
this is, recall that the typical viscosity of water at ambi- 
ent temperature is of the order of 10~ 2 Pa.s. How long 
would one have to wait to drink a glass of water with a 
viscosity 10 14 times large r, or how long w ould it take for 
cathedral glasses to flow ( Zanotto. 19981) ? 

As a matter of fact, the temperature at which the liquid 
does not flow anymore and becomes an amorphous solid, 
called a 'glass', is protocol-dependent. It depends on the 
cooling rate and on the patience of the people carrying 
out t he experiment: solidi ty is a timescale-dependent no- 
tion ( Sausset et al. 2O10h . Pragmatically, T g is defined 
as the temperature at which the shear viscosity is equal 
to 10 13 Poise (also 10 12 Pa.s). 



2. Kinetic fragility 

The increase of the relaxation timescale of supercooled 
liquids is remarkable not only because of the large num- 
ber of decades involved but also because of its tempera- 
ture dependence. This is vividly demonstrated by plot- 
ting the logarithm of the viscosity (or the relaxation 
time) as a functio n of T a /T, as in Fig. [2] This is called 
the 'Angell' plot (jAngell. 19951) and is very helpful in 
classifying supercooled liquids. A liquid is called strong 
or fragile depending on its position in the Angell plot. 
Straight lines correspond to 'strong' glass-formers and to 
an Arrhenius behaviour, 



, exp 



E 



(i) 



In this case, one can extract from the plot an effective ac- 
tivation energy, E, suggesting quite a simple mechanism 
for relaxation by 'breaking' locally for instance a chemical 
bond. The typical relaxation time is then dominated by 
the energy barrier to activate this process and, hence, has 
an Arrhenius behaviour. Window glasses generically fall 
in this category. The terminology 'strong' and 'fragile' is 
not related to the mechanical properties of the glass but 
was introduced in relation to the evolution of the short- 
range order close to T g . Strong liquids, such as Si02, 
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FIG. 2 Arrhenius plot of the relaxation time of several glass- 
forming liquids approaching the glass temperature T g . For 
'strong' glasses, r Q increases in an Arrhenius manner as tem- 
perature is decreased, logr a ~ _E/(fcsT), where E is an ac- 
tivation energy and the plot is a straight line. For 'fragile' 
liquids, the plot is bent and the effective activation energy 
increases when T is decreased towards T g . Timescales acces- 
sible to numerical simulations are indicated. BKS: numerical 
model of silica; LJ: numerical model of a binary Lennard- 
Jones mixture; GLY: glycerol; OTP: ortho-terphenyl; SAL: 
salol; PC: propylene carbonate; DEC: decaline. 



typically have a locally tetrahedric structure which per- 
sists both below and above the glass transition, contrary 
to fragile liquids. Nowadays, the notion of fragility ap- 
pears in connection with the temperature evolution of 
transport properties. 

If one tries to define an effective activation energy 
for fragile glass-formers using the slope of the curve in 
Fig. [2j then one finds that this energy scale increases 
when the temperature decreases, a 'super- Arrhenius' be- 
haviour. This increase of energy barriers immediately 
suggests that the glass formation is a collective phe- 
nomenon for fragile supercooled liquids. Support for this 
interpretation is provided by the fact that a good fit of 
the relaxation time or the viscosity is given by the Vogel- 
Fulcher-Tamman law (VFT): 



t q = t exp 



DTn 



[T-T 



(2) 



which suggests a divergence of the relaxation time, and 
therefore perhaps a phase transition of some kind, at a 
finite temperature To . A smaller D in the VFT law corre- 
sponds to a more fragile glass. Note that there are other 
com parably good f its of these curves, such as the Bassler 
law ([Bassler. 19871) . 



r a = t exp 



I — 

T 



(3) 



that only lead to a divergence at zero temperature. An- 
other fit proposed and tested recently consists in replac- 
ing l/T by l/T - l/T on in the Bassler law. This form, 



Substance 


OTP 


2-MTH 


N-PROP 


3-BP 


12PD 


T g 


246 


91 


97 


108 


172 


T 


202.4 


69.6 


70.2 


82.9 


114 


T K 


204.2 


69.3 


72.2 


82.5 


127 


Tk/T 


1.009 


0.996 


1.028 


0.995 


1.11 



TABLE I Values of glass transition temperature, VFT singu- 
larity and Kauzmann tempera tures for five supercooled liq- 
uids (|Richert and Angell. 199Sh . OTP: o-terphenyl; 2-MTH: 
2-methyltetra-hydrofuran; N-PROP: n-propanol; 3-BP: 3- 
bromopentane. 12PD: 1-2 prop-diol. 



where T on represents some onset temperature for slow 
dynamic s, is motivated in par t by dynamical facilitation 
models (|Elmatad et ai. 20091) . see Sec. IIV.C.2I Actu- 
ally, although the relaxation time increases by 14 orders 
of magnitude, the increase of its logarithm, and therefore 
of the effective activation energy is very modest, and ex- 
perimental data do not allow one to unambiguously de- 
termine the true u nderlying functional l aw without any 
reasonable doubt ( Hecksher et al. 20081 ). Several com- 
parisons between fits aimed at determining which one 
works better can be found in the l iterature, see e.g. the 
recent work ( Elmatad et al. 20101 ). Although these are 
certainly useful, one should not forget that since the in- 
crease of the effective activation energy is very modest, 
pre-asymptotic effects to the 'true' limiting behavior can 
play an important role. Hence, comparisons of simple 
fits with a small number of parameters could be mislead- 
ing. For this and other reasons, physical interpretations 
in terms of a finite temperature phase transition must 
always be taken with a grain of salt. It is recommended 
to use the same grain of salt to deal with fits supposedly 
demons trating the absence o f finite temperature s ingu- 
larities (jHecksher et al, 20081 ; [Elmatad et al. 20091 ). 



3. Thermodynamic aspects 

There are additional experimental facts that shed some 
light and might reinforce the interpretation of data in 
terms of a finite temperature singularity. Among them, is 
an empirical connection found between kinetic and ther- 
modynamic behaviours. Consider the part of the entropy 
of the liquids, S exc , which is in excess compared to the 
entropy of the corresponding crystal. Once this quantity, 
normalized by its value at the melting temperature, is 
plotted as a function of T, a remarkable connection with 
the dynamics, in particular the VFT law, emerges (see 
([Martinez and Ange ll. 2001 ) for a compilation of e xperi- 



mental data and ( Debenedetti and Stillinger. 200ll ) for a 
discussion). As for the relaxation time, one cannot follow 
this curve below T g in thermal equilibrium. However, ex- 
trapolating the curve below T g apparently indicates that 
the excess entropy vanishes linearly at some finite tem- 
perature, called T/f, which is very close to zero for strong 
glasses and, generically, very close to To, the temperature 
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at which a VFT fit diverges. This coincidence is quite 
remarkable: for materials with glass transition temper- 
atures that vary from 50 K to 1000 K the ratio Tk/Tq 
remains c lose to 1. Some examples are provided in Ta- 
ble |H see (jRichert and Angell. 19981) for a more extensive 
list. For the majority of liquids the ratio is close to one 
up to few percent. Note, however, there are some liquids 
where Tk and To differ by as much as 20%, and so a 
perfect correlation between the two temp eratures is not 
established experimentally ( Tanaka. 20031 ). 

The chosen subscrip t for Tk stands for Kauz- 
mann ( Kauzmann. 19481 ) who recognized Tk as a very 
important temperature for the physics of glasses. Kauz- 
mann further suggested that some change of behaviour 
(phase transition, crystal nucleation, etc.) must take 
place above Tk, because below Tk the entropy of the 
liquid, a disordered state of matter, becomes less than 
the entropy of the crystal, an ordered state of matter. 
This situation, that seemed perhaps paradoxical at that 
time, is in fact not a serious problem. There is no gen- 
eral principle that would constrain the entropy of the 
liquid to be larger than that of the crystal. As a mat- 
ter of fact, the crystallization transition for hard spheres 
takes place precisely because the crystal becomes the 
state with the largest entropy at sufficiently high den- 
sity ( Alder and Wainwright. 19621 ). 

On the other hand, the importance of Tk stands, par- 
tially because it is experimentally very close to To. Ad- 
ditionally, the quantity 5 CXC which vanishes at Tk, is 
thought to be a proxy for the so-called configurational 
entropy, S c , which quantifies the number of metastable 
states (actually, its logarithm, see below). A popu lar 
physical picture due to Goldstein ( Goldstein. 19691 ) is 
that close to T g the system explores a part of the en- 
ergy landscape (or configuration space) which is full of 
minima separated by barriers that increase when tem- 
perature decreases. The dynamic evolution in the energy 
landscape would then consist in a rather short equili- 
bration inside the minima followed by 'jumps' between 
different minima that are well-separated in time. At T g 
the barriers have become so large that the system re- 
mains trapped in one minimum, identified as one of the 
possible microscopic amorphous configurations of a glass. 
Following this interpretation, one can split the entropy 
into two parts. A first contribution is due to the fast re- 
laxation inside one minimum, a second one counting the 
number of metastable states, S c = log Ar meta stabic, which 
is called the 'configurational' entropy. Assuming that 
the contribution to the entropy due to the 'vibrations' 
around an amorphous glass configuration is not very dif- 
ferent from the entropy of the crystal, one finds that 
S^xc ~ S c . Within this approximation Tk corresponds 
to the temperature at which the configurational entropy 
vanishes. Since the configurational contribution to the 
specific heat is given by TdS oxc /dT, a linear vanishing 
of 5 xc near Tk (suggested both by experimental obser- 
vations and theoretical arguments, see Sec. IIV[) would 
lead to a discontinuity (a downward jump) of the spe- 



cific heat and thus to a thermodynamic phase transition. 
Note, however, that the above assumptions should not be 
tak en for granted, see for ins t ance the recent discussion s 
in (j Angell and Borick. 2002 1 : iDvre. 20061: IWvart. 2009T) . 
Furthermore, locating the transition temperature re- 
quires an extrapolation that is not well controlled, as 
is also required for the relaxation time. 



4. Static and dynamic correlation functions 

At this point the reader might have reached — despite 
our numerous warnings — the conclusion that the glass 
transition may not be such a difficult problem: there are 
experimental indications of a diverging timcscale and a 
concomitant singularity in the thermodynamics. It sim- 
ply remains to find static correlation functions displaying 
a diverging correlation length related to the emergence of 
'amorphous order', which would indeed classify the glass 
transition as a standard second order phase transition. 
Remarkably, this remains an open and debated question 
despite several decades of research. Simple static corre- 
lation function are quite featureless in the supercooled 
regime, notwithstanding the dramatic changes in the dy- 
namics. A simple static quantity is the structure factor 
defined by 



S(q) 



PqP-q 



where the Fourier component of the density reads 



Pi 



N 

E- 

3=1 



(4) 



(5) 



with N is the number of particles, V the volume, and Yj is 
the position of particle j. The structure factor measures 
the spatial correlations of particle positions, but it does 
not show any diverging peak in contrast to what happens, 
for example, at the liquid-gas critical point where there 
is a divergence at small q. A snapshot of a supercooled 
liquid configuration in fact just looks like a glass config- 
uration, despite their widely different dynamic proper- 
ties. More co mplicated static cor relation functions have 
been studied ( Debenedetti. 19961 ). especially in numeri- 
cal work, but until now there are no strong indications of 
a diverging static lengthscale, although th i s issue is con- 
stantly debated dMenon and Nagel 19951 iNelson. 2002t 
iFernandez. 20061: ICavagna. 20091: iTanaka et ai. 20101) 
Recent results do suggest that it is possible to identify 
some growing static lengthscales. We will come back to 
this point in Sec. IVII and in other sections. 

The difficulty in finding a signature of the glass transi- 
tion disappears if one focuses on dynamic correlations or 
response functions. For instance, a dynamic observable 
studied in light and neutron scattering experiments is the 
intermediate scattering function, 



(6) 



F(q,t) = (-p q (t)p_ q (0) 





t (nsec) 

FIG. 3 Temperature evolution of the intermediate scatter- 
ing function normaliz ed by its value at ti me equal to zero for 
supercooled glycerol (|Wuttke et al., 19961 ). Temperatures de- 
crease from 413 K to 270 K from left to right. The solid lines 
are fit with a stretched exponential with exponent /3 = 0.7. 
The dotted line represents another fit with /3 = 0.82. 
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FIG. 4 Temperature evolution of the dielectric susceptibility 
of the glass-former benzopheno ne measured over more than 
10 decades of relaxation times (|Pardo et al, 2007t ). Dynam- 
ics slows down dramatically as temperature is decreased and 
relaxation spectra become very broad at low temperature and 
reveal the existence of additional 'secondary' relaxation pro- 
cesses. 



Different F(q, t) me asured by neutron s cattering in su- 
percooled glycerol ( Wuttke et al., 1996h are shown for 
different temperatures in Fig. [3J These curves suggest 
a first, rather fast (and hence not accessible in this ex- 
periment), relaxation to a plateau followed by a second, 
much slower, relaxation. The plateau is due to the frac- 
tion of density fluctuations that are frozen on intermedi- 
ate timescales, but eventually relax during the second re- 
laxation. The latter is called 'alpha- relaxation', and cor- 
responds to the structural relaxation of the liquid. This 
plateau is akin to the Edwards-Anderson order param- 
eter, qEA, defined for spin glasses which measures the 
fraction of frozen spin fluctuations ( Mezard et al., 1988t 
iBinder and KobT 20051) . Note that qEA continuously in- 
creases from zero below the spin glass transition. Instead, 
for structural glasses, a finite plateau already appears 
above any putative transition. Figures [4j [6] and [13] be- 
low contain more illustration of the different time regimes 
observed in time correlators. 

The intermediate scattering function can be probed 
only on a relatively small regime of temperatures. In 
order to track the dynamic slowing down from micro- 
scopic to macroscopic timescales, other correlators have 
been studied. A popular one is obtained from the dielec- 
tric susceptibility, which is related by the fluctuation- 
dissipation theorem to the time correlation of polariza- 
tion fluctuations. It is generally admitted that different 
dynamic probes reveal similar temperature de pendencies 
for the relaxation time ( Jakobsen et al., 20051 ) . The tem- 
perature evolution of the imaginary part of the dielectric 
susceptibility, e"(w), is shown in Fig. @] for glass- former 
benzophenone, where a a very wi de temperature win- 
dow is covered (jPardo et al., 2007h . At high tempera- 
ture, a good representation of the data is given by a 
Debye law, e(w) = e(oo) + Ae/(1 + ioJT a ), which cor- 



responds to an exponential relaxation in the time do- 
main. When temperature is decreased, however, the re- 
laxation spectra become very broad and strongly non- 
Debye, which is the frequency analog of the strech- 
ing of the relaxation observed in the time-domain. In- 
deed one particularly well-known feature of the spec- 
tra is that they are well fitted, in the time domain, 
for times corresponding to the alpha-relaxation with a 
stretched exponential, exp(— (t/r a )^). In the Fourier do- 
main, forms such as the Havriliak-Ncgami law are used, 
e(ui) = e(oo) + Ae/(1 + (iajr Q ) Q ) 7 , which generalizes 
the Debye law. The exponents /3, a and 7 depend in 
general on temperature and on the particular dynamic 
probe chosen, but they capture the fact that relaxation 
is increasingly non-exponential when T decreases to- 
wards T g . A connection was empirically established be- 
tween fragility and degree of non-exponentiality, more 
fragile li quids being characterized by broa der relaxation 
spectra ( Debenedetti and Stillinger. 200ll ) , although th e 
correlation is — again — not very solid ( Heuer. 2008|) . 
The data of Fig. @] also make it clear that the 
relaxation spectra are actually quite complex and 
characterized by one or several secondary processes, 
that h ave been quite extensively s tudied experimen- 
tally dJohari and Goldstein. 19701 : iDixon et al., 19901 
iBlochowicz et al., 20061) . 

To sum up, there are many remarkable phenomena 
that take place when a supercooled liquid approaches 
the glass transition. Striking ones have been pre- 
sented, but many others have been left out for lack 
of space (lAngell. 19951: iDebenedetti and Stillinger. 20011: 



iDebenedetti. 19961: IBinder and Kob. 20051) . We havedis- 
cussed physical behaviours, relationships or empirical 
correlations observed in a broad class of materials. This 
is quite remarkable and suggests that there is some inter- 
esting physics to be studied in the problem of the glass 
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transition, which we see as a collective (critical?) phe- 
nomenon which is relatively independent of microscopic 
details. 



B. Other 'glasses' in science 

We now introduce some other systems whose phe- 
nomenological behaviour is close or, at least, related, to 
the one of glass-forming liquids, showing that glassiness is 
truly ubiquitous. It does not only appear in many differ- 
ent physical situations but also in more abstract contexts, 
such as computer science. 



1. The jamming transition of colloids and grains 

Colloidal suspensions co nsist of big particles suspended 
in a solvent ( Larson. 19991 ) . The typical radii of the parti- 
cles are in the range R = 1 — 500 nm. The solvent, which 
is at equilibrium at temperature T, acts as a source of 
noise on the particles whose short-time dynamics is bet- 
ter described as being Brownian rather than Newtonian. 
The microscopic timcscale for this diffusion is given by 
t = R 2 /D where D is the short-time self-diffusion co- 
efficient. Typical values are of the order of r ~ 1 ms, 
and thus are much larger than the ones for molecular liq- 
uids (in the picosecond regime). The interaction poten- 
tial between particles depends on the systems, and this 
large tunability makes colloids very attractive objects for 
technical applications. 

A particularly relevant case, on which we will focus in 
the following, is a purely hard sphere potential, which is 
zero when particles do not overlap and infinite otherwise. 
In this case the temperature becomes irrelevant, apart 
from a trivial rescaling of the microscopic timcscale. Col- 
loid al hard sphere s systems have been intensively stud- 
ied ( Larson. 19991 ) in experiments, simulations and the- 
ory varying their density p, or their volume fraction 4> = 
T-nR 3 p. Hard spheres display a fluid phase from = to 
intermediate volume fractions, a freezing-crystallization 
transition at (f> ~ 0.494, and a melting transition at 
<f> ~ 0.545. Above this latter value the system can 
be compressed until the close packing point <fi ~ 0.74, 
which corresponds to the FCC crystal. Interestingly for 
our purposes, a small amount of polydispersity (particles 
with slightly different sizes) efficiently prevents crystal- 
lization. In this case, the system can be more easily 'su- 
percompressed' above the freezing transition without nu- 
cleating the crystal, at least on experimental timescales. 
In this regime the relaxation timescale in creases very 
rapidly with (j) ( Pusev and van Megen. 19861 ). At a pack- 
ing fraction <f> g ~ 0.57 — 0.59 it becomes so large com- 
pared to typical experimental timescales that the system 
does not relax anymore: it is 'jammed'. This 'jamming' 
transition is obviously reminiscent of the glass transition 
of molecular systems. In particular, the location <p g of 
the colloidal glass transition is as ill-defined as the glass 



temperature T g . 

Actually, the phenomena that take place increasing the 
volume fraction are analo gous to the ones seen in molec - 
ular supercooled liquid (IPusev and van Meeen. 19861) : 
the viscosity increases very rapidly and can be fit- 
ted ( Cheng et al, 20021 ) by a VFT law in density as 
in Eq. ([2]); dynamical correlation functions display a 
broad spectrum of timescal es and develop a plateau 
( Pusev and van Megen. 19871 ) . no static growing corre- 
lation length has been found, etc. Also the phenomenon 
of dynamic heterogeneity that will be addressed in 
Sec. IIII| seems similar in c o lloids and atomic systems 
(jKegel and Blaaderen. 200(1 IWeeks et al, 2000f) . How- 
ever, it is important to underline a major difference: be- 
cause the microscopic timescale for colloids is so large, 
experiments can only track at best the first 5-6 decades 
of slowing down ( Brambilla et al. 20091 ). A major conse- 
quence is that the comparison between the glass and col- 
loidal transitions must be performed by focusing in both 
cases on the first 5 decades of the slowing down, which 
corresponds to relatively high temperatures in molecular 
liquids. Understanding how much and to what extent the 
glassiness of colloidal suspensions is related to the one of 
molecular liquids remains an active domain of research. 

The glassiness of driven granular media has recently 
been thoroughly analyzed. Grains are macroscopic 
objects and, as a consequence, do not have any thermal 
motion. A granular material is therefore frozen in a 
given c onfiguration if no energy is injected into the 
system ( Jaeger et al, 19961 ). However, it can be forced 
in a steady state by an external drive, such as shearing 
or tapping. The dynamics in this steady state shows 
remarkable similarities (and differences) with simple 
fluids. The physics of granula r materials is a very 
wide subject (j Jaeger et al, 19961 ). In the following we 
only address briefly what happens to a polydisperse 
granular fluid at very high packing fractions. As 
for colloids, the timescales for relaxation or diffusion 
increase very fast when density is increased, without 
any noticeable change in structural properties. Again 
it is now established (jD'Anna and Gremaud. 2001 



Marty and Dauchot. 200 



Keys et al 2007 



Chaudhuri et al, 20071 : ICandelier et al, 2009bl ) that 
many phenomcnological properties of the granular glass 
transition also occur in granular assemblies. As for col- 
loids, going beyond the mere analogy and understanding 
how much these different physica l systems are related is 
a very active domain of research ( Liu and Nagel. 200 1[) . 

This very quest ion has been asked in a visual manner 
by Liu and Nagel ( Liu and Nagel. 19981 ) who rephrased it 
in a single picture, now known as a 'jamming phase dia- 
gram'. By building a common phase diagram for glasses, 
colloids and grains, they asked whether the glass and 
jamming transitions of molecular liquids, colloids and 
granular media are different facets of the same 'jammed' 
phase. In this unifying phase diagram, a jammed 'phase' 
(or jammed phases) can be reached either by lowering 
the temperature as in molecular liquids, or increasing 
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the packing fraction or decreasing the external drive in 
colloids and granular media. 



2. More 'glasses' in physics and beyond 

There are many other physi cal contexts in which glassi- 
ness plays an important role (lYoung. 19981) . One of the 
most famous examples is the field of spin glasses. Exper- 
imentally, spin glasses are composed of magnetic impu- 
rities interacting by quenched random couplings. At low 
temperatures, their dynamics become extremely slow and 
they freeze in an amorphous spin configuration dubbed 
a 'spin glass' by P. W. Anderson. Experiments on spin 
glasses, in particular aging studies, have played an impor- 
tant role in the context of amorphous materials. There 
are many other physical systems, often characterized 
by quenched disorder, that show glassy behaviour, like 
Coulomb glasses, Bose glasses, frustrated magnets, etc. 
In many cases, however, one finds quite a different physics 
from structural glasses: the similarity between these sys- 
tems is therefore only superficial from the phenomeno- 
logical point of view, but the theoretical techniques and 
ideas developed in particular in the field of spin glasses 
are highly relevant in thoeretical studies of the glass tran- 
sition. 

Finally, and quite remarkabl y, glassiness emerges even 
in other branches of science ( Mezard et al.. 20071 ). In 
particular, it has been discovered recently that concepts 
and techniques developed for glassy systems turn out to 
apply and be very useful tools in the field of neural net- 
works and computer science. Problems like combinato- 
rial optimization display phenomena completely analo- 
gous to phase transitions, actually, to glassy phase tran- 
sitions. A posteriori, this is quite natural, because a typ- 
ical optimization problem consists in finding a solution 
in a presence of a large number of constraints. This can 
be defined, for instance, as a set of N Boolean variables 
that satisfies M constraints. For N and M very large 
at fixed a = M/N , this problem very much resembles 
finding a ground state in a statistical mechanics prob- 
lem with quenched disorder. Indeed one can define an 
energy function (a Hamiltonian) as the number of un- 
satisfied constraints, that has to be minimized, as in a 
T = statistical mechanics problem. The connection 
with glassy systems lies in the fact that in both cases the 
energy landscape is extremely complicated, full of min- 
ima and saddles. The fraction of constraints per degree 
of freedom, a, plays a role similar to the density in a hard 
sphere system. A detailed presentation of the relation- 
ship between optimization problems and glassy systems 
is clearly out of the scope of the present review. We sim- 
ply illustrate it pointing out that a central problem in 
optimization, random fc-satisfiability, has been shown to 
undergo a glass transition when a increases that is remi- 
niscent of the one of structural g lasses and can be trea ted 
analytically using similar tools (Krzakal a et al.. 2007? ). 
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FIG. 5 Mean-squared displacements of individual particles in 
a simple model of a glass- forming liquid composed of Lennard- 
Jones particles observed on a wide time window. When tem- 
perature decreases (from left to right), the particle displace- 
ments become increasingly slow with several distinct time 
regimes corresponding, in this order, to ballistic, localized, 
and diffusive regimes. 



C. Computer simulations of molecular glass-formers 

Studying the glass transition of molecular liquids at a 
microscopic level is in principle straightforward since one 
has to answer a very simple question: how do particles 
move in a liquid close to T g ? It is of course a daunting 
task to attempt answering this question experimentally 
because one should then resolve the dynamics of single 
molecules to be able to follow the trajectories of objects 
that are a few Angstroms large on timescales of tens or 
hundreds of seconds, which sounds like eternity when 
compared to typical molecular dynamics usually lying 
in the picosecond regime. In recent years, such direct 
experimental investigations have been attempted using 
time and space r esolved techniques such as ato mic 
force microscopy ( Vidal Russel and Israeloff. 2000D or 
single molecule spec troscopy (jAdhikari et al.. 2007t 



iMackowiak et al. 20091 ) , but this remains a very difficult 
task. 

In numerical simulations, by contrast, the trajectory 
of each particle in the system can, by construction, be 
followed at all times. This allows one to quantify easily 
single particle dynamics, as proved in Fig. [5] where the 
averaged mean-squared displacement A(t) measured in a 
simple Lennard- Jones glass-former is shown. The mean- 
squared displacement is defined by 



A(i) 



' 1 N 



ri(0)| s 



(7) 



where (t) represents the position of particle i at time t 
in a system composed of TV particles; the brackets indi- 
cate an ensemble average over initial conditions weighted 
with the Boltzmann distribution. The main observa- 
tion from the data shown in Fig. [5] is that single par- 
ticle displacements considerably slow down when T is 
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decreased. This can be quantified by measuring the 
self-diffusion constant, D s , formally defined as D s = 
lirrif-voo A(f)/(6f). The data in Fig. [5] show that D s de- 
creases by orders of magnitude when temperature de- 
creases, and thus mirrors the behaviour of the (inverse 
of the) viscosity shown in Fig. [2] for real systems. There- 
fore, to explain the phenomenon of the glass transition, 
one must equivalently explain why molecular motions be- 
come so slow at low temperatures. 

Additionally, a very rich dynamics is observed in Fig [51 
with a plateau regime at intermediate timescales, corre- 
sponding to an extended time window during which par- 
ticles vibrate around their initial positions, as in a crys- 
talline solid. The difference with a crystal is of course 
that this transient localization does not correspond to a 
well-defined position in an ordered structure, and it is 
only transient so that all particles eventually escape and, 
concomitantly, the structure relaxes at long times. De- 
scribing the molecular motions responsible for this very 
broad spectrum of relaxation timescales is a challenge in 
this field. 

In recent years, computer experiments have played 
an i ncreasingly imp ortant role in glass transition stud- 
ies ( Andersen. 20051 ). It could almost be said that par- 
ticle trajectories in numerical work have been studied 
under so many different angles that probably very little 
remains to be learnt from such studies in the regime that 
is currently accessible using present day computers. Un- 
fortunately, this does not imply complete knowledge of 
the physics of supercooled liquids. As shown in Fig. [SJ 
it is presently possible to follow the dynamics of a sim- 
ple glass-forming liquid over more than eight decades of 
time, and over a temperature window in which average 
relaxation timescales increase by more than five decades. 
This might sound impressive, but a quick look at Fig. [2] 
shows, however, that at the lowest temperatures studied 
in the computer, the relaxation timescales arc still orders 
of magnitude faster than in experiments performed close 
to the glass transition temperature. Simulations can be 
directly compared to experiments performed in this high 
temperature regime, but this also implies that simula- 
tions focus on a relaxation regime that is about eight 
to ten decades of times faster than in experiments per- 
formed close to T g . Whether numerical works arc useful 
to understand the glass transition itself at all is therefore 
an open, widely debated, question. We believe that it is 
now possible to numerically access temperatures which 
are low enough that many features associated to the glass 
transition physics can be observed: strong decoupling 
phenomena (see Sec. Illlj) . clear deviations from fits to the 
mode-coupling theory (which are experimentally known 
to hold only at high temperatures, see Sec. IIV.B.2[) . and 
crossovers towards truly activated dynamics. 

Classical computer simulations of supercooled liquids 
usually proceed by solving a cleverly discretized version 
of Hamilton's equations for the particles' positions and 
momenta and a given potential interaction between par- 



ticles (|Allen and Tildeslev. 19871) : 



dYi dH dpi dH 
dt dpi ' dt dri ' 



where 



N 



(8) 



(9) 



is the system's Hamiltonian composed of a kinetic part 
and an interaction term V({rj}). We have written 
Eqs. (|HJl and ((9J in terms of the center of mass trajec- 
tories, as is appropriate for atoms although, of course, 
numerical simulation s can deal with molecula r degrees 
of freedom as well (| Allen and Tildeslev. 1987h . bmce 
the equations of motion are energy conserving, they de- 
scribe the dynamics of atomistic systems in the micro- 
canonical ensemble. Constant temperature or constant 
pressure schemes have been developed allowing simula- 
tions to_j2e_rjerfOT7ned_jn_any desired statistical ensem- 
ble ( Allen and Tildeslev. 1987T ). Similarly, noncquilib- 
rium simulation techniques exist that allow, for instance, 
computer studies of the aging dynamics or the nonlinea r 
rheology of supercooled fluids (jEvans and Morris. 20081 ) , 
see also Sec. |V] 

If quantitative agreement with experimental data on 
an existing specific material is sought, the interaction 
must be carefully chosen in order to reproduce reality, 
for instance by combining classical to ab-initio simula- 
tions. From the more fundamental perspective adopted 
here, one rather seeks the simplest model that is still 
able to reproduce qualitatively the phenomenology of 
real glass-formers, while being considerably simpler to 
study. The implicit, but quite strong, hypothesis is that 
molecular details are not needed to explain the behaviour 
of supercooled liquids, so that the glass transition is in- 
deed a topic for statistical mechanics, with little influence 
from chemical details. A considerable amount of work 
has therefore been dedicated to studying models where 
point particles interact via a simple pair potential such 
as Lennard- Jones interactions: 
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where r,j = jr,— Vj |, e and a represent an energy scale and 
the particle diameter, respectively. Other popular models 
are soft spheres, where only the steep short range repul- 
sion in Eq. (|10[) is considered, or even hard spheres where 
the repulsion is made infinitely steep. If the system is too 
simple, such as the one defined in (flU|) . the glass transi- 
tion cannot be studied because crystallization takes place 
when temperature is lowered. Some frustration must be 
introduced to devise numerical models with good glass- 
forming abilities. A common solution, inspired by exper- 
imental studies of metallic glasses, is to use mixtures of 
different atoms, as in the popular model devised by Kob 
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FIG. 6 Intermediate scattering function at wavevector q = 
1.7 A" 1 for the Si atoms from T = 6100 to T — 2750 K, 
obtained from molecular dynamics simu lations of silica. Data 
adapted from (|Horbach and Kob. 200 ll ). 



and Andersen (|Kob and Andersen. 1994h which uses a 
non-additive binary mixture of Lennard- Jones particles. 

More realistic materials arc also studied focusing for in- 
stance on the physics of network forming materials (such 
as silica, Si02, the main component of window glasses), 
multi-component ones, anisotropic particles, or molecules 
with internal degrees of freedom. Connections to ex- 
perimental work can be made by computing quantities 
that are experimentally accessible such as the interme- 
diate scattering function, Eq. (|6]), static structure fac- 
tors, S(q), Eq. ((4]), or thermodynamic quantities such 
as specific heat or configurational entropy, which are di- 
rectly obtained from particle trajectories and can be mea- 
sured in experiments as well. As an example we show in 
Fig. [S] the intermediate scattering function F(q,t) ob- 
tained from a molecular dynamics simulation of a classi- 
cal model for Si02 as a function of t ime for different tem- 
peratures ( Horbach and Kob. 200ll) . The numerical data 
compare very favourably with the experimental results 
obtained from neutron scattering shown in Fig. [3J Note 
that they actually access the dynamics over a broader 
range of timescales and temperatures. 

An important role is played by simulations also because 
a large variety of dynamic and static quantities can be si- 
multaneously measured in a single model system. As we 
shall discuss later, there exist scores of different theoret- 
ical approaches to describe the physics of glass-formers, 
and they sometimes have their own set of predictions 
that can be readily tested by numerical work. For in- 
stance, quite a large number of numerical papers have 
been dedicated to testing in detail the predictions formu- 
lated by the mode-coup ling theory of the glass tra nsition, 
as reviewed recently in ( Gotze. 19 99; Got ze. 20081) . Here, 
computer simulations are particularly well-suited as the 
theory specifically addresses the relatively high temper- 
ature window that is studied in computer simulations. 

While Newtonian dynamics is mainly used in numer- 
ical work on supercooled liquids, a most appropriate 



choice for atomistic materials, it can be interesting to 
consider alternative dynamics that are not determinis- 
tic, or which do not conserve the energy. In colloidal 
glasses and physical gels, for instance, particles undergo 
Brownian motion arising from collisions with molecules 
in the solvent, and a stochastic dynamics is more appro- 
priate. Theoretical considerations might also suggest the 
study of different sorts of dynamics for a given interac- 
tion between particles, for instance, to assess the role of 
conservation laws and structural information. Of course, 
if a given dynamics satisfies detailed balance with re- 
spect to the Boltzmann distribution, all structural quan- 
tities remain unchanged, but the resulting dynamical be- 
haviour might be very different. More generally one can 
ask the question: how universal is the glass transition 
phenomenon? Docs it depend upon the specific micro- 
scopic dynamics? 

Several papers have studied in detail the influence of 
the dynamics on the resulting dynamical behaviour in 
glass-formers using different types of microscopic dynam- 
ics. Gleim et al. studied 'stochastic dynamics' which gen- 
eralizes Newton2s_e^iiations to include non-deterministic 
forces (|GIeim et al. 19981 ): 



1 dt 2 



dr. 
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^ dt 



(ii) 



Specifically, a friction term proportional to the veloc- 
ity with a damping constant £ is added to the right 
hand side, as well as a Gaussian distributed white 
noise rji, whose correlations are related to the damping 
via the fluctuation-dissipation theorem, (Vii^Vji^)) = 
6kBT£5ij6(t — t'), so that the equilibrium distribu- 
tion at temperature T is indeed recovered. When £ 
gets large the dynamics becomes similar to a purely 
Brownian dynamics, as rec ently studied for instance in 
(|Szamel and Flenner. 2"oo3 ): 



^ dt 



dV({rj}) 



(12) 



In that case, particles are described by their positions 
only, and momenta play therefore no role. A simi- 
lar ty pe of description, althou g h numerically m ore effi- 
cient (jBerthier and Kob. 20071: iBerthier. 2007al) . is ob- 
tained using a standard Monte Carlo approach where 
the change in potential energy between two configura- 
tions is used to accept or reject a trial move. In both 
cases of Brownian and Monte Carlo dynamics, particles 
have diffusive (rather than ballistic) behaviour at very 
short times where differences between the different types 
of dynamics can therefore be expected. 

Quite surprisingly, however, the equivalence between 
these three types of stochastic dynamics and the orig- 
inally studied Newtonian dynamics was quantitatively 
established at the level of the averaged dynamical cor- 
relators for all three types of dynamics mentioned above, 
except at very short times where obvious differences 
are indeed expected. This equivalence can probably 
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FIG. 7 Spatial map of single particle displacements in the 
simulation of a binary mixture of Lennard- Jones mixture in 
two dimensions. Arrows show the displacement of each parti- 
cle in a trajectory of length comparable to the structural re- 
laxation time. The map reveals the existence of particles with 
different mobilities during relaxation, but also the existence 
of spatial correlations between these dynamic fluctuations. 
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FIG. 8 Time resolved squared displacements of individual 
particles in a simple model of a glass-forming liquid composed 
of Lennard- Jones particles near the fitted mode-coupling crit- 
ical temperature. The average is shown as a smooth full 
line. Trajectories are composed of long periods of time dur- 
ing which particles vibrate around well-defined positions, sep- 
arated by rapid jumps that are widely distributed in time, 
underlying the importance of dynamic fluctuations. 



III. DYNAMIC HETEROGENEITY 

A. Existence of spatio-temporal dynamic fluctuations 



be traced back to the existence of fast and slow de- 
grees of freedom. It is reasonable to think that the 
former act as an effective thermal bath for the latter 
thus making the three types of dynamics equivalent on 
long timescales. However, this interpretation has to be 
taken with a grain of salt since important differences were 
found whe n fluctuations of dyna mical correlators were 
considered ( Berthier et ai, 2007al) . even in the long-time 
regime corresponding to the structural relaxation. 

Another crucial advantage of molecular simulations is 
illustrated in Fig. [7] This figure shows a spatial map of 
single particle displacements recorded during the simu- 
lation of a binary Lennard- Jones mixture in two dimen- 
sions. This type of measurement, out of reach of most ex- 
perimental techniques that study the liquid state, reveals 
that dynamics might be very different from one particle 
to another. More importantly, Fig.[7]also unambiguously 
reveals the existence of spatial correlations between these 
dynamic fluctuations. The presence of non-trivial spatio- 
temporal fluctuations i n supercooled liquids is called 'dy- 
namic heterogeneity' ( Ediger. 2000[ ). The phenomenon 
has become a substantial component of the field of the 
glass transition, and computer simulations have naturally 
played an important role since they reveal the heteroge- 
neous nature of the relaxation much more directly than 
experiments. We discuss the phenomenon of dynamic 
heterogeneity in more detail in the next section. 



A new facet of the relaxational behaviour of super- 
cooled liquids has emerged in the last decade thanks 
to a considerable experimental and theoretical effort. 
It is called dynamic heterogeneity, and plays now 
a ce ntral role in modern descriptions of glassy liq- 
uids ( Ediger. 200(i ). As anticipated in the discussion 
of Fig. [7] in the previous section, the phenomenon 
of dynamic heterogeneity is related to the spatio- 
temporal fluctuations of the dynamics. Initial motiva- 
tions stemmed from the search for an explanation of 
the non-exponentiality of relaxation processes in super- 
cooled liquids, related to the existence of a broad relax- 
ation spectrum. Two natural, but fundamentally differ- 
ent, explanations can be put forward. (1) The relax- 
ation is locally exponential, but the typical relaxation 
timescale varies spatially. Hence, global correlation or 
response functions become non-exponential upon spa- 
tial averaging over this spatial distribution of relaxation 
times. (2) The relaxation is complicated and inherently 
non-expon ential, even loc ally. Experimental and theoret- 
ical works (lEdieer. 2000Tl suggest that both mechanisms 
arc likely at play, but definitely conclude that relaxation 
is spatially heterogeneous, with regions that are faster 
and slower than the average. Since supercooled liquids 
are ergodic systems, a slow region will eventually become 
fast, and vice- versa. A physical characterization of dy- 
namic heterogeneity entails the determination of the typ- 
ical lifetime of the heterogeneities, as well as their typical 
lcngthscale. 

A clearer and more direct confirmation of the hetero- 
geneous character of the dynamics also stems from sim- 
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ulation studies. For example, whereas the simulated av- 
erage mean-squared displacements are smooth functions 
of time (see Fig. [5]) , time signals for individual particles 
clearly exhibit specific features that are not observed un- 
less dynamics is resolved both in space and time. These 
features are displayed in Fig. [8j What do we see? We 
mainly observe that particle trajectories are not smooth 
but rather composed of a succession of long periods of 
time where particles simply vibrate around well-defined 
locations, separated by rapid 'jumps'. Vibrations were 
previously inferred from the plateau observed at interme- 
diate times in the mean-squared displacements of Fig. [SJ 
but the existence of jumps that are clearly statistically 
widely distributed in time cannot be guessed from av- 
eraged quantities only. The fluctuations in Fig. [8J sug- 
gest, and direct measurements confirm, the importance 
played by fluctuations around the averaged dynamical 
behaviour. 

A simple type of such fluctuations has been studied in 
much detail. When looking at Fig. [8j it is indeed natural 
to ask, for any given time, what is the distribution of 
particle displacements. This is quantified by the self-part 
of the van-Hove function defined as 
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FIG. 9 Decoupling between viscosity (full line) and 
self-diffusi on coefficient (sym bols) in supercooled ortho- 
terphenyl (|Mapes et al.. 20061 ) . The dashed line shows a fit 
with a 'fractional' Stokes-Einstein relation, D 3 ~ (T/r)p with 
£ ~ 0.82 instead of the 'normal' value £ = 1 which only holds 
at high temperatures. 



C.(M) = (^E«('-[ri(t)-r<(0)]) 
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For an isotropic Gaussian diffusive process, one gets 
G s (r,t) = exp(-\r\ 2 /{4D 8 t))/(4irD s t) 3 / 2 . Simulations 
reveal instead strong deviations from Gaussian be- 
havio ur on the timesca les relevant for structural relax- 
ation ( Kob et al.. 19 97). In particular they reveal tails 



in the distributions that are 'fat', in the sense that they 
are much wider than expected from the Gaussian ap- 
proximation. These tails are in fact well described by 
an exponential, rather than Gaussian, decay in a wide 
time window comprising the st ructural relaxation, suc h 
that G s (r, t) ~ exp (— |r|/A(*)) (IChaudhuri et al.. 2007|) . 
Thus, they reflect the existence of a population of parti- 
cles that moves distinctively further than the rest and 
appears therefore to be much more mobile. The ex- 
ponential form of the tail originates from the inter- 
mittent nature of the particle trajectories exposed in 
Fig. |8j m ade of a succession of 'j umps' separated by vi- 
brations (jChaudhuri et al.. 20071) . Actually, such a tail 
would be present in simple j ump models for diffusion 
( Hansen and McDonald. 2006f ). This observation implies 
that relaxation in a viscous liquid differs qualitatively 
from that of a normal liquid where diffusion is close to 
Gaussian, and that a non-trivial statistics of single par- 
ticle displacements exists. 

A long series of questions immediately follows this 
seemingly simple observation. Answering them has been 
the main occupation of many workers in this field over 
the last decade. What are the particles corresponding 
to the tails effectively doing? Why are they faster than 
the rest? Are they located randomly in space or do they 
cluster? What is the geometry, time and temperature 



evolution of the clusters? Are these spatial fluctuations 
correlated to geometric or thermodynamic properties of 
the liquids? Do similar correlations occur in all glassy 
materials? Can one predict these fluctuations theoreti- 
cally? Can one understand glassy phenomenology using 
fluctuation-based arguments? How can these fluctuations 
be detected experimentally? 

Another influential phenomenon that was related 
early on to the existence of dynamic heterogeneity 
is the decoupling of self-diffusion (D s ) and viscos- 
ity (77). In the high temperature liquid self-diffusion 
and viscosity are related by the Stokes-Einstein rela- 
tion ([Hansen and McDonald. 2006|) . D s r\jT = const. 
For a large particle moving in a fluid the constant is equal 
to 1/ (67ri?) where R is the particle radius. Physically, the 
Stokes-Einstein relation means that two different mea- 
sures of the relaxation time R 2 /D s and nR 3 /T lead to 
the same timcscalc up to a constant factor. In super- 
cooled liquids this phcnomcnological law breaks down, a s 
shown in Fig. [9] for ortho-terphenyl ( Mapes et al.. 2006h . 
It is commonly found that Dj 1 does not increase as fast 
as rj/T so that, at T g , the product D s n/T has increased 
by 2-3 orders of magnitude as compared to its Stokes- 
Einstein value. This phenomenon, although less spectac- 
ular than the overall change of viscosity, is a significant 
indication that different ways to measure relaxation times 
lead to different answers and, thus, is a strong hint of the 
existence of broad 'distributions' of relaxation timescales. 

Indeed, a natural explanation of this effect is that dif- 
ferent observables probe di fferently the un derlying distri- 
bution of relaxation times ( Ediger. 2000l ). For example, 
the self-diffusion coefficient of tracer particles is dom- 
inated by the more mobile particles whereas the vis- 
cosity or other measures of structural relaxation probe 
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the timescale needed for every particle to move. An 
unrealistic but instructive example is a model where 
there is a small, non-percolative subset of particles that 
are blocked forever, coexisting with a majority of mo- 
bile particles. In this case, the structure never relaxes 
but the self-diffusion coefficient is non-zero because of 
the mobile particles. Of course, in reality all particles 
move, eventually, but this shows how different observ- 
ables are likely to probe different moments of the distri- 



theoretical frameworks I 


Stillinger and Hodgdon. 1994 


Tarius and Kivelson. 1995 


;|jung et al.. 2004ri. 



The phenomena described above, clearly show that the 
dynamics is spatially heterogeneous. However, they are 
in principle not able to probe whether this is related to 
purely local fluctuations or if there are instead increas- 
ingly spatially correlated fluctuations. This is, however, 
a fundamental issue both from the experimental and the- 
oretical points of view. How large are the regions that 
are faster or slower than the average? How does their 
size depend on temperature? Are these regions com- 
pact or fractal? These important questions were first 
addres sed in pioneering works using four-dime nsional 
NMR (jTracht et al, 19981: iReinsberg et al, 200ll ). or by 



directly probing fluctuations at the nanoscopic scale us- 
ing microscopy techniques. In particular, Vidal Rus- 
sel and Israeloff using Atomic Force M icroscopy tech- 
niques ( Vidal Russel and Israeloff. 200fj ) measured the 



polarization fluctuations in a volume of size of few tens 
of nanometers in a supercooled polymeric liquid (PVAc) 
close to T g . In this spatially resolved measurement, the 
hope is to probe a small enough number of dynamically 
correlated regions, and detect their dynamics. Indeed , 
the signals shown in ( Vidal Russel and Israeloff. 2000h 
reveal that the dynamics is very intermittent in time: 
it switches between moments with intense activity, and 
moments with no dynamics at all, suggesting that ex- 
tended regions of space indeed transiently behave as fast 
and slow regions. A much smoother signal would have 
been measured if such dynamically correlated 'domains' 
were not present. Spatially resolved studies such as NMR 
experiments or atomic force microscopy are quite dif- 
ficult to perform. They give undisputed information 
about the typical lifetime of the dynamic heterogeneity, 
but their determination of a dynamic correlation length- 
scale is rather indirect and performed on a small num- 
ber of liquids in a small temperature window. Never- 
theless, the outcome is that a non-trivial dynamic cor- 
relation length emerges at the glass transition, where it 
reaches a value of the order of about ~ 10 molecule di- 
ameters ( Ediger. 2000f ). 



B. Multi-point correlation functions 

Recently, substantial progress in characterizing spatio- 
temporal dynamical fluctuations was obtained. In par- 
ticular, it is now understood that dynamical fluctuations 



can be measured and characterized through the use of 
high-order correlation and response functions. These 
multi-point functions can be seen as generalizations of the 
spin glass susceptibility measuring the extent of amor- 
phous long-range order in spin glasses. We now intro- 
duce these dynamical functions and summarize the main 
results obtained from their study. 



1. Why four-point correlators? The spin glass case 

No simple static correlation has yet been found to re- 
veal any notable feature upon approaching the glass tran- 
sition (see Sec. IVII for recent theoretical progress). As a 
consequence, it is quite natural to investigate whether a 
growing lengthscale associated with the slowing down of 
the system is hidden in some dynamic correlation func- 
tion. 

Spin glass theory faced a similar conundrum, whose 
solution we briefly recall because it has been instru- 
mental in the developments for glass-forming liquids. 
We know that some hidden long-range order develop s 
at the spin glass transition ( Binder and Young. 19861) . 
However, also for spin glasses, conventional two-point 
functions are useless. Even if spins s x and s x + y have 



non-zero static correlations (s x s x + y ) in the spin glass 
phase, the average over space, [• • • ] x = J dx ■ ■ ■ , 
for a given distance |y| vanishes because the pair- 
wise correlations randomly change sign whenever x 
changes. The insight of Edwards and Anderson is that 
one sh ould first square (s x s x + v ) be fore averaging over 
space ( Edwards and Anderson. 1975f ). In this case, the 
resulting four-spin correlation function indeed develops 
long-range tails in the spin glass phase. This cor- 
relation in fact decays so slowly that its volume in- 
tegral, related to the non-linear magnetic susceptibil- 
ity of the material, diverges i n the whole spin glass 
phase ( Binder and Young. 1986h . 

The Edwards- Anderson idea can in fact be understood 
from a dynamical point of view, which is important for 
understanding both the physics of the spin glass just 
above the transition, and its generalization to structural 
glasses. Consider, in the language of spins, the following 
four-point correlation function: 

G 4 (y, t) = [(s x (t = 0)s x+y (t = 0)s x (t)« x+y (t)}] x . (14) 

Suppose that spins s x and s x+y develop static corre- 
lations (s x s x+y ) within the glass phase. In this case, 
G<4(y, t — > oo) will clearly converge to the spin glass cor- 
relation [(s x s x+y ) 2 ] x . More generally, Gi{y,t) for finite 
t is able to detect transient tendencies to spin glass or- 
der, for example slightly above the spin glass transition 
temperature T c . Close to the spin glass transition, both 
the persistence time and the dynamic length diverge in 
a critical way: 



G 4 (y,i)«j/ 2 - d -"G(|,^, 



(15) 
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where £ ~ (T - T c )~ v and r ~ (T - r c )~ zl/ . As men- 
tioned above, the static non-linear susceptibility diverges 
as / dyG\{y,t — > oo) ~ £ 2 ~''. More generally, one can 
define a time-dependent dynamic susceptibility as: 



(16) 



Xi{t) = / dy G 4 (y,i), 



which defines, provided G4(0, i) is a number of order 1, a 
correlation volume, i.e. the typical number of spins corre- 
lated in dynamic events taking place over the time scale t. 
As we shall discuss below, Xi(t) can a l so be interpreted 
as a quantitative measure of the dynamic fluctuations. 
Note however that the precise relation between \4 an d £ 
depends on the value of the exponent 77, which is physi- 
cally controlled by the detailed spatial structure encoded 
in G4: 



Xi(t = r) oc £ 



2-r, 



(17) 



Therefore, spin glasses offer a precise example of a sys- 
tem which gets slower and slower upon approaching T c 
but without any detectable long-range order appearing 
in two-point correlation functions. Only more compli- 
cated four-point functions are sensitive to the genuine 
amorphous long-range order that sets in at T c and give 
non-trivial information even above T c . In the case of spin 
glasses it is well established that the transition is related 
to the emergence of a low temperature spin glass phase. 

In the case of the glass transition of viscous liquids 
the situation is much less clear. First, unlike spin 
glasses where the disorder is quenched, glass-formers 
tend to freeze in an amorphous state where disorder is 
instead self-induced. Second, there might be no true 
phase transition toward a low temperature amorphous 
phase. It is nevertheless still reasonable to expect that 
the dramatic increase of the relaxation time is due to 
a transient amorphous order that sets in and whose 
range grows approaching the glass transition. Grow- 
ing timescal es should be somehow related to growing 
lengthscales ( Montanari and Semerjian. 2006b| ). A good 
candidate to unveil the existence of this phenomenon is 
the function G 4 (y, i) introduced previously, since noth- 
ing in the above arguments was specific to systems with 
quenched disorder. The only difference is that although 
transient order is detected in G 4 (y, t) or its volume inte- 
gral Xi(t) f° r times of the order of the relaxation time, 
in the long-time limit these two functions may not, and 
indeed do not in the case of supercooled liquids, show 
long-range amorphous order. This roots back to the dif- 
ferent nature of the glass and spin glass transitions. 



2. Four-point functions in supercooled liquids 

In the case of liquids, we may consider a certain space 
dependent observable o(x, t) , such as, for example, the 
local excess density 6p(x.,t) = p(x, t) ~ po, where po is 



the average density of the liquid, or the local dipole mo- 
ment, the excess energy, etc. We will assume in the fol- 
lowing that the average of o(x, t) is equal to zero, and 
the variance of o(x, t) normalized to unity. The dynamic 
two-point correlation is defined as: 



G (r, t) = [o(x, t = 0)o(x + r, t)] 



(18) 



where the normalization ensures that G D (r = 
0,t = 0) = 1. The Fourier transform of C D (r,t) 
defines a generalized dynamic s tructure factor 
S Q (k,t) (jHansen and McDonald. 20061) . All experi- 
mental and numerical results known to date suggest 
that as the glass transition is approached, no spatial 
anomaly of any kind appears in C (r,t) (or in 5 (k, t)) 
although of course there could still be some signal 
which is perhaps too small to be measurable. The 
only remarkable feature is that the slowing down of the 
two-point correlation functions often obeys, to a good 
approximation, "time-temperature superposition" in the 
a- relaxation regime t ~ T a , i.e.: 



C (r,t) ««o(r)/ 



ra(T) 



(19) 



where q is often called the non-ergodicity (or Edwards- 
Anderson) parameter, and the scaling function f(x) de- 
pends only weakly on temperature. This property will be 
used to simplify the following discussions, but it is not a 
needed ingredient. 

The spatial correlations of the relaxation process can 
be probed studying the distribution (over dynamical his- 
tories) of the correlation G (r, t), in particular its covari- 
ance. Quite generally, one expects that since C (r,t) is 
defined as an average over some large volume V, its vari- 
ance T,q is of order ^ 2 ^ TI /V, where £ is the lengthscale 
over which C (r,t) is significantly correlated. More pre- 
cisely we define: 

G 4 (y,t) = [o(x,0)o(x + r,t)o(x + y,0)o(x + y + r,7i)] x 
-[o(x,f = 0)o(x + r,t)]2, (20) 



and its space integral, 



V 



(21) 



which is nothing but the variance of the spontaneous fluc- 
tuations of G D (r, t) averaged over a volume V. This vari- 
ance can thus be expressed as an integral over space of a 
four-point correlation function, which measures the spa- 
tial correlation of the temporal correlation. This integral 
over space is also the Fourier transform of G 4 (y 7 t) with 
respect to y at the wavevector q equal to zero. We want 
to insist at this stage that r and y in the above equa- 
tions play very different roles: the former enters the very 
definition of the correlator we are interested in Eq. (|18l) , 
whereas the latter is associated with the scale over which 
the dynamics is potentially correlated. Correspondingly, 
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great care should be devoted to distinguish the wavevec- 
tor k, conjugate to r, and q conjugate to y. 

Specializing to the case r = (local dynamics), one 
finally obtains: 



Xi (t) EE NY? C . 



(22) 



The analogy with spin glasses developed above sug- 
gests that this quantity reveals the emergence of tran- 
sient amorphous long-range order. Although, as we 
shall discuss, the situation is more complicated than 
what was originally surmised, this analogy was indeed 
the main motivation for the first numerical investigatio n 
of xn(t) in a supercooled liquid (|Dasgupta et al, 199lh . 
It was later realized that Xi{t) is m fact the nat- 
ural diverging susceptibility in the context of Jo- 
spin descriptions of supercooled liquids, where a 
true dynamical pha se transition occurs at a certain 
critical temperature ([Kirkpatrick and Thirumalai. 19881 : 
iFranz and Parish 2000l ). However, since in real systems 
no true phase transition can be observed, one expects 
Xi{t) to grow until t rj r a and decay back to zero there- 
after. Until r ai there cannot be strong differences be- 
tween a system with quenched disorder and a system 
where disorder is dynamically self-induced. 

Measuring the 'local' relaxation suggests to follow the 
displacement of single particles over distances typically 
corresponding to the interparticlc distance. Therefore, 
X4{t) can be accessed cither by measuring the fluctua- 
tions of the Fourier transform of C (r,i) evaluated at 
a wave-vector, fcp, of the order of the first peak in the 
structure factor (jBerthier. 2004 ) , or by performing a spa- 
tial average J drC (r,t)w(r) where w(r) is an overlap 
function equal to o ne for lengths of the order of 2n/ko 
and zero otherwise (jLacevic et al. 20031) or a Gaussia n 
function with a suitable width ( Dauchot et al, 2005f h 
The dependence of dynamical correlations on the 
coarse-graining lengt h has been thoroug hly studied, 
both in simulations ([Chandler et al., 20061) and exper- 
imcnts dDauchot et al, 20051: iDuri and Cipelletti. 20061: 
I Abate and Durian. 2007 ). showing that the dynamics be- 
comes homogeneous when the coarse-graining is made 
too small (where dynamics is dominated by trivial ther- 
mal vibrations), or too large (because dynamics results 
in this limit from a succession of several uncorrelated re- 
arrangements). 

For supercooled liquids, the function X 4(t) 
has been measured by molecular dynamics, 
Brownian and M onte Carlo sim ulations in dif- 
ferent liquids (IFranz et al., 19991: 



Bennemann et al, 19991: 



Parish 1999: 



Lacevic et al, 2003: 



Glotzer. 2000t iBerthier 2004 IVogel and Glotzer. 2004 : 



Berthier. 2007al: iParsaeian and Castillo. 2008 : 

El Masri et al., 20101 ). Moreover, its behaviour has 
been theoretically investigated using various perspec- 
tives, as described below in Sec. IIV1 

As suggested by its definition in Eq. (|2"2")l . the four- 
point dynamic susceptibility is measured in practice 
through the repeated measurement of a chosen local 
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FIG. 10 Time dependence of the self-intermediate scatter- 
ing function (top), and its spontaneous fluctutations (bot- 
tom), for different temperatures decreasing from left to right 
in a Lennard- Jones supercooled liquid in Monte Carlo simu- 
lations. The lowest temperature is highlighted with symbols. 
For each temperature, Xi{t) has a maximum near the relax- 
ation time r a , which shifts to larger times and has a larger 
value when T is decreased, revealing the increasing lengthscale 
of dynamic heterogeneity in supercooled liquids approaching 
the glass transition. Moreover, the time dependence of X i(t) 
is characterized by several distinct time regimes, correspond- 
ing to microscopic, early and late /3 regimes, and structural 
relaxation, as indicated by the vertical lines. 



time-correlation function, Xi(t) being nothing but the 
variance of the statistical fluctuations between different 
measurements. 

A typical example of such a measurement is shown in 
Fig. [TQl taken from Monte Carlo simulations of a sim- 
ple Lennard-Jones supercooled liquid. A similar qualita- 
tive behaviour is found in nearly all cases, as detailed in 
( Toninelli et al, 2005bl ). At a given temperature, Xi(f) 
is an increasing function of time at short times reflect- 
ing the fact that dynamic heterogeneities slowly build 
up with time. It then has a peak on a timescale that 
tracks the structural relaxation timescale r a , and finally 
it decreases back to zero when t — > oo. When temper- 
ature is decreased the observed time evolution becomes 
slower, mimicking the overall slowing down of the dynam- 
ics also seen in averaged two-time correlators, see Fig.lTOl 
The decrease at long times constitutes the above men- 
tioned major difference with spin glasses. In a spin glass, 
Xi{t) would be a monotonically increasing function of 
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time whose long-time limit coincides with the static spin 
glass susceptibility. 

The most important information extracted from the 
temperature evolution of Xi(t) is that, at least in the 
range available to numerical simulations, the value of the 
peak at r Q increases typically from a high temperature 
value of order one, and increases by at most 2 orders of 
magnitude down to the lowest temperature at which the 
system can be equilibrated, suggesting that dynamics be- 
comes spatially increasingly correlated when T decreases. 

As shown in Fig. 1101 the time and temperature be- 
haviour of Xi(t) i s very rich. The growth of Xi(t) towards 
its peak value is composed of several time regimes, closely 
reflecting the broad spectrum of relaxation processes 
characterizing time correlation functions, see Fig. [TO] 
The short-time dynamics corresponding to the approach 
and departure from the plateau in time correlation func- 
tions is also reflected in the time dependence of X4(i), 
which grows with distinct temporal power laws in the 
early and late j3 regimes. Additionally, the temperature 
evolution of X4(£), and in particular the peak height, can 
be quantitatively studied. This peak value, xti measures 
the volume on which the dynamical processes relevant to 
structural relaxation at f w t q , are correlated. It is found 
to increase when the temperature decreases and the dy- 
namics slows down. In the temperature regime above the 
mode-coupling temperature, the growth is well-described 
by an algebraic relation between the peak amplitude, x\i 
and the relaxation time, r a , as shown in Fig. 1101 The 
different exponents intr oduced in Fig. [TU] are discussed in 
(|Toninelli et al, 2005bl) . 

More direct evidences of a growing dynamical cor- 
relation length can be obtained by measuring directly 
Giiy,t). It has been checked that the increase of the 
peak of %4(£) corresponds, as expected, to a growing dy 
n amic lengthscale £ (jYamamoto and Onuki. 1998a ; 



Bennem ann et al, 199 



Doliwa and Heuer. 200 



Lacevic et al, 2003 



Donati et al., 1999: 



Lacevic et al., 2002: 



Berthier et al, 2007al ). 



IWhitelam et al, 2004 ; 
However, these measurements 
are much harder than the ones of X4> because very large 
syste ms need to be simulated to determine £ unambigu- 
ously jStein and Andersen. 20081 : iKarmakar et al., 20091: 
iFlenner and Szamel. 20101) . Note that if the dynamically 
correlated regions were compact, the peak of xa would be 
proportional to £ 3 in three dimensions, directly relating 
X4 measurements to that of the relevant lengthscale of 
dynamic heterogeneity. 

The study of the growth laws of X4(£)> £ an d the evolu- 
tion of G4(y, t) contains very useful information to unveil 
the complexity of the rel axation processes and t o contrast 
theoretical approaches ( Toninelli et al., 2005bl ). In fact, 
many theories of the glass transition assume or predict, 
in one way or another, that the dynamics slows down 
because there are increasingly larger regions over which 
particles have to relax in a correlated or cooperative way, 
see Sec. [TV] 

Furthermore, the growth of £ suggests that the glass 



transition is indeed a collective phenomenon character- 
ized by growing timescales and lengthscales, reminiscent 
of critical phenomena. A clear and conclusive under- 
standing of the relationship between the lengthscale ob- 
tained from G4(y, i) and the relaxation timescale is there- 
fore still the focus of an intense research activity. 



3. Three-point correlation and response functions 

Although readily accessible in numerical simula- 
tions, Xi(t) is m general very small and difficult 
to measure directly in experiments, except when 
the range of the dynamic correlation is macroscopic , 
as in granular materials dMartv and DauchotT 20051: 
iDauchot et al, 2005 : Keys et al, 20071) or in soft glassy 
materials (j Weeks et al, 20071 ). where it can reach 
in some cases the micrometer and even millime- 
ter range dMaver et al. 2004 iDuri and Cipelletti. 200"6t 
M accarrone et al, 20101 ). To access Xi(t) m molecular 
liquids, one should perform time-resolved dynamic mea- 
surements probing very small volumes, with a linear size 
of the order of a few nanometers. Although doable, such 
experiments remain to be performed with the needed ac- 
curacy. 

It was rece ntly realized that simp ler alternative pro- 
cedures exist ( Berthier et al, 200 5a). The central idea 
underpinning these results is that induced dynamic fluc- 
tuations are in general more easily accessible than spon- 
taneous ones, and both types of fluctuations can be re- 
lated to one another by fluctuation-dissipation theorems. 
The physical motivation is that while four-point correla- 
tions offer a direct probe of the dynamic heterogeneities, 
other multi-point correlation functions give very useful 
and direct information about the microscopic mecha- 
nisms leading to these heterogeneities. For example, one 
might expect that the slow part of a local enthalpy (or 
energy, density,...) fluctuation per unit volume, 5h, at 
position x and time t — triggers or eases the dynamics 
in its surroundings, leading to a systematic correlation 
between 5/i(x, t = 0) and o(x + y, t = 0)o(x + y + r, r a ). 
This physical intuition suggests the definition of a family 
of three-point correlation functions that relate thermo- 
dynamic or structural fluctuations to dynamical ones. 
Interestingly, and crucially, some of these three-point 
correlations are both experimentally accessible and give 
bounds or approximations to the four-point dynamic cor- 
relations, as we now detail. 

In the same way that the space integral of the four- 
point correlation function is the variance of the two-point 
correlation, the space integral of the above three-point 
correlation is the covariance of the dynamic correlation 
with energy fluctuations: 



VN J dxdx '°( x ' + r,i)o(x',0)<5/i(x,0) 
1 j dy[ O (x + y + r,t) O (x + y,0)^(x,0)] x . (23) 
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Note that for the enthalpy we use the notation H(t = 
0) = jf J dxh(x., t = 0), so that h is an enthalpy per unit 
volume. Hence, using the fact that the enthalpy fluctua- 
tions per particle are of order ^fcpksT (where cp is the 
specific heat in kp units), the quantity NY*ch / \[cp^bT 
defines the number of particles over which enthalpy and 
dynamics are correlated. Of course, analogous identities 
can be derived for the covariance with, say, energy or 
density fluctuations. 

Although interesting in itself, the covariance T,qh is 
just as hard (or even harder) to measure experimentally 
as X4- However, Y*ch can be related, using linear re- 
sponse theory, to a response function which is much eas- 
ier to access in experiments. Let us prove this by con- 
sidering a system in the grand-canonical NPT ensemble. 
The probability of a given configuration C is given by the 
Boltzmann weight exp(-/3H[C])/Z, where (3 = l/k B T 
and Z is the grand-partition function. Suppose one stud- 
ies a static observable O with the following properties: (i) 
O only depends on the current microscopic configuration 
C of the system and (ii) O can be written as a sum of 
local contributions: 



O = - I dy o(y). 



(24) 



In this case, a well-known static fluctuation -dissipation 
theorem holds ( Hansen and McDonald. 2006[) : 



9(0) 
d/3 



dy (o(y)Sh(O)) = -NE 



OH, 



(25) 



where we decompo sed the enthalpy in a sum of lo cal con- 
tributions as well (jHansen and McDonald. 20061) . 

Interestingly, in the case of deterministic Hamiltonian 
dynamics, the value of any local observable o(x, t) can 
be seen as a highly complicated function of the initial 
configuration at time t = 0. Therefore, a two-time corre- 
lation function, now averaged over both space and initial 
conditions, can be rewritten as a thermodynamical aver- 



age: 



C (r, t; T) = J dx o(x + r, f)o(x, t = 0) 



x exp 



(3 / dx h(x,t = 0) 



(26) 



Hence, the derivative of the correlation with respect to 
temperature (at fixed pressure) directly leads, in the 
case of purely conservative Hamiltonian dynamics, to the 
covariance between initial energy fluctuations and the 
dynamical correlation, in direct analogy with Eq. ()25l) . 
Defining 

G T (y,t) = (o(y + v,t)o(y,0)6h(0,0)}, (27) 



one finds: 

dC (v,t;T) 

XtM) = 



dT 



kpT 2 



dyG T (y,t). (28) 



Hence, the sensitivity of a two-time dynamical func- 
tion to temperature, \t, is directly related to a 
three-point spatial correlation function. The above 
result in Eq. (|28p is extremely general and ap- 
plies to many different situations. However, it does 
not apply when the dynamics is not Newtonian, 
as for instance for Brownian particles or in Monte- 
Carlo numerical simulations ( Doliwa and Heuer. 2000: 
iBerthier and Kob. 20071) . The reason is that in these 
cases, not only the initial probability but also the transi- 
tion probability from the initial to the final configuration 
itself explicitly depends on temperature. In Brownian 
dynamics, for example, the noise in the Langevin equa- 
tion depends on temperature. Hence, dC (r,t;T)/dT 
receives extra contributions from the whole trajectory, 
that depend on the explicit choice of dynamics. 

The equality (|2"5)) . although in a sense a trivial result 
obtained from linear response theory, has a deep physical 
consequence, which is the growth of a dynamical length 
upon cooling in glassy systems, as we show now. Define 
r a (T) such that C o (0,t = r a ;T) = e' 1 (say). Differenti- 
ating this definition with respect to T gives 



_ dT a dC o (0,t = Ta ;T) | dC o (0,t = r a ;T) 



dT 



dt 



dT 



Since C o (0,i;T) decays from 1 to zero over a time scale 
r Q , one finds that generically, using Eq. 



dy 



(o(y,t 



)o(y, 0)^(0,0)) 



T d\nr a 



P y/Cp~k B T 



Pi) 



c P dT 

(30) 

Now, since 5h is of order po^/cpkpT and (o 2 ) is normal- 
ized to unity, the quantity xo = Gt(0, r a )/ po^/cpksT is 
not expected to appreciably exceed unity. The above in- 
tegral can be written as Xo v t, which defines a volume vt 
over which enthalpy fluctuations and dynamics are ap- 
preciably correlated. Note that the interpretation of vt 
as a true correlation volume requires that %o be of order 
one, and its increase is only significant if xo is essentially 
temperature independent. If this is not the case, then the 
integral defined in Eq. (|30[) could grow due to a growing 
Xo and not a growing length, which would obviate the 
notion that Vt is a correlation volume. 

Assuming xo < 1, a divergence of the right hand side 
of the equality (|28|) necessarily requires the growth of vt- 
More precisely, as soon as r a increases faster than an in- 
verse power of temperature, the slowing down of a Hamil- 
tonian system must necessarily be accompanied by the 
growth of a dynamic correlation length, which is there- 
fore a general, powerful consequence of the use of linear 
response theory. This result is thus directly relevant to 
supercooled liquids, where T a typically increases in an ac- 
tivated manner, with, possibly, a finite temperature dy- 
namic singularity. From a theoretical perspective, it also 
implies that any theory predicting a dynamic singularity 
necessarily contains a prediction for diverging dynamic 
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lcngthscales accompanying the glass transition 1 . 

The study of these three-point correlation and re- 
sponse functions is therefore a useful path to character- 
ize dynamical heterogeneity and dynamical correlations. 
Quantitative results can be obtained studying experi- 
mentally xt- This has been done in connection with 
an inequality on xa that we shall describe in the fol- 
lowing section. Another interesting development that 
will be discussed later on consists in focusing on re- 
sponse functions, like xt, but where the perturbing field 
i s spatially depend ent, e.g. with an oscillatory shape 
(jBiroli et al., 20061) . This allows one to probe directly 
the size and the shape of the dynamically correlated re- 
gions. 

Before concluding, let us stress that we have consid- 
ered the response of time correlations to a temperature 
change, but other perturbing fields may also be relevant, 
such as density, pressure, concentration of species in the 
case of mixtures, etc. For example, for hard-sphere col- 
loids, temperature plays very little role whereas small 
changes of de nsity can lead to enormous ch anges in relax- 
ation times ( Pusev and van Megen. 19861 ). The deriva- 
tion of this section can be straight forwardly extended to 
these cases ( Berthier et ai. 2007al) . 



enthalpy fluctuates), one obtains: 



4. Inequalities on xa an d experimental measurements 

In previous subsections, we argued that Xi(t) is a 
fundamental quantity in order to understand dynamic 
heterogeneities in supercooled liquids, but we then pro- 
ceeded to describe a series of alternative multi-point sus- 
ceptibilities, in particular XT{t), which contain alterna- 
tive information on heterogeneities. We now close the 
loop and show that both types of susceptibilities are in 
fact not independent from one another, but closely re- 
lated. 

This can be done using the general fo rmalism devel- 
oped long ago in ( Lebowitz et al.. 1967t ). which gives 
expressions for the strength of fluctuations of physical 
observables measured in distinct statistical ensembles. 
Applied to the spontaneous fluctuations of global two- 
time correlation functions, and considering transforma- 
tion from NPH (where enthalpy is fixed and tempera- 
ture fluctuates) to NPT (where temperature is fixed but 



1 An intriguing case, which is not fully understood, is the exam- 
ple of systems with an Arrhcnius behavior at low temperature. 
The general considerations laid out in the text suggest that these 
systems are characterized by a dynamical correlation length di- 
verging at zero temperature, which contrasts with the idea that 
relaxation in Arrhenius systems is a simple, locally activated pro- 
cess. However, the present results only hold for energy conserv- 
ing systems for which t hermal activation m ay be more collective 
than usually surmised l lTarzia et al. 20201) . 
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where Xi PH (*) is the variance of the correlation function 
in the NPH ensemble where enthalpy does not fluctuate, 
a manifestly non-negative quantity. This allows us to 
obtain the following inequality: 

T 2 X 2 T (t) _ T 2 fdC o (0,t;T) 
cp cp 



dT 



(32) 



Note that there is a simpler way to obtain the above 
inequality. In the previous section, xrif) was shown to 
be related to the covariance of enthalpy and dynamic 
fluctuations, Y>ch ■ Since X4 is related to the covariance 
of dynamic fluctuations T*cc, one can easily check that 
Eq. (|32p is just a rewriting of the Cauchy-Schwarz bound: 



< 



Y? C Y? H , where Tl 2 H is the variance of the enthalpy 



fluctuations, equal to c P (k B T) 2 /N in the NPT ensem- 
ble. 

The inequality (f3"2"j) is very interesting because the 
right hand side is an experimentally measurable quan- 
tity which therefore provides a direct lower bound on Xi- 
Thus, if T 2 x\(t) I cp increases substantially above one, 
Xa has to increase as least as much if not more. In par- 
ticular, as soon as xt increases faster than T _1 at low 
temperatures, Xi w iH eventually exceed unity; since Xi 
is the space integral of a quantity bounded from above, 
this again means that the lengthscale over which the 
four-point correlation G&(y,t) extends has to grow as 
the system gets slower and slower. Again, more quan- 
titative statements require information on the amplitude 
and shape of G^y, t), which has to be provided from 
theoretical or numerical calculations. 

Equation (|3ip makes precise the intuition that dy- 
namic fluctuations are partly induced by the fluctua- 
tions of quantities that physically a ffect the dynamic be- 
havior (|Ediger. 19981: iDonth. 200ll) . in that case the en- 
thalpy. The inequality (|32p provides a correct estimate 
of Xi if there are no "hidden" variables which also con- 
tribute to the dynamic fluctuations. It is however quite 
difficult to determine whether such additional contribu- 
tions exist. Theoretical investigations in the context of 
approximate models for the glass transition described in 
Sec. IIV1 and detailed numerical calculations where all 
terms in Eq. (|3"Tj) can be separately evaluated, greatly 
clarified this issue. The central conclusion is that the ex- 
perimentally accessible response functions which quan- 
tify the sensitivity of average correlation functions to an 
infinitesimal change in control parameters can be used 
in Eq. (|32[) not only to yield lower bounds to X<i(*)j but. 
in fact to provide useful quantitative approximations to 
four-point functions. Although the relative precision on 
Xi{t) provided by the bound at a given temperature is 
modest, the rate of growth is accurately reproduced. 

The results discussed above have opened the way to ex- 
perimental characterization of the growth of Xi i n molec- 
ular glass- formers close to T g . In order to make use of the 
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FIG. 11 Dynamic scaling relation between number of dynam- 
ically correlated particles, JVcorr,4i an d relaxation timescale, 
T a , for a number of glass-formers, determined using Eq. (|32|l . 
For all materials, a similar trend is found, with a rapid 
initial increase of AT C orr,4 near the onset of slow dynam- 
ics, followed by a slower, presumably logarithmic beahvior, 
closer to the laboratory glass temperature. Adapted from 
(jDalle-Ferrier et al. 20071 '). 



Figure [IT] also indicates that N co „^ grows faster when 
r Q is not very large, close to the onset of slow dynamics, 
and a power law relationship -/V CO rr,4 ~ i~® (see Fig. [TO]) 
is a good description in this regime (T a /To < 10 4 ). 
The growth of iV corr 4 becomes much slower closer to 
T g . A change of 6 decades in time corresponds to a 
mere increase of a factor about 4 of iV corri 4, suggesting 
logarithmic rather than power law growth of dynamic 
correlations. This is in agreement with several theories 
of the glass transition which are based on activated 
dynamic scaling, see Sec. IIVI 

Note that all the results above can be general- 
ized to cases where the order parameter inducing 
the glass transition is the densi ty. This happens for 
colloids and granular media dBerthier et al. 2005al 
iLechenault et al. 2008at iBrambilla et al. 20091 ). One 
can obtain a new inequality where density fluctuations 
play the same role of enthalpy fluctuations and that 
provide a lower bound to \i m terms of the derivative 
of the correlation function with respect to the den- 
sity. The different inequalities yielding experimentally 
accessible ways to quantify the strength of dynamic 
heterogeneity certainly have the appeal of simplicity, 
and are by now routinely used in many different systems 
(IStevenson and Wolvnes. 20061 : ICapaccioli et al. 2008 : 



Dalle- Ferrier ef al. 200*1 



Koppensteiner et al. 2010): 



Chen et al. 2009: 



Maggi et al. 201fj iRoland et al. 20i0f ) 



Gainaru et al. 201C ; 



inequality (p?2"j) one must be able to detect time correla- 
tion functions at constant temperature, Co(t; T) with suf- 
ficient precision that dynamics at nearby temperatures, 
T and T + AT can be resolved, when AT is small enough 
that linear response holds: 



dC (t;T) C (t;T) - C (t;T + AT) 



dT 



A simpler 
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discussed in 
C (t;T) with 



(jDalle-Ferrier et al 
an empirical form containing a few number of parame- 
ters, and then take the temperature derivative of these 
parameters to indirectly estimate XT(t)- Combining 
these different methods, the lower bound (|32|) to \i 
has been computed for many molecular glass-forming 
liquids. It is easy to convince oneself that the lower 
bound has also the correct time dependence, with a 
peak developing for times of the order of the relaxation 
timescale. The value of this peak is therefore a lower 
bound to the peak of \i and, hence, to the number 
of dynamically correlated molecules, denoted iV corri 4. 
We show in Fig. [TTJ the evolution of iV corri 4 for many 
different glass-formers in the entire supercooled regime 
(jDalle-Ferrier et al. 2007t ) as a function of the relaxation 
timescale r Q . Note that, actually, iV corrj 4 is expected 
to be equal to the number of dynamically correlated 
molecules up to a proportionality constant which is not 
known from experiments, probably explaining why the 
high temperature values of A^ corr .4 are smaller than one. 



C. Current status of dynamic heterogeneity studies 

The present section on dynamic heterogeneity was 
a very brief summary of a collective research effort of 
very large amplitude that lasted about 20 years, and 
which already forms the core of several re cent reviews 
' Sillescu. 199a lEdieer. 20001 : iRichert. 2002h . and a book 



Berthier et al. 201 if ). Progress to characterize, visual- 



ize, and quantify dynamic heterogeneity as well as an ex- 
ploration of its detailed physical consequences have been 
truly dramatic in recent years. The impact of this re- 
search is such that tools developed to study dynamic het- 
erogeneity in liquids are now routinely used in scores of 
different systems, and 'dynamic heterogeneity' is a con- 
cept that is commonly employed in a broad range of sit- 
uations, much beyond the physics of the glass transition 
([Berthier et al. 201lh . 

Despite this progress, several key questions are still 
unanswered. Our discussion above has been focused on 
the issue of the characterization of the spatial fluctua- 
tions involved in the phenomenon of dynamic hetero- 
geneity. This is justified because direct measurements 
of growing dynamic lengthscales have provided the long- 
sought evidence in favor of the collective nature of the 
glass transition itself. This fact being now established, it 
remains to understand more precisely and quantify the 
connection between these growing lengthscales and the 
increasing viscosity of liquids approaching the glass tran- 
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sition, which appears as a topic for research in the coming 
years. 

We have described in some detail the physical con- 
tent of multi-point dynamic susceptibility such as X±(P)- 
These functions have played a major role in the above 
story, but we now understand that they contain also a 
number of embarrassing features. For instance, we men- 
tioned how Xi(t) retains a dependence on the statisti- 
cal ensemble where it is measured, as in Eq. (|31[) . and 
that it also depends on the microscopic equations of mo- 
tion for the system (Newtonian versus Brownian) . These 
subtle issues make the analysis of four-point susceptibil- 
ities somewhat ambiguous, especially when estimates for 
lengthscales are sought. We shall describe below, for in- 
stance in Sec. lIV.B.3[ that alternative dynamic functions 
now exist that should be easier to analyze, but these have 
not all been measured in simulations or experiments yet. 
Thus, more detailed studies of a larger family of dynamic 
susceptibilities are certainly most wanted in the future. 
A very promising avenu e of research consists in stu dy- 
ing non-linear responses (jBouchaud and Biroli. 20051) . A 
first pioneering experimental measurement of non-linear 
dielectr ic susceptibility for glass-form ers appeared re- 
cently (ICrauste-Thibierge et ai. 20101 ). 

Additionally, direct experimental measurements of dy- 
namic lengthscales are still not available for molecular 
glass-formers, and are scarce even for colloidal materials. 
Thus, it would be useful to develop new experimental 
tools to resolve the dynamics of molecular glass-formers 
on small lengthscales and longer timcscalcs to obtain this 
much needed information. It is not yet clear whether 
molecular dynamics simulations of model systems have 
covered a broad enough range of timescales and are thus 
relevant to understand the physics of real glass-formers 
near the experimental glass transition temperature. We 
also believe that further work should be devoted to a bet- 
ter characterization of the geometry (and not only typical 
lengthscale) of the dynamically heterogeneous regions, 
since contradicting results are available in the literature 
dDonati et ai. 19981: lAppignanesi et ai. 20061) . 



IV. SOME MODELS AND THEORETICAL 
APPROACHES 

A. A few key questions 

We now present some theoretical approaches to the 
glass transition. It is impossible to cover all of them in 
this review, simply because there are way too many of 
them. This is perhaps the clearest indication that the 
glass transition remains an open theoretical problem. 

We have chosen to present in some detail those ap- 
proaches that, we believe, contain keystone ideas and 
at the same time have a solid statistical mechanics ba- 
sis. Loosely speaking, they have a Hamiltonian, can 
be simulated numerically, or studied analytically with 
tools from statistical mechanics. Of course, the choice 



of Hamiltonian is crucial and contains very important 
assumptions about the nature of the glass transition. 
All the approaches we present have given rise to unex- 
pected results. One finds more in them than what was 
supposed at the beginning, which leads to new, testable 
predictions. Furthermore, with models that are precise 
enough, one can test (and hopefully falsify!) these ap- 
proaches by working out all their predictions in great de- 
tail, and comparing the outcome to experimental data. 
Such detailed analysis is often not possible with 'phys- 
ical pictures', or simpler phenomenological modeling of 
the problem. Our drastic choice of theories leaves be- 
hind many other approaches that, although interesting, 
could not be covered without increasing the length of 
this review beyond reasonable limits. Recent reviews 
are av ailable on these and we refer the interested read- 
ers to (iDebenedetti. 199f3 iDonth. 20011: IScior tino. 20051: 
iDvre. 200fA IChen et ai. 20101) . 

Before going into the models and theories, we would 
like to formulate a few important questions that theoreti- 
cians seek to address and that will guide our presentation 
of theories below: 

(1) Why do the relaxation time and the viscosity in- 
crease when T g is approached? Why is this dramatic 
growth different from an Arrhenius law? 

(2) Can one understand and describe quantitatively 
the broad relaxation spectra characterizing the dy- 
namical behaviour of supercooled liquids, in particular 
non-exponential relaxations, and their evolution with 
fragility? 

(3) Is there a deep relation between kinetics and ther- 
modynamics (such as To ~ Tk), and why? 

(4) Can one understand and describe quantitatively 
the spatio-temporal fluctuations of the dynamics? How 
and why are these fluctuations related to the dynamic 
slowing down? 

(5) Is the glass transition a collective phenomenon? If 
yes, of which kind? What is the correct 'order parame- 
ter', and the nature of the associated transition? 

(6) Is the experimental glass transition at T g the mani- 
festation of a finite or zero temperature phase transition, 
sometimes called the 'ideal glass transition' ? Or is there 
instead an avoided, hidden, or inaccessible transition? 

(7) Is there a geometric, real space explanation for the 
dynamic slowing down that takes into account molecular 
degrees of freedom? 

(8) Are there simplified (e.g. lattice) glass models 
which essentially capture the physics of the glass tran- 
sition of molecular liquids? 

Before embarking in detailed theoretical explanations 
it is important to stress that the glass transition appears 
as a kind of 'intermediate coupling' problem, since for 
instance typical correlation lengthscales are found to be 
at most a few tens of particle long close to T g . As a con- 
sequence, recognizing and validating 'the' correct theory 
is extremely difficult since key signatures could be buried 
(and probably are) under preasymptotic, microscopic de- 
tails. These are probably useful also to make incor- 



21 



rect theories appear reasonable. To obtain quantitative, 
testable predictions, one must therefore be able to work 
out also these preasymptotic effects. This is a particu- 
larly difficult task, especially in cases where the asymp- 
totic theory itself has not satisfactorily been worked out 
yet. Therefore, at this time, theories are mainly judged 
by their overall predictive power and theoretical consis- 
tency. 



B. Mean-field free energy landscapes and Random First 
Order Transition (RFOT) theory 

1. Mean-field glass theory and complex free energy landscapes 

In the last two decades, th ree independent lines o f 
research, Adam-Gibbs t heory (lAdam and Gibbs. 19651) . 
mode-co upling theory dGotze. 2008 ) and spin glass 
theory ijMezard et ai. 19881 ) . have been merged in 
a common framework to produce a theoretical en- 
semble that now goes under the name of Ran- 
dom First Order Transition theory (RFOT), a ter- 
minology introduced by Kirkpat rick, Thirumalai and 
Wolynes ( Kirkpatrick et ai, 1987t) who played, among 
many other researchers, a key role in its development. 
Here, we do not follow the rambling developments as 
they took place, but summarize RFOT theory in a more 
modern and unified way. Note that our use of the name 
'RFOT' is much broader than the more common, but 
much narrower meaning often implied in the literature. 
Reviews dedicated to di fferent aspects of RFOT theory 
have appeared recently dLubchenko and Wolynes. 2007 ; 
Gotze. 20081 ICavagna. 2009t iBiroli and Bouchaud. 2009 ; 



Mezard and Parisi. 201QT) 

As discussed previously, two hallmarks of the dynam- 
ics of glass-formers are that (i) close to T g a liquid remain 
stuck for a very long time in amorphous configurations 
and (ii) the number of these configurations is exponen- 
tially large in the system size. RFOT theory starts as a 
mean-field approach to these phenomena. As such it has 
to be able to capture the right kind of symmetry breaking 
and deal with an exponential number of states. 

Broadly speaking, mean-field theory is based on the 
study of the free energy landscape as a function of the 
order parameter. For example, for ferromagnets, in the 
Curie- Weiss approach, one computes the free-energy as a 
function of the global magnetization by a mean-field ap- 
proximation. This yields direct access to the nature and 
properties of the ferromagnetic phase transition, and a 
simple description of the low temperature phase since 
the two minima of the free energy energy correspond to 
the two ferromagnetic states. However, computing the 
free energy as a function of the global energy or density 
is not enough for the glass transition, because one must 
deal with the existence of many different amorphous con- 
figurations. As a consequence one is forced to compute 
the free-energy F as a function of the entire density field 
(instead of a single variable as in the ferromagnetic tran- 



sition), F being defined through the Legendre transform. 

Consider for simplicity an interacting particle lat- 
tice model, the generalization to continuum systems is 
straightforward. In the lattice case a given configuration 
is determined by the number of particles, ri£, on each site 
i. In order to define F, one first introduces the thermo- 
dynamic 'potential' 

W({m}) = -ilog^ exp l-(3H({n t })+J2P^nA , 

{fit} \ i J 

where H({rii}) is the Hamiltonian. The free energy func- 
tion F({pi}) is defined as 



(34) 



where the p*s satisfy the equations + p, L = and, 
hence, are functions of {pi}, which specifies an aver- 
aged density profile. Note that this construction can be 
generalized to spin systems replacing the positive inte- 
ger m by ±1 variables s,. In this context F is called 
the 'TAP' free energy since its introduction by Thoulcss, 
Anderso n and Palmer in the context of mean-field spin 
glasses ( Thouless et ai, 19771 ). The generalization to 
continuum system can be also performed by replacing the 
discrete variable n, by a continuum densi ty field p(x). I n 
this case F is called 'density functional' ( Oxtobv. 1990T) . 

The free energy landscape is the hyper-surface gener- 
ated by scanning F over all possible values of {pi}- Its 
critical points, in particular the minima, play a crucial 
role. In fact, by deriving Eq. (|34j) with respect to pi one 
finds 



dF 
dp% 



Mi- 



(35) 



Thus, when there are no external fields (or local chemical 
potentials) the solutions of these equations are all the 
stationary points of the free energy landscape 2 . 

What are the main features of F for a system ap- 
proaching the glass transition? Unfortunately, this ques- 
tion cannot be answered exactly for a realistic three- 
dimensional system. One has either to make use of ap- 
proximations (as in the Curie- Weiss description of ferro- 
magnets) or focus on simplified geometries, such as mean- 
field Bethe lattices, which, hopefully, provide a good ap- 
proximation to finite dimensional ones. 

A quite large number of such studies have led to sim- 
ilar results, and thus to a consistent mean-field picture. 
The free energy landscape becomes 'rugged' at low tem- 
perature and characterized by many minima and saddle 
points. Actually, the number of minima is exponentially 



2 For particle systems there is always a global chemical potential 
(j, fixing the number of particles. In this case, one includes the 
global term fi • n; in the definition of F so that all fi* are zero. 
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large in the system size, which suggests the definition of 
an entropy, called 'configurational entropy' or 'complex- 
ity': 



s c = —logN(f), 



(36) 



where Af(f) is the number of free-energy minima with 
a given free energy density / (per unit of free energy 
density). The density profile corresponding to one given 
minimum is amorphous and lacks any type of periodic 
long-range order, and different mimima are very different. 
This is a very welcome theoretical result, as real glasses 
can be found in a large number of different amorphous 
configurations, which can be interpreted as mean-field 
free energy minima. 

Assuming that all minima are mutually accessible, one 
can compute the thermodynamic properties, i.e. the par- 
tition function by summing over all states with their 
Boltzmann weights: 



Z = e 



-pw 



(37) 



where the sum runs over the minima. Formally, one can 
introduce a free-energy dependent complexity, s c (f,T), 
that counts the number of free-energy minima with free- 
energy density / at temperature T. The partition func- 
tion of the system then reads: 



Z(T) 



df exp 



T 



Ns c (f,T) 



(38) 



For large N, one can as usual perform a saddle-point 
estimate of this integral, that fixes the dominant value of 
/, noted f*(T): 



ds c (f,T) 



df 



1=1* (T) 



1 

T' 



(39) 



The temperature dependent complexity is therefore 
s c (T) ee s c (f*(T),T). The total free energy of the sys- 
tem is f p (T) = f* - Ts c (T). A typical shape of the 
configurational entropy as a function of / and a graphic 
solution of Eq. ([3"5f are plotted in Fig [TJ] The analysis 
of the configurational entropy, or complexity s c (T), re- 
veals that s c (T) decreases when temperature decreases, 
as long as T is above a critical temperature, Tk, below 
which s c (T) vanishes. There exists also a second, higher 
temperature, Tmct > Tk, above which s c (T) drops dis- 
continuously to zero again. We use the symbol Tmct 
(as in mode-coupling theory) on purpose, and justify our 
choice below. Although the complexity vanishes in both 
regimes, the situations below Tk and above Tmct are 
very different. At these two temperatures the part of 
the free-energy landscape relevant for the thermodynam- 
ics changes drastically in two very different ways. Above 
Tmct, there is just one minimum dominating the equilib- 
rium measure corresponding to the homogeneous density 
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FIG. 12 Typical shape of the configurational entropy, s c , as 
a function of free energy density, / in the range Tk < T < 
Tmct for random first order landscapes. A graphic solution 
of Eq. (|39[) is obtained by finding the value of / at which the 
slope of the curve is j3. Note that s c is also a function of 
temperature, so this curve in fact changes with T. 



profile of the high temperature liquid. At Tmct the ho- 
mogeneous liquid state becomes fragmented in an expo- 
nential number of states, or minima. At Tk the number 
of minima becomes sub-exponential in the system size, 
such that s c (T <T K ) = 0. 

Surprisingly the total free energy f p (T) is not sin- 
gular at Tmct- This is one of the most unexpected 
result consistently emerging from analytical solutions. 
This suggests that at Tmct the liquid state fractures 
into an exponential number of amorphous states, but 
that this transition has no thermodynamical counter- 
part, and is therefore a purely dynamical phenomenon. 
At Tk instead, a thermodynamic phase transition takes 
place since the contribution to the entropy coming from 
the configurational entropy disappears, typically linearly, 
s c (T) ~ (T — Tk)- Therefore, the specific heat is found 
to make a sharp downward jump at Tk, thus providing 
an exact realization of the 'e ntropy vanishing' mechanism 
conjectured by Kauzmann ( Kauzmann. 194q ). This is a 
second welcome result: the thermodynamic signature of 
this mean-field transition mirrors the basic experimental 
finding that the specific heat is nearly discontinuous at 
the experimental glass temperature T g . 

This rich physical behavior can be derived from a 
number of perspectives. A first concrete example is 
given by 'lattice glass models' ( Biroli and Mezard. 200l[ ) 
solved by the Bethe approximat io n or on Bethe lat- 
tices dPica Ciamarra et ai. 20031 : iRivoire et ai, 200"3 : 
Iliansen-Goos and Weigt. 2005 ). Lattice glass models 
contain hard particles sitting on the sites of a lattice. 
The Hamiltonian is infinite if there is more than one 
particle on a site and, more crucially, if the number of 
occupied neighbors of an occupied site is larger than 
a fixed parameter, m. The Hamiltonian is zero other- 
wise. Tuning the parameter m, changing the type of 
lattice, in particular its connectivity, can yield differ- 
ent models. Lattice glasses are simple statistical me- 
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chanical models mimicking the physics of hard sphere 
systems. Numerical simulations on cubic lattices have 
shown t hat they seem to behave as bona fide glas s- 
formers rtPica Ciamarra et ai. 20031: iDarst et al., 20101) . 

Alternatively, a density functional theory analy- 
sis of t he free energy la n dscape yields very simila r 
results (|Singh et al. 19851 : IDasgupta and Vails. 2000l ). 
This is a more realistic microscopic starting point, but 
it inevitably contains some approximations, in partic- 
ular related to the specific fo rm of the free energ y 
functional ( Singh et al. 19851 IDasgupta et al.. 199lh . 
The adopted form is the Ramakrishnan-Youssouf den- 
sity functional and most studies focused on hard 
spher e systems. In the first of these investiga- 
tions ( Singh et al.. 1985T ) a particular amorphous profile, 
whose only free parameter was the cage radius over which 
particles are free to vibrate was plugged in the density 
functional. Minimization with respect to the cage radius 
revealed that amorphous structures become stable, in a 
variational sense, at high enough density. More recent in- 
vestigations performed a full minimization and reached 
qualitatively similar, but much more detailed con clusions 
( Dasgupta et al.. 19911 : IChaudhuri et al.. 2008bD . 

Finally, other very popular models, are the ones in- 
troduced in the spin glass literature. Probably the most 
studied example of such spin glasses is the p-spin mod el, 
defined by the Hamiltonian ( Gross and Mezard. 1984 ) 



H 



^ ] Jil i p Si 1 ...Si p , 



(40) 



where the SiS are Ising or spherical spins, p > 2 and 
Ji!,...,i p quenched random couplings with zero mean and 
variance p\/(2N p - 1 ). 

These models are certainly not realistic in terms of 
modeling microscopic degrees of freedom in a fluid, but 
they are representative of the class of systems with a ran- 
dom first order transition and have the advantage that 
a variety of computations can be performed without any 
approximation, and both their dynamic and static prop- 
erties can be investigated analytically in full detail, again 
yielding results as described above. Their dynamics can 
be studied in full detail, including various nonequilibrium 
conditions as described in SccfV] Another fruitful result 
concerns the interpretation of the nature of the low tem- 
perature phase in terms of replica symmetry breaking, 
so that conne ction with the field o f disordered systems 
can be made (jMezard et al. 19881 : iParisi. 20031) . Tech- 
nically, the thermodynamics of the p-spin can be solved, 
for p > 2, using a one-step replica symm etry breaking 
ansatz, see ( Castellani and Cavagna. 20051 ) for a review. 
This means that the probability distribution function of 
the overlap between states, the Parisi function P(q), has 
two peaks below Tk, one at q(T) > which quantifies 
the self-overlap within the states, and another one at 
g = implying that differ ent states are totally uncorre- 
cted (jMezard et al.. 1988f) . 

Let us now discuss the dynamical behaviour which re- 
sults from the above analysis of the free energy land- 
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FIG. 13 Correlation C(t) as a function of time for the p- 
spin model with p = 3 for several temperatures approaching 
Tmct, obtained from numerical solution of Eq. (|42|l . The 
curves show the appearance of a two-step decay characterized 
by several scaling laws discussed in the text. 



scape. Below Tmct, the system is in a metastable state 
from which it cannot esca pe, because free ener gy barriers 
diverge with system size ( Barrat et ai. 19961 ). This di- 
vergence is a direct consequence of the mean-field nature 
of the present set of approximations. Therefore, the re- 
laxation time diverges, within mean- field, at Tmct- The 
stability of these states can be analyzed by computing the 
free energy Hessian in the minima. One finds that states 
become more 'fragile' when T — > T MCTl are marginally 
stable at T = Tmct, unstable for T > Tmct- As a 
consequence, one expects that the dynamics slows down 
approaching Tmct fr om above as the landscape b ecomes 
more and more 'flat' ( Kurchan and Laloux. 1996T ). 

Indeed the dynam ics of many of th ese models can 
be analyzed exactly ( Cugliandolo. 20031) . In particular, 
mean-field p-spin models have been analyzed in great de- 
tail and provide a paradigm for mean field glassy dynam- 
ics. The equations of motion considered in the literature 
arc Langcvin equations, 



d Si {t) 
dt 



-n(t)si(t) 



dH 

dsi(t) 



ili(t), 



(41) 



where rji(t) is a Gaussian thermal noise of zero mean and 
variance IT given by the fluctuation-dissipation theorem. 

We focus on the spherical version of the model, 
and on the time autocorrelation function G(t) = 
7r Si( s *W s i(0))- Note that n(t) is the Lagrange mul- 
tiplier enforcing the constraint C(0) = 1. The equation 
of motion for C(t) at thermal equilibrium reads: 



dC{t) 
dt 



-TC(t) 



V_ 

2T 



dt'C p - x (t-t') 



dC(f) 
dt' 



(42) 



We will meet this equation again in the next section 
about mode-coupling theory. We shall then postpone a 
detailed study and just anticipate some results that will 
be derived later. At high temperature, the correlation 
function decays quickly to zero. Decreasing the temper- 
ature, the relaxation timescale increases and a two-step 
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FIG. 14 Evolution of V(q)/T = -^ F iq v - log(l - q) - q as a 
function of q for different temperatures across the dynamical 
transition. The solution with q > appears discontinuously 
at Tmct- 



relaxation emerges, see Fig H~3l where we have plotted the 
numerical solution of the previous equation. At Tmct 
the timescale r a corresponding to the slow relaxation di- 
verges algebraically, 



(T — Tj 



MCT) 



(43) 



where 7 is a critical exponent. The value of the plateau 
q = lirtij^oo C(t), called Edwards- Anderson parameter in 
the spin glass literature, satisfies a simple equation that 
can be obtained taking the infinite time limit of Eq. (|42|) : 



Tmct (iKirkpatrick et q/T 19891: IStevenson et al., 20061: 
ICavagna. 20091 ). 

Let us recap the overall picture arising from a mean- 
field analysis of the properties of the free energy land- 
scape of glasses. At high temperature the dynamics is 
fast and the system is in the liquid state. Approaching 
Tmct the dynamics slows down because of the appear- 
ance of incipient stable states. Decreasing the tempera- 
ture to Tmct, it takes a longer time to find an unstable 
direction, and thus to relax. Below Tmct there is an 
exponential (in the system size) number of states. The 
partition function and the thermodynamics are obtained 
by summing over all of them their corresponding Boltz- 
mann weight. This procedure is justified by the fact that, 
in a real finite dimensional system, the barriers between 
states should actually become finite. In the regime below 
Tmct, there is a competition between single state free 
energies that would favor the lowest free energy states, 
and configurational entropy that would favor the highest 
ones that are more numerous. Lowering the tempera- 
ture disfavors the entropy term and at Tx the system 
undergoes a phase transition where the sum in Eq. (|37|) 
is again dominated by only few terms corresponding to 
states with free energy density fx given by s c (fx,T) = 0. 
This transition corresponds to a bona fide 'entropy crisis' 
mechanism. 



1 



2T 2y 



(44) 



A graphical solution of this 
scntcd in Fig. [Ml where 



equation is pre- 
we plot V(q)/T = 



Jo dq' j3^7 — ^!fzq' p 1 ■ The minima of V are the 
solution of Eq. (14"4"1) . Clearly the minimum at q = 
is always present. Another solution, qEA, appears 
at T M ct = \Jv{v - 2)P- 2 {p - iy-P/2 and it can be 
interpreted as the long time limit of the correlation 
function inside one typical state. Since the states have 
an infinite lifetime (in mean-field theory) the system 
remains trapped forever into the one it started from. 
It is important to remark that qEA is discontinuous at 
transition, which leads to the two-step behavior shown 
in Fig. 1131 By contrast, qEA is continuously growing 
from at the transition in spin glasses. 

Note from Fig. [14] that at T MCT the Edwards- 
Anderson parameter is zero and, concomitantly, V(q) 
has a vanishin g second derivative at qEA- It is pos- 
sible to show ( Franz and Parish 19981 ) . that this is in- 
deed related to the fact that the free energy Hes- 
sian of the states below Tmct develops zero modes 
at Tmct- This behavior resembles very much the 
one of a spinodal transition. In fact, this anal- 
ogy has been fruitfully explored, for instance to de- 
scribe real space features of dynamic heterogeneity near 



2. Mode-coupling theory (MCT) 

The dynamical transition appearing upon approach- 
ing Tmct in RFOT landscapes is mathematically anal- 
ogous to the one predicted to occur in supercooled 
liquids by the Mode-Coupling Theory (MCT) of the 
glass transition, although the latter has a priori no di- 
rect interpretation in terms of a free energy landscape. 
This theory was introduced separately by Lcuthesscr 
(jLeuthesser. 19841) and Bengtzelius, G otze, Sjolandcr 
and coworkers ( Bengtzelius et al., 1984[ ). It has been 
used to describe and predict the average dynamics, in 
particular the dynamical structure factor and the self- 
diffusion, for moderately supercooled liquids and colloids. 
Recently, it has been generalized to describe dynamical 
correlations and some aspects of dynamic heterogeneity, 
as described in Sec. IIV.B.3I In Sec. IIV.B.4[ we discuss 
successes and limitations of MCT. 

Originally, MCT was developed using the pro- 
jection operator formalism ( Leuthesser. 19841 : 
iBengtzelius et al., 19841 ). A good introduction to 
this method c an be found in the book by Zwanzig 
( Zwanzig. 200lh . The starting point of the method is the 
derivation of the following equation for the dynamical 
structure factor F[k,t) of a single component atomic 
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liquid: 



d 2 F{k,t) k 2 k B T 
dfi 



F(k,t) + 



mS(k) 

/ d,TM(k,T)—F(k,t-T)=0. 
Jo at 



(45) 



Generalizations to mixtures and non-atomic liquid are 
also available. This is an exact equation whose in- 
puts arc the static structure factor S(k) = F(k,0) 
and the memory kernel M(k,r) for a given particle 
mass m, and temperature, T. In a second, crucial 
step MCT suggests a self-consistent approximation for 
the memory kernel M(k,r). It is possible to show 
that the memory kernel corresponds to the variance 
of the rand om force acting on the density fi eld, sec 
the review ( Reichman and Charbonneau. 20051 ). Thus, 
M(k, t) captures the effect of all degrees of freedom other 
than the density field on the density field itself. The phys- 
ical idea motivating MCT is to focus on the slow part of 
the random force. Technically, the path is in principle 
straightforward: one should identify the dominant slow 
modes, project the random force onto them, and derive 
the dynamical equations for their correlation functions. 
Of course, in practice this remains difficult because the 
number of slow modes is infinite. 

Within MCT, only the bilinear density products con- 
tribute to the slow part of M . After projection, the mem- 
ory kernel is expressed in terms of a four-point function. 
In a final approximation, this function is factorized as 
the product of two-point density functions F(k,t). This 
leads to the MCT self-consistent equations: 

d 2 F(k,t) , lS dF{k,t) k 2 k B T „, , 
= + v{k) — ^ + * , F(k, t) + 



dt 2 



dt mS(k) 
dF(k,r) 



dTM M CT(k,t-r) 

dr 

pk B T 



(46) 



M MC T(k,t) 



dk'|V k -k',k'r x 



167r 3 m 
xF(k',t)F(\k'-k\,t), 

Mc-k',k< = {(k-k')c(fc') + k- (k-k')c(|k-k'|)j47) 

where we have rewritten the result using the direct cor- 
relation function c(k) = (1 — 1/S(k)) /p. The effective 
friction term represents the effect of the fast degrees of 
freedom. 

This final expression clearly shows that MCT is 
a particular closure of the equations on dynamical 
correlation functions. It is similar in spirit to sev- 
eral other closure schemes used in physics, such as 
Kraichnan's 'Direct Interaction Approximation' for 
turb ulence, or various larg e-iV field theoretical meth- 
ods (jBouchaud et al., 19961 ). Indeed, field theoretical 
derivations of MCT have long been sought, but this 
is in fact still a very active area. The first pioneering 
works were published short l y after th e original MCT 



started from stochastic equations for the slow degrees of 
freedom of a liquid, the so-called non-linear fluctuating 
hydrodynamics, and rederived the MCT equations as 
a self-consistent, one- loop approximation. Motivations 
for the field-theoretical approach are that it provides a 
complementary way to derive MCT which is in principle 
more suitable to nonequilibrium generalizations, and 
perhaps to systematic improvement. Unfortunately this 
approach is plagued by difficulties related to the preser- 
vation of time- reversal symmetry in self-consistent 
loop expansions ( Mivazaki and Reichman. 20051: 
lAndreanov et at, 20061) . Very recent work aimed 
at getting fully consistent field-theoretical deriva- 
tions of MCT equations ( Kim and Kawasaki. 20081 
iNishino and Havakawa. 20081) . but this is technically 
more involved than could naively be anticipated. 

Within MCT, dynamical correlation functions 
are obtained by numerical integration, once the 
static structure factor of the liquid is known. This 
kind of analysis was performed on a large num- 
ber of different glassy liq uids such as Lennard - 
Jones (jBengtzelius. 19861: iNauroth and Kob. 1997 ) 
hard spheres sy stems (iBarrat et al, 198£ ; 



iBarrat and Latz. 19901: iFoffi et al., 2003[ ). or silica 
melts ( Sciortino and Kob. 200ll) . and it has revealed 
a common behavior. It is found that small changes 
in S(k) lead, at high density or low temperature, to a 
great variation in F(k,t) that resembles very much the 
one shown in Fig. [13] for the correlation function of the 
p-spin model. The MCT equations display a dynamical 
glass transition to a phase where the average density of 
the liquid remains frozen in an amorphous profile. 

The similarity with the p-spin model is not casual, 
and there is indeed a deep connection between the 
two. A first, somewhat technical, way to unveil it con- 
sists in simplifying the wavevector dependence of the 
equations assuming that the integral over k is domi- 
nated by values close to kg where the structure fac- 
tor has a str ong peak. This s o-called 'schematic' ap- 
proximation (ILeuthesser. 19841 : IBengtzelius et a~, 19841: 
iGotze. 19991) . yields a simplified equation of motion for 
F(k ,t) that reads: 



v{k 



dF{kp,t) 
' dt 



k 2 k B T 
mS(ko) 

kpA 2 

8tt 2 P 



F(k ,t) 



dt'F 2 {k ,t-t') 



(48) 

dF(ko,t') 
dt' 



works were published short l y alter th e original M(J1 
derivation ( Das et al, 19851 : iDas. 1990l) . The authors 



where A is the area under the peak of S(k) at 
fco- A simple change of variables maps this equa- 
tion to that of the 3-spin model, Eq. (|4"2"|) . This re- 
lation with fully-connected models suggests that MCT 
should be interpreted as a mean-field approximation. 
Note that this does not imply that MCT becomes 
exact in the limit of large spatial dimensionality, as 
shown by recent calculation s ( Ikeda and Mivazaki. 2010t 
ISchmid and Schilling. 2oTol) . 

The solution of MCT equations displays a very rich 
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phenomenology as seen in Fig. [13] There are three time- 
regimes. A fast relaxation toward a plateau, whose value 
depends on k, a slow relaxation close to the plateau, 
called /3-regime, and finally the s tructural relaxation on 
the timescale of the a-regime; see (jKob. 2003tfDas. 20041) 
for more details. In the following we shall denote e the 
relative distance from the transition. For molecular liq- 
uids the control parameter is the temperature and, hence, 
e = (T — T c ) /T c ; for colloids the control parameter is the 
density and so e = (<f> c — 4>)/4> c . 

Keeping track of the wavevector dependence, the de- 
tailed properties of the dynamics in the three regimes are 
as follows: 

(1) Fast relaxation - Some degrees of freedom relax on 
timcscalcs of the order of tq, even close to the transi- 
tion. This regime is identified taking the limit e — > and 
keeping t finite. In this case F(k, t) approaches a plateau 
at long times whose value is denoted fkS(k). The non- 
ergodic parameter, is the fraction of density fluctua- 
tions that become frozen at the transition. At large times 
the behavior of F(k,t) is: 



F{h,t)*f k + 



h(k) 
t a 



t > TQ. 



The exponent a satisfies the equation 

r 2 (i- c 



T(l -2a) 



A. 



(49) 



(50) 



where A is a number that can be computed using the 
structure fact or only. Its expression is complicated, see 
( Gotze. 19851 ). The previous equation implies < a < 
1/2. 

(2) ^-regime - In this sector, the timescale diverges 
as Tfj ~ e" 1 / 2 " and the dynamical structure factor scales 
as: 



F(M)« fk + V~zh(k)g(t/Tf}), 



(51) 



where g(x) oc x~ a for x -C 1 and g(x) oc x b for x 3> 1. 
Note that the previous expression implies that all the 
k dependence factorizes and is contained in h(k) only, 
the so-called 'factorization property'. The exponent b 
satisfies the equation 



r 2 (i + ) 
r(i + 26) 



A. 



(52) 



which implies that < b < 1. 

(3) a-regime - In this sector the timescale diverges as 
r a ~ e -7 , where 7 = l/2a + 1/26. The factorization 
property does not hold anymore except for small t/r a 
because the solution has to match the one found in the 
/3-rcgime. 

We refer again the reader to Fig. [T3] for a visual illus- 
tration of the different time regimes predicted by MCT 
for dynamic structure factors in supercooled liquids. 

Mode-coupling theory provides predictions also for 
other correlators such as the self-intermediate scattering 
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FIG. 15 Fit of the evolution of the equilibrium structural re- 
laxation time of a Lennard-Jones liquid (temperature is var- 
ied) and a hard sphere fluid (density is varied) with the pre- 
dicted MCT algebraic divergence, r a ~ e~ 7 , where e is the 
reduced distance to the transition. The plateau at low e shows 
that T a remains finite at e = 0, and the data directly indi- 
cate that the transition can actually be crossed at thermal 
equilibrium and is thus avoided. 



function from which the mean-squared displacements and 
thus the self-diffusion coefficient can be obtained. The 
previous properties remain essentially unaltered and all 
correlators display quite similar scalings. 

All these predictions have been tested in great de- 
tail in numerical simulations and in experiments both 
on molecular liquids and in colloids. It has also been 
shown that adding corrections to MCT does not spoil 
the main predicti ons and the universalit y of MCT has 
been e stablished (lAndreanov et ai, 2 009). Different re- 
views (jGotze. 19991: iKob 20031: iDas. 2004| ) have already 
appeared on these tests. When fitting data using MCT, 
a central difficulty arises from the fact that the actual 
transition is not present, as expected from its mean-field 
nature. The absence of a genuine mode-coupling singu- 
larity is undisputed for molecular liq uids both in simula- 
tions and experiments ( Gotze. 2008[ ). Recent numerical 
and experimental works sugg est that the same situation 
holds in hard sphere systems dSanten and Krauth. 20001 : 
iBrambilla et al, 20091: iFlenner and Szamel. 20101) . This 
is illustrated in Fig. [T5] where the predicted MCT alge- 
braic divergence of the structural relaxation time for a 
Lennard-Jones liquid and a hard sphere fluid are super- 
imposed on numerical data. While the fit is accurate 
over a window of 2 to 3 decades, it clearly fails to cap- 
ture the low temperature or large density regimes of these 
systems. Clearly, therefore, mean-field concepts cannot 
directly be applied to understand the glass transition, 
and a more refined analysis is needed. 

In conclusion, MCT predicts a transition where the 
system has a dynamical arrest, particles stop to diffuse 
and the density becomes frozen around an amorphous 
profile. Additionally, MCT yields non-trivial predictions 
for the behavior of dynamical correlators that serve as a 
guideline in the study of moderately supercooled liquids. 
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3. Dynamical correlations within Mode Coupling Theory 

We introduced MCT as a dynamical theory for two- 
time correlation functions. However, the recent surge of 
interest on high-order correlation functions as probe of 
dynamic correlations and dynamic heterogeneity suggests 
that it could be interesting to develop an MCT approach 
also for multi-point correlation functions. 

The MCT transition was originally described as a 
local phenomenon, the sel f-consistent blocking of the 
particles in their 'cages ' dGotze. 19991) . On general 
grounds, a diverging relaxation time is expected to arise 
from processes involving an infinite number of parti- 
cles (leaving aside the case of quenched obstacles), as 
discus sed in Sec. IIII.BI and established by rigorous re- 
sults (jMontanari and Semeriian. 2006bl) . Actually, even 
the cage mechanism requires some kind of collective be- 
haviour: in order to be blocked by one's neighbors, the 
neighbors themselves must be blocked by their neighbors 
and so on until a certain scale that, intuitively, should set 
the relaxation timescale of the system. So one expects 
that even within MCT, 'cages' should in fact be de scribed 
as spatially correlated objects ( Biroli et al.. 20061 ). It has 
indeed been shown that all susceptibilities and correla- 
tions defined in Sec. IIII.BI such as X4 and xt, diverge at 
the MCT dynamic singularity. 

Historically, the 'local-cage' point of view 
was challenged in the c ontext of mean-field 
disordered system s in dFranz and Parish 2000 ) 
see (IKirkpatrick and Thirumalai. 198E ; 



IKirkpatrick et al.. 19891) for early results), since these 
models are analogous to schematic MCT equations. 
Franz and Parisi argued that a dynamical suscep- 
tibility similar to Xi{t) diverges at the dynamical 
mode-coupling trans ition. The first full MCT analysis 
was developed in (jBiroli and Bouchaud. 20041 ) . using 
a field-theoretical approach. This clearly showed the 
existence of a diverging lengthscale within MCT. 
Later, a different suscep tibility, Xqo (li 1 was intro- 
duced ( Biroli et al. 20061 ) . which quantifies the response 



of the dynamical structure factor to a static oscilla- 
tory external perturbation with wavevector qo. For a 
perturbation localized at the origin, £/(x) — UqS(x), 
one finds SF( qi ,y,t) = U J d^e^x^!^)- This 
susceptibility is akin (although not exactly related) to a 
three-point density correlation function in the absence 
of the perturbation. Although quite different from 
the four-point functions considered previously in the 
literature, Xqo (li ; i s expected to reveal the existence 
of a dynamical correlation length of the homogeneous 
liquid as well. Physically indeed, Xqo( ( lii^) measures 
the influence of a density fluctuation at a given point 
in space on the dynamics elsewhere. Additionally its 
scaling form is not affected by complications due to con- 
servation laws, as is the case for four-point correlators, 
see Sec. IIII.CI Recent work on Kac glassy models also 
found a similar diverging len gth at the MCT transition 
( Franz and Montanari. 20071 ) using replica techniques. 



Let us su mmarize the critic al properties found for 
X qo (qi,£) in ( Biroli et at. 20061 ). As for the dynamical 
structure factor there is a different critical behavior in the 
/?- and a-regimes, although there is a unique diverging 
correlation lengthscale £ ex e -1 / 4 . 

(1) In the /3-regime, one has: 



Xqo(qi>i) 



r«5 



S{qi)h(qi)gp 



L L 



(53) 



where gp is a scaling function and T a posi- 
tive constant, whi ch are both obtained quantitatively 
( Biroli et al.. 20061 ). In particular, gp behaves as a ~ t a 
for small t/rp, and as ~ t b for large t/rp. 

(2) In the a-regime the critical behavior of Xqo( ( li)^) 



Xqo(qi 5 i) 



s(r 9o 2 /^) 
V~<V^ + rq 2 ) 



9a,q t [ 



(54) 



with 3 a certain regular function with S(0) ^ and 
S(d 1) ~ \ jv such that Xq behaves as q§ for large qo, 
independently of e. Also, g am {u C 1) = S(qi)h(q{)u b , 
to match the /3-regime, and g a {u 2> l,<?i) — > 0. 

The behavior of four-point quantities, such as G4 
and X4j is more complicated because of ensemble de- 
pendencies and the in fluence of conservation laws, see 
( Berthier et al.. 2007al ) for a detailed discussion. It is 
found that G4 and X4 should have a similar critical be- 
havior, but not too close to the MCT transition, cross- 
ing over to a distinct behaviour in its vicinity, due to 
conservation laws. Therefore, for decisive tests of MCT 
predictions regarding multi-point functions, the quantity 
X qo (qi,£) should be preferred, but no such study has 
been reported yet. 

The study of critical properties of dynamical 
correlations for glassy liquids and, hence, the com- 
parison with MCT pr edicti ons is still in its infancy 
(iToninelli et al.. 2005bl: ISzamel and Flenner. 2006 ; 
Berthier et al. 2007at iBerthier et al.. 20"07b : 

Stein and Andersen. 20081 : M enner and Szamel. 20101 ). 



At the time of this writing many questions remain 
open and deserve further studies. The determination 
of the dynamic correlation length is subtle and because 
of the complication brought about by the existence of 
conserved quantities G4 could be difficult to analyze. 
Certainly a numerical study of Xqo (li j *) would be 
very im portant. Furth er more the role of finit e size 
effects jBerthier. 2003bt iKarmakar et al. 20091 ) , the 
dependence of the results on the observable used to 
probe dynamical correlations, the universality of the 
results, the pos sible anisotropic characte r of dynamic 
heterogeneity ( Flenner and Szamel. 20071 ) all figure 
among open questions needing further investigation. 



4. Current status of MCT 

The first and most important drawback of mean-field 
dynamics and MCT is that the MCT transition it de- 
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scribes is not observed in real materials. Additionally, 
a comparison between Eq. ^ and (pi!?]) makes it clear 
that MCT cannot be used to describe viscosity data 
close to the experimental glass transition T g , since it 
does not even predict thermally activated behaviour. 
Worse, MCT predicts a transition at which the system 
freezes completely: not only a fraction of the density fluc- 
tuations get frozen but also self-diffusion gets arrested 
( Fuchs et al.. 19981 ) . This i s a theoretic al artifact as it 
can be rigorously proven ( Osada. 19981 ) that the self- 



diffusion coefficient cannot go continuously to zero at 
thermal equilibrium and finite temperature and pressure 
(this excludes the case of hard particles at close packing). 

With all these major failures, then, why on earth 
should one continue to study and talk about MCT? We 
provide several reasons. 

It is now recognized that the MCT transition must 
be interpreted as an approximate theory of a crossover 
taking place in the dynamics. Actually, there are 
many other phy sical examples, such as the Kondo mode l 
( Hewson. 1997 ) or spinodal points ( Debenedetti. 19961) . 
where crossovers become transition in approximate self- 
consistent theory. The fact that the transition is sharp 
in the theoretical treatment allowed the derivation of a 
variety of scaling laws which are as many predictions that 
can be tested in real materials, or, more prosaically, sim- 
ple formula that can be used to 'fit the data'. Part of the 
success of MCT has been its ability to propose such clear- 
cut predictions (and fits), and graphical representations 
along which data could be analyzed. 

Indeed, one finds in the literature scores of papers 
where an MCT analysis of data is performed. Since the 
transition is (at best) only 'avoided', its 'crossover' na- 
ture offers quite a lot of flexibility for fitting, and judging 
the quality of the fits is often a difficult (and subjective) 
exercise, while 'negative' results can always be attributed 
to 'preasymptotic' corrections rather than deficiencies of 
the theory itself. This has led to many controversies in 
the literature which seem to persist to this day. 

A major achievement of MCT is the possibility to ap- 
ply the same formalism to different materials and theoret- 
ical models, basically starting from the microscopic inter- 
actions between atoms or molecules. This is again a rea- 
son for the success of MCT: each time a new model with 
a different kind of interactions or chemical composition 
is defined, MCT can be used to analyze its behaviour, 
and possibly predict new qualitative trends for the time 
dependence of correlation functions and their evolution, 
even in cases where several control parameters are rele- 
vant (e.g. mixtures, attractive vs. repulsive interactions, 
complex fluids, etc.). Striking and successful results and 
predictions, later confirmed by simulations and experi- 
ments, have been obtained in several cases: for example 
for the dynamical phase diagram of attractive colloids, 
for the behavior of the non-ergodic parameter as a func- 
tion of the wave- vector for glassy liquids, for the treat- 
ment of molecules or particles with non-spherical shapes , 
see the reviews (iKob 20031: iDas. 2004 iGotze. 20081 ). 



Thus, MCT is always a useful starting point when a new 
system with unknown behaviour arises. 

Although the MCT predictions are limited to a modest 
time window of 2-3 decades corresponding to the onset 
of glassy dynamics in molecular liquids, this time win- 
dow is actually the most experimentally rele vant in col- 
loidal materials ( Cipelletti and Ramos. 20051 ) , since typ- 
ical microscopic timescales for colloidal particles are in 
the millisecond range (instead of the nanosecond range 
for atoms). Thus, MCT performs much better in soft 
matter systems, to the point that actually observing de- 
viations from MCT predictions in an experiment can 
represent a challenge. Even for the canonical and well 
studied system of colloidal hard spheres, it was only re- 
cently suggested experimentally that the MCT transition 
i s avoided in the same manner as in molecular liquids 
(jBrambilla et al. 20091 ). Thus we believe MCT will con- 
tinue to be a useful tool, given the current rapid devel- 
opment in the synthesis of new colloidal particles. 

On the theoretical side, it is now clear that MCT ha s 
a status of a mean- field theory ( Andreanov et al. 20091) . 
As such, one expects major changes once fluctuations will 
be taken into account. As we pointed out already some 
kinds of fluctuations wipe out the sharp MCT transi- 
tion and makes it a cross-over. Moreover, exactly as for 
critical phenomena, one expects critical fluctuations to 
change the value of the MCT exponents below an upper 
cri tical dimension that was determined to be eq ual to 
8 (jBiroli and Bouchaud. 20071 iFranz et al. 20101 ). The 
role of fluctuations on the MCT transition is a current 
topic of research. 

As a conclusion, MCT has clear and well-understood 
limitations and it will be never possible to test its pre- 
dictions in a very sharp way because it is not related to a 
true phase transition but, likely, just to a crossover. Still, 
its overall efficiency and flexibility, its ability to deliver- 
ing actual predictions makes it useful. For this reason 
it continues to be developed, applied and generalized to 
study many different physical systems and situations, in- 
cluding aging systems and non-linear rheology of glassy 
materials, see Sec. IVl 



5. Quantitative computations using replica 



In Sec. IIV.B.11 we presented the theoretical picture 
emerging from solving mean-field models (or geometries) 
of the glass transition, and we found the resulting sce- 
nario rich and encouraging: it generically supports the 
existence of a configurational entropy vanishing transi- 
tion associated, at a higher temperature, with a dynam- 
ical transition, a la MCT, which corresponds to the ap- 
pearance of incipient metastable states. 

Still, to make a connection with experimental results 
at least two main issues need to be addressed. First, one 
has to transform this set of mean-field ideas into a work- 
ing tool able to produce quantitative calculations for the 
case of supercooled liquids. Second, dynamics within a 
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rugged landscape must be revisited in order to explain 
the crossover nature of the MCT transition, and the ex- 
istence of a regime where dynamics is thermally activated 
and the viscosity increases in a super- Arrhenius manner 
which is incompatible with the algebraic divergence pre- 
dicted by mode-coupling theory. In the present section 
and the next, we shall review briefly these two lines of 
research. 

First, we shall focus on the approaches that have 
been devised to obtain quantitative microscopic pre- 
dictions. In fact, MCT can be seen already as such 
a tool, but it is limited to the regime above Tmct- 
Below this temperature it cannot be applied any- 
more. At the time of this writing, a quantitative 
microscopic theory of the dynamics valid below Tmct 
is still lacking, see ( Bhattacharva et al., 20051) for a 



possible attempt. An alternative strategy is to leave 
dynamics for a while, and to turn to thermodynamics. 
The idea is to compute approximately the configura- 
tional entropy, the Kauzmann temperature, and the 
plateau value of dynamical correlation functions in 



the glass phase dMonasson. 1995; Franz et al, 1998a 


Mezard and Parisi. 1999 


; Parisi and Zamponi. 2005 


Mezard and Parisi. 201C 


). 



All the approaches developed to compute static prop- 
erties of systems characterized with a rug ged landscap e 
make use of replica in one way or another ( Parisi. 2003|) . 



Physically, the reason is that one aims at describing (or 
at least counting) mctastable states which all have amor- 
phous density profiles. This is similar to the spin glass 
case where the amorphous order is not revealed by look- 
ing at magnetization profiles, as discussed in Sec. IIII.B.ll 
Inspired once more by the physics of disordered systems, 
the idea is again to let the system itself indicate what 
are these states, and project distinct copies of the sys- 
tem to 'recognize' the metastable states. In the absence 
of quenched disorder, however, it is not enough to repli- 
cate the system, one is also forced to physically 'couple' 
the different copies of the system using an appropriate 
field. 

It is useful to fi rst implement this idea for two copies 
( Monasson. 19951) . 1 and 2, of the system, characterized 
by density profiles p\ (x) and pi (x) . Using notations from 
Sec. HV.B.ll we write the free energy density of a single 
copy of the system as: 
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(55) 



where H [pi] is the microscopic Hamiltonian. Let us now 
use the second copy of the system to scan the locally sta- 
ble configurations of the first one. To do so, we introduce 
a quadratic coupling of strength g > between the two 
copies and compute the new free energy 
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The free energy /2[pi] will be small when p\ is a config- 
uration which corresponds to a mctastable state. There- 



fore, sampling all configurations of p\ weighted with 
cxp(— f3f2[pi]) is a procedure to scan all metastable 
states, so that 
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(57) 



is the average free energy of all metastable states; here 
Z = J Dpi exp(— /3f2[pi]) is a normalization. When it ex- 
ists, the difference between the total free energy density 
of the system, /, and the free energy density of the min- 
ima, /mota, is related to the number of metastable state 
available to the system, see also Eq. (|38|) . As a conse- 
quence, this free energy cost is in fact equal to —Ts c . 
The lesson to be learnt from this example is that the in- 
troduction of identical copies of the system allows one 
to com pute directly properties of the free energ y land- 
scape dMonasson. 19951 : iFranz and Parisi. 1998 ). It is 



important to remember, however, that the limit g — ¥ 
in Eq. (|57| has a strong mean-field flavour, since 
genuine meta stable states only exist in this limit, see 
( Parisi. 20031 ) for a more precise discussion. 

Several quantitative approaches have been 
developed and are based, in one way or an- 
other, on pro cedures similar to the one outlined 
see dMonasson. 1995 : 



above, 



Franz et al, 1998at 



iMezard and Parisi. 19991: Parisi and Zamponi. 2005 ) 



In the scheme of ( Monasson. 1995[) . one introduces m 
copies of the system constrained to be in the same 
free energy minimum. Technically, this corresponds to 
take m copies of the density configuration p\. This 
generalizes the partition function in Eq. (|38[) to 



Z m = lim — 

g^o Z 
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df exp 
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associated to the replicated free energy ip(m,t) = 
— Tlog Z m . Note that m only enters the first term, as all 
systems are identical and are characterized by the same 
metastable states. Repeating the saddle point calcula- 
tion yields 



/ = 



dip(m,T) 
dm 



m 2 d(ip(m,T)/m) 



T 



dm 



(58) 



which shows that the configurational entropy ,s c can be 
accessed by computing the thermodynamic properties of 
a system of m coupled replicas, provided, as is usual 
within replicas calculations, that the number of copies of 
the system is analytically continued to non-integer val- 
ues, as implicitly assumed in Eq. (l58l) . 



A remarkable achievement ( Mezard and Parisi. 1999t ) 
is that not only the properties of the liquid state above 
T K can be computed analytically from Eq. (|S"8"]) . but 
also the ones of the 'ideal' glass below Tr-. Recall 
from Fig. [T^] that the Kauzmann temperature is defined 
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by Tk = {ds c /df) 1 for the value of / correspond- 
ing to s c = 0. For the replicated system, one gets 
= m{ds c /df)^ 1 and the location of the transition 
depends on to. For values of to < 1, one typically finds 
that Tjf < Tk- This means that the properties of the 
replicated liquid with to ^ 1 are deeply connected to the 
ones of the non-replicated glass with m = 1. Indeed, as- 
suming that the transition has a nature similar to the one 
found in mean- field models, one can use the fact that the 
free energy is continuous at Tk and glass states below 
Tk have s c (T < Tk) — to obtain the free energy of the 
(non-replicated) glass as 



/gla SS (T) = ^TO*(T),T)/TO*(T), 



(59) 



where m*(T < Tk) < 1 is self-consistently defined by 
T (m*) = T and ^,( TOj r) is the energy of the replicated 
liquid defined above. 

To summarize, starting from the hypothesis that the 
free energy landscape of supercooled liquids resembles 
the picture gained from mean-field models and geome- 
tries, one introduces replicas as a useful mathematical 
tool to probe the thermodynamics of systems both above 
and below Tk- One ends up with an additional varia- 
tional parameter, m*(T), to describe the low tempera- 
ture phase, which is formally strictly equivalent to the 
additional parameter to entering the one-step replica 
symmetry schem e needed to solve mean-field models 
( Monasson. 19951 ). Note that in this approach, the nature 
of the broken symmetry is a starting hypothesis rather 
than a natural outcome of the calculations. 

In practice, of course, some approximations must be 
made to compute the free energy of the replicated liquid 
at low temperature, which make heavy use of liquid 
state theory and which might well depend on the studied 
system. The main outcomes are the calculation of the lo- 
cation of the Kauzmann transition, the thermodynamic 
properties of the liquid (in particular the configurational 
entropy), and the glass (ground state energy, specific 
heat, structure). These microscopic computations 
have been developed for a variety of glass-forming 
liquids, such as m onoatomic or mixtu res of Lennard- 
Jones particles dColuzzi et al. 20001). soft -spheres 
dMezard and Parisi. 19991: IColuzzi et al, 19991 ). hard 



spheres ( Franz et al, 1998at Parisi and Zamponi. 2005f) . 
sticky hard spher es ( Velenich et al. 20061 ). silica 
( Coluzzi et al, 20021 ). Perhaps the most impres- 
sive achievement is the detailed description of the 
large volume fraction b e haviour of hard sphere s 



( Parisi and Zamponi. 20051 : iParisi and Zamponi. 2010() : 
the glass transition was located, the equation of state 
and structure of the glass obtained up to the 'Glass Close 
Packing' density, which can be seen as a firm definitio n 
of the notion of random close packing ( Bernal. 19591) . 
whose location can t hen be predicted accurate ly in any 
spatial dimension (IParisi and Zamponi. 20061 ) or for 
binary mixtures ( Biazzo et al., 2009( ). with excellent 
agreement with experimental results and simulations. 



Although this quantitative side of RFOT theory 
is a most desirable feature, assessing quantitatively 
the quality of the results is not easy as experiments 
and simulations typically fail to approach the Kauz- 
mann transition. Even then, cases are known where 
a transiti on is predicted i n a regime where n one is 
expected ( Thalmann. 20021 : IColuzzi et al., 20"02l) . sug- 
gesting that the quality of the approximations used 
to obtain quantitative res ults plays an important role 
( Ikeda and Mivazaki. 2010l) . Thus, when accurate re- 
sults are sought, the problem might well becom e tech- 
nically quite involved ( Parisi and Zamponi. 20~10h . 



6. Scaling arguments beyond mean-field theory and 
point- to-set lengthscale 

The quantitative calculations described in the previous 
section remain mean-field in nature because they com- 
pute the properties of the supercooled liquid state, as if 
it were formed by a collection of states with infinite life- 
time, as is the case in mean-field models. This is clearly 
incorrect: two thermodynamically stable states cannot 
coexist and have different intensive free energy at finite 
temperature, otherwise the system would nucleate from 
the one with highest free energy to the one with low- 
est, showing that the highest is in fact not a stable state. 
Furthermore, an exponential number (in the system size) 
of stable states seems impossible: for large sizes there 
would not be enough boundary co nditions to select one 
state from the other ( Fisher. 20031 ). 

Additionally these calculations cannot address 
the connection to dynamical properties, a cru- 
cial missing ingredient for a theory of the glass 
transition. Presently, t here only exist phe- 

nomenological argum ent s dKirkpatrick et al. . 19891: 
IXia and Wolvnes. 20001: IBouchaud and Biroli. 200l. 
backed by m icros copic computations jDzero et al, 20051: 
iFranz. 20061 iFranz and Montanari. 20071) . that 
yield a possible scenario for the dynamics , 
dubbed 'mosaic state' in ( Kirkpatrick et al, 19891) . 
Since this aspect o f RFOT theory was re- 
viewed recently iLubchenko and Wolvnes. 2007t 
iBiroli and Bouchaud. 20091) . we shall be brief. Schemat- 
ically, the mosaic picture states that, in the regime 
Tk < T < Tmct, the liquid is composed of domains of 
linear size £. Physically, the lengthscale £ represents the 
lengthscale above which it does not make sense to talk 
about metastable states anymore. 

The way to measure, or even to define precisely, the 
mosaic lengthscale £ was not clear from the way it 
was initially introduced. Recently, the so-called 'point- 
to-set' correlation lengthscale was defined both in the 
context of RFOT t heory as a practical measu re of the 
mosaic lengthscale (IBouchaud and Biroli. 2004), and in 
more general settings ( Montanari and Semeriian. 2006bt 



iMezard and Montanari. 20061 ). The point-to-set length 
is a measure of the spatial extent over which the effect of 
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equilibrium amorphous boundary conditions propagate. 

To understand the origin of this lengthscale, we con- 
sider the following 'experiment' where, starting from an 
equilibrium configuration of the system, we freeze the 
position of all particles outside a cavity of radius R. We 
then let the system evolve in the presence of this con- 
straint, which acts as a pinning field. The point-to-set 
lengthscale is defined as the minimal cavity radius such 
that the pinning field has no influence at the center of 
the cavity. As such, it is a measure of the many-body 
correlations between a point (the center of the cavity) 
and a set of points (the pinned boundary). It is also 
important to emphasize the similarity between this cav- 
ity procedure, and the coupling between copies applied 
in Eq. ([56"|) . While the latter was homogeneous, the 
former is spatially inhomogeneous, and quantifies how 
far in space the coupling b etween states can propagate 
( Kurchan and Levine. 20091 ). 

The constraint on the boundary of the cavity in fact 
acts as a 'template' for the particles inside the cavity, 
whose effect can be evaluated as follows. By selecting a 
different state and deforming it only close to the bound- 
ary to satisfy the constraint, one would get a free energy 
cost of the order of TR e with 9 < 2. However, doing so, 
one would also gain entropy as the system could explore a 
multiplicity of different states, giving rise to a free-energy 
contribution —Ts c R 3 . Entropy obviously gains on large 
lcngthscales, R 3> 1, while interface cost dominates at 
small R. Therefore, the crossover in the competition be- 
tween these two terms takes place for a lengthscale R = £ 
obtained by equating the two terms, 
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(60) 



In a real liquid, where there is no cavity, one can conjec- 
ture that there is a self-generated dynamical boundary 
condition acting on each patch of lengthscale f . 

The dynamical counterpart to this argument is as fol- 
lows. Dynamically, the configurational entropy on scales 
smaller than £ is too small to stir the configurations effi- 
ciently and loses against the dynamically generated pin- 
ning field due to the environment. By contrast, ergodic- 
ity is restored at larger lengthscale. Hence, the relaxation 
time of the system is the relaxation time, t(£), of a finite 
size regions of the system. It is only after this long, but 
plausible scries of arguments, that barriers encountered 
during relaxation finally become finite and involve a fi- 
nite number of particles, unlike in the original mean-field 
treatment of the landscape where barriers diverge with 
system size. 

Now, assuming thermal activation over energy bar- 
riers which are supposed to grow with size as <f^, 
with ij) > 9, one predicts fi nally using Eq. ([60]) 
that (|Bouchaud and Biroli. 2001 
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where c is a constant. This argument is rather 
generic and therefore not very predictive. Re- 
cent micro scopic c omputations (|Dzero et ai. 2005 ; 
Franz. 20061: Montanari and Semeriian. 2006a : 



Franz and Montanari. 20071 ) attempted the compu- 
tation of the exponents 9 and tp, directly addressing 
analytically the problem of the cavity described above. 
Equations for the overlap profile between the initial 
'template' configuration and configurations thcrmalizcd 
in the presence of the pinning boundaries were obtained. 
These calculations confirmed the existence of a non- 
trivial crossover lengthscale above which the overlap 
inside the cavity vanishes, indicating that order does not 
propagate on lengthscales much larger than £. 

The results are unfortunately not yet conclusive. Al- 
though part of the computations can be jus tified and con - 
trollcd by using for instance Kac models ( Franz. 20061) . 
other parts involve uncontrolled replica symmetry break- 
ing schemes. The calculations provided the estimate 9 = 
2. Note that some other phe nomenological argument s 
suggest the value of 9 = 3/2 (jKirkpatrick et ai. 1989( ). 
There are no detailed com putation available for ijj, only 
the suggestion that ip = 9 (jKirkpatrick et ai, 19891) . 

Note that using the value 9 = 3/2 with 9 = 
tp simplifies Eq. (|61|l into a form that is well- 
known experimentally and relates logr Q directly to 
l/Sc, which is the celeb rated Adam-Gibbs rela- 
tion (lAdam and Gibbs. 19651) between relaxation time 
and configurational entropy. As discussed in Sec. III.A[ 
such a relation is in rather good quant itative agree- 
ment with many experim ental results ( Angell. 1997t 
iHodee. 19971 IJohari 2000l) . RFOT theory, therefore, re- 
formulates an d generalizes the mecha nism suggested by 
Adam-Gibbs ( Xia and Wolvnes. 2000h . Furthermore, us- 
ing the fact that the configurational entropy vanishes lin- 
early at Tk, Eq. (|6 1 1) becomes similar to the VFT di- 
vergence of Eq. ([2]), with the identification between two 
important temperatures, 



To = Tk, 



(62) 



(61) 



which embodies the deep connection between thermody- 
namics and dynamics characterizing RFOT theory. The 
above equality between two temperatures that are com- 
monly used in the description of experimental data con- 
stitutes a central achievement of RFOT theory, since it 
accounts for the empirical relation found between the ki- 
netics and the thermodynamics of supercooled liquids. 
It should be kept in mind, however, that experiments 
have not established its validity beyond any doubts, as 
discussed in detail in Sec. III. At 

Wolyncs and co-workers have obtained several other re- 
sults using phenomenological arguments based on RFOT. 
Two remarkable ones are: the relation between fragility 
and specific heat jump at the glass transition and stretch- 
ing exponent /? of time-dependent correlation functions, 
and the speeding up of th e dynamics close to a free sur- 
face (recently observed in ( Ashtekar et ai, 2010h ) . These 
predictions and several others are discussed in detail in 
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(|Lubchenko and Wolvnes. 20071 ). 



7. Current status of RFOT theory 

We described RFOT theory 'patchwork' (not 

to say 'mosaic') of apparently distinct theoretical ap- 
proaches to the glass problem: Adam-Gibbs theory, 
mode-coupling theory, mean-field spin glass theory and 
replica approaches, and the mosaic state scaling picture. 
As such, RFOT theory is clearly an impressive theoret- 
ical piece of work, which gives a very consistent overall 
scenario for the glass transition and nonequilibrium phe- 
nomena related to the glassy state, based on peculiar 
features of the free energy landscape as well as tools to 
perform microscopic calculations. 

Coming from high temperatures, dynamics primarily 
slows down because there appear incipient metastablc 
states, in a restricted temperature window described in 
full microscopic detail by mode-coupling theory. Decreas- 
ing further the temperature, the dynamics becomes dom- 
inated by the thermally activated barrier to be crossed 
from one metastablc state to another, in a way consis- 
tent with the deep relation between dynamical correla- 
tion length and timescale discussed in Sec. IIIII In this 
regime, the thermodynamic behaviour can also be de- 
scribed at the thermodynamic level using replica calcu- 
lations which predict the existence of a hnitc temperature 
thermodynamic transition towards a genuine 'ideal' glass 
state. A description of the dynamics near the transition 
exists, but contains several steps that heavily rely on em- 
pirical scaling concepts. 

There are of course several weaknesses in this con- 
struction. First, although we attempted here to give 
a unified view, the theory is still pretty much made 
of distinct pieces that do not necessarily smoothly ht 
together and have perhaps no strict boundaries. For 
instance, the details of the crossover between MCT 
and activated regimes are not well understood. In 
early works, attempts were made to include 'hopping 
effects' in mode-coupling theory, deriving expressions 
for the memory kernels which t ransform the sharp MCT 
transition into a cross over ( Das and Mazenko. 19861 : 
iGotze and Sjogren. 19871) . How ever, the status of these 
'extended' MCT is debated (jAndreanov et ai, 20061 
ICates and Ramaswamv. 20061 ). Moreover, there are 
also strong indications that thermal activation is 
in fact already at play in the t emperature regime 
usually de scr ibed by MCT (iDennv et ai, 20031: 
iHeuer. 20081 iBerthier and Garrahan. 2003bl ). Fi- 
nally, recent works attempted to include non- 
pertur bative processes i n the MCT descrip - 
tion (jMaver et ai. 20061: iBhattacharva et al., 200a) . 



complete and solid version of RFOT theory remains to be 
worked out. This is especially true for the 'mosaic' part of 
RFOT theory which yield dynamic predictions, namely a 
Vogel-Fulcher divergence of the relaxation time. This im- 
plies that opponents can criticize RFOT theory because 
it contains too many uncontrolled assumptions, while 
supporters can always hide a weak point as resulting from 
some approximation rather than from the approach itself. 
We do not see how this issue can be resolved, as data 
themselves do not allow clear-cut conclusions. Although 
the ultimate consequences of the theory are sometimes 
in very good agreement with experiments, as Eq. (|62p . 
one should not conclude too fast that RFOT theory is 
correct. In this context, a pertinent line of investigation, 
allowed by numerical simulations, is to more directly test 
the microscopic mechanisms underpinning the deriva- 
tion of the mosaic picture ( Jack and Garrahan. 20051 : 
ICavagna et ai, 20071 ). In particular, recent works have 
shown that the static 'point-to-set' correlation length- 
scale de scribed in Eg. (IBU1) does increase upon super- 
cooling ( Biroli et al., 2008f h Furthermore, first measure- 
ments of the exponent and ip of the mosaic the- 
ory have been obtained with the somewhat surprising 
results 9 ~ 2 > ij} ~ 1 ( Cammarota et al., 2009at 
ICammarota et al., 2009 bt ). At present, direct tests of the 
mosaic picture are quite involved (even conceptually), 
and hence are rare and not yet conclusive. This line 
of research appears nevertheless very promising to estab- 
lish, disprove, o r develop further the mosaic picture, see 
(jCavagna. 20091) for more details. 



but no treatment generically applicable to liquids is yet 
available. 

Second, when applied to three-dimensional liquids, 
RFOT theory relics on several, sometimes distinct, types 
of approximations. For finite dimensional systems, a 



The fact that these approximations become exact 
in several mean-field settings (fully connected models, 
Bethe lattices) suggests that RFOT theory might have 
a status similar to the Curie- Weiss theory for the fer- 
romagnetic systems which docs contain correct elements 
of the real theory. Current research can thus be seen 
as an attempt to understand and describe better iron 
mean-field effects. Going beyond mean-field theory is 
not only technically but also conceptually difficult. For 
instance, Stillinger claimed that the configurational en- 
tropy vanishes only at zero temperatur e, thus suggestin g 
that no entropy crisis can take place (jStillinger. 19881) . 
This was related to his identification of metastable states 
with energy minima, which is now recognized to be 
an incorrect approximation, e ven in mean-field mod- 
els (jBiroli and Monasson. 2000h . where the well-defined 
metastable states are in general made of a large number 
of energy minima, in the spirit of the 'metabasins', some- 
times described in numerical wo rks as a large assembly of 
inherent stuctures ( Heuer. 2008}) . The same criticism ap- 
plies to recent work who similarly claimed having demon- 
strated the abse nce of a glass transit ion in bidimensional 
binary mixtures ( Donev et al., 20071 ). However, a correct 
definition of metast able states beyond mean-fi eld theory 
is still lacking, see ( Biroli and Kurchan. 20 01) for a dis- 
cussion and a first attempt. Thus, going beyond mean- 
field theory is clearly a difficult, but quite exciting and 
important task, from which new results can be expected 



33 



in the future. 

A further source of concern for RFOT theory is the 
fact that quite often the mean-field models from which it 
largely originates seem to behave q uite differently when 
studied in finite dimensions (lAlvarez et al.. 1996[ 
Brangian et al 2002t iMoore and Drossel. 2002 ; 



Moore. 20061 : iMoore and Yeo. 2006^ In fact, there 
does not yet exist a theoretical model in finite dimen- 
sion, for which th e RFOT theory sc en ario can be shown 
to apply, see (ISarlat et al.. 20091 ; iLiers et al., 2010t 
ICastellana et al.. 20101 ) for recent efforts. From a theo- 
retical perspective, such a discovery would be a highly 
decisive step, even if the model were very abstract and 
not obviously connected to experimental systems. 

Thus, there is hope that in the next few years, joint 
theoretical and experimental efforts will drive RFOT the- 
ory into a corner, where its status can be made precise. 
At the time of writing, one can state that RFOT theory 
is still imperfect, but the broadness of its scope and pre- 
dictions, the number of distinct approaches that system- 
atically give back at least some piece of it, the generality 
of the concepts it uses, makes one believe that it contains 
at least some useful seeds to construct a 'final' theory of 
the glass transition-if such a thing exists. 



is far from a trivial assumption since it implicitly uses the 
fact that there exists a temperature below which the ma- 
terial is nearly jammed, so that a description in terms of 
sparse mobility is valid. This approach can thus not self- 
consistcntly capture the microscopic origin of the dynam- 
ical slowing down in supercooled liquids. Instead, it can 
possibly become a very effective and useful description of 
structural relaxation occurring near the glass transition 
(jPalmer et al.. 1984 ). 

At present, it is still unclear whether this main assump- 
tion of mobility conservation is correct, only approxi- 
mate, or whether it plays the central role suggested by 
the facilitation approach. An important point is the fact 
that a large number of theoretical models have been de- 
fined based solely on the idea of kinetic constraints. They 
are called 'Kinetically Constrained Models' (KCMs) and 
they all display a phenomenology which is strikingly con- 
sistent to the on e described in Sec . III. Al for molecu- 
lar glass- formers (jRitort and Sollich. 20031 ). In the last 
decade or so, several KCMs have been defined and stud- 
ied in quite some detail. In turn, this large body of the- 
oretical work has produced new results and predictions, 
and has thus also triggered more research on competing 
theories. 



C. Kinetically constrained models and dynamic facilitation 

1. The physical picture 

Another approach to the glass transition problem 
i s based on the concept of dy namical facilitation, see 
( Chandler and Garrahan. 2010f) for a review. In short, 
'facilitation' captures the physical idea that since viscous 
liquids arc almost solid, mobility is so sparse at any given 
time that any local relaxation event is likely to trigger, 
or 'facilitate' the relaxation of nearby molecules after a 
time which is short compared to the macroscopic relax- 
ation time but large compared to the microscopic one, 
so that mobility can propagate throughout the sample. 
Thus, the focus is less on molecules than on their mobility 
( Fredrickson and Andersen. 1984ah . 

There is no doubt that at least some degree of facil- 
itation is present in nearly jammed materials, but the 
theoretical approach described in this section goes well 
beyond this simple observation and posits that the en- 
tire dynamics of the system is mainly due to facilitation 
effects. This means that typically a mobile region of the 
sample can become unjammed and thus mobile only if 
i t is adjacent to a region whi ch is already unjammed 
( Garrahan and Chandler. 2003?) . This is equivalent to 
postulating that, except for very rare events, mobility 
in a viscous liquid cannot spontaneously arise in an im- 
mobile region in space, nor can it spontaneously decay. 
This is obviously a very strong assumption. 

This constraint is conjectured to become effective 
roughly below Tm ct and to dominate the dynamic al evo- 
lution close to T g ( Garrahan and Chandler. 2003f ). This 



2. Kinetically constrained models (KCMs) 

Over the last 25 years, KCMs have been central to 
the development of the facilitation approach to the glass 
transition. 

The fi rst example we shall prese nt is the Kob- Andersen 
model ( Kob and Andersen. 19931 ). Interestingly, this 
model is quite similar to the lattice glass models de- 
scribed in Sec. IIV.BI The model is again an attempt 
to capture the physics of a hard sphere system, and the 
fact that dynamics becomes slow at high density because 
the environment of each particle is very crowded. The 
Kob-Andersen model is a lattice gas, with occupation 
number on each site i of a regular lattice n; . There is no 
interaction apart from the hard-core constraint, and the 
Hamiltonian is thus trivial: 



H[{m}] = o, 



o, 1. 



(63) 



Geometric frustration is introduced at the level of the ki- 
netic rules, that are defined as constrained local moves. 
Namely, a particle can jump to a nearest neighbor site 
only if (i) that site is empty, to satisfy the hard-core con- 
straint; (ii) the sites occupied before and after the move 
have less than m neighbors, m being an adjustable pa- 
rameter. Kob and Andersen studied the case m — 4 for a 
cubic lattice in d = 3 (m = 6 corresponds to the uncon- 
strained lattice gas), a nd the model displays glas sy dy- 
namics at large density (|Kob and Andersen. 19931) . Such 
kinetically constrained lattice gases have been studied 
in various spatial dimensions, for different values of m, 
for di fferent constraints, or ev en different lattice geome- 
tries (jRitort and Sollich. 20031 ). They can be thought of 
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as models capturing the idea of a 'cage' effect in a strict 
sense, since a particle with a dense neighbor shell cannot 
diffuse. Dynamic facilitation i s thus a direct consequen ce 
of steric effects in this model ( Chaudhuri et ai. 2008al) . 

The crucial difference with lattice glass models is that 
here all configurations are allowed by the Hamiltonian, 
but their kinetics depends on geometry, while in lattice 
glass models kinetic rules ign ore the geometry but no t 
all configurations are allowed ( Biroli and Mezard. 200lh . 
This is actually a deep difference: KCMs assume that 
geometrical constraints act at the level of kinetic rules 
with no thermodynamic counterpart and no reference to 
the interaction which are responsible for them. 

In this lattice gas picture, the connection with glass- 
formers is not obvious because density, rather than tem- 
perature controls the dynamics. Thermal models with 
similar features can in fact be defined by focusing on 
holes rather than particles. This points toward a model 
with a small concentration of mobile regions, which move 
by creation and annihilation. This is actually very rem- 
iniscent of Glarum 'defects' theory where the relaxation 
proceeds via simple diffusion of a l ow concentration of 
independent defects ( Glarum. 19601 ). Using the conju- 
gated ideas of kinetic constraints, facilitation and rare 
defects, Frcdrickson and Andersen defined and studied 
a f amily of kinetic Ising models for t he glass transition 
in ( Fredrickson and Andersen. 1984al) . This last article 
is a seminal paper that opened a whole new research 
avenue, making possible the study of phenomena asso- 
ciated to the glass transition via simple kinetic models. 
The Frcdrickson- Andersen models consist in an assembly 
of non-interacting 'defects', or 'spins', 



/v 



H[{m}] = J^ 



0, f, 



(64) 



where J is an energy scale for creation of mobility, and 
rii = 1 represent the mobility 'state' at site i, whose 
averaged concentration becomes exponentially small at 
low temperature, (n,) « exp(— J/T). As for the Kob- 
Andcrsen lattice gas, the non-trivial ingredient lies in 
the chosen rates for the transition between states. The 
kinetic rules stipulate that a transition at site i can hap- 
pen with a usual Glauber rate, but only if site i is sur- 
rounded by at least k defects (k = corresponds to the 
unconstrained limit). 

As for kinetically constrained lattice gases, these mod- 
els have been studied in different spatial dimensions, on 
different lattices, and using a number of distinct def- 
initions of the kinetic rules, yielding a large number 
of possible glassy b ehaviours ( Ritort and Sollich. 20031 : 
iLeonard et ai, 20071 ). The similarity between those spin 
facilitated models and the kinetically constrained lattice 
gases is striking. Altogether, they now form a large fam- 
ily of models generically called KCMs. 

These models can be divided into several classes. 
'Non-cooperative' models, such as the Fredrickson- 
Andcrsen model for k = 1 (the least constrained 



model) display Arrhcnius dynamic slowing down and 
are thus reminiscent of strong glass-formers. They 
are well-described by simple diffusion of point mo- 
bility defects. 'Cooperative' models display a super- 
Arrhenius dependence of the structural relaxation time. 
This is the case for the Frcdrickson- Andersen model 
with a stronger kinetic constraint, k > f. An- 
other example is the 'East' model where the kinetic 
constraint with k = 1 has a directional character: 



can facilitate the dynamics ( 


Jackie and Eiseinger, 19911; 


iBerthier and Garrahan. 2005 


), which can be rational- 



ized on the basis that displacements in liquids have 
a vectorial character that cou ld extend to facilitation 
( Garrahan and Chandler. 20031) . For an overwhelming 



majority of KCMs, relaxation times diverge only in the 
limit of zero temperature, as even one single defect can 
diffuse and relax the entire system (at least in the sim- 
plest models) , and the defect concentration only vanishes 
at T = 0. However, KCMs can also be defined with 
kinetic rules and geometries for which the existence of 
a finite temperature glass transition can be established 



(Toninelli et ai, 200 



Let us remark that the facilitation approach, and in 
particular KCMs, encode in a new and more microscopic 
way the older, but well-known, concept of free volume. 
Free volume models are among the most widely used 
models to analyze experimental data, especially in poly- 
meric ^s^teins^IlTBvJiave^ be- 
fore (jCohen and Grest. 19821 iDebenedetti. 19961 ). Here, 



the main idea is that dynamic slowing down occurs be- 
cause the 'free volume to move' available to each particle, 
Vf, vanishes at some temperature To as v/ « a(T — To), 
a relation which connects volume to temperature. Sta- 
tistical arguments then relate relaxation timescales to 
free volume assuming that movement is possible if lo- 
cally there is 'enough' free volume available, more than 
a typical value vq. A VFT divergence is then predicted: 



To 



exp 7 



Vf 



exp 



TVo/a 
[T - To}' 



(65) 



where 7 is a numerical factor and /1 = 1. Predictions 
such as Eq. (|65")) are used to justify the wide use of free 
volume approaches, despite the many (justified) criti- 
cisms that have been raised. The physics at work is 
obviously strongly reminiscent of the above description 
of the Kob- Andersen kinetically constrained lattice gas. 
In fact, the analogy is even semi-quantitative in some 
cases, since cooperative KCMs with a finite tempera- 
ture glass transition have been shown rigorously to be 
characterized by divergence of the relaxation time as in 
Eg. (163]) . albeit with a v alue of fj, different from one 



(|Toninelli and Biroli. 20081) . 

As with free volume approaches, it is not exactly clear 
what is meant by 'mobility defects' within KCMs in 
terms of the original interacting system they seek to de- 
scribe, nor how kinetic constraint can truly emerge from 
the unconstrained dynamics of a many-body system. A 
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good news on this front is that at least a proof of prin- 
ciple that kineti c rules can emerge has been obtained 
(|Garrahan. 20021 ). Several examples are available but 
here we only mention the simple case of the bidimen- 
sional plaquette model defined by a Hamiltonian of a 
p-spin type, but in two dimensions on a square lattice of 
linear size L, 



L-l L-l 

H = — J SijSi+ijSij+iSi+ij+i, 
i=i j=i 



(66) 



where Si.j — ±1 is an Ising variable lying at node 
of the lattice. Contrary to KCMs, the Hamiltonian in 
Eq. (|66[) contains genuine interactions, which are no less 
(or no more) physical than p-spin models discussed in 
Sec. lIV.Bl to which it actually strongly resembles. Inter- 
estingly the dynamics of this system is (trivially) mapped 
onto that of a KCM by analyzing its behaviour in terms of 



plaquette variables, pi 



Si ,j Si + 1 ( j iJ i ! j -|_ 



such 



that the Hamiltonian becomes a non-interacting one, 
H = —Jj^i iPi,ji as m Eq. f64|l - More interestingly, the 



analogy also applies to the dynamics ( Garrahan. 20021) . 
The fundamental moves arc spin-flips, but when a single 
spin is flipped the states of the four plaqucttcs surround- 
ing that spin change. Considering the different types 
of moves, one quickly realizes that excited plaquettes, 
Pij = +1, act as sources of mobility, since the energetic 
barriers to spin flips are smaller in those regions. This 
observation allows to identify the excited plaquettes as 
defects, by analogy with KCMs. Similarly to KCMs, also, 
spatially heterogeneous dynamics, diverging lengthscalcs 
accompanying diverging timescales and scaling behaviour 
sufficiently close to T = (see below) ca n be estab- 
lished by further analysis ([Jack et ai, 20051 ). providing 
a simple, but concrete example, of how an interacting 
many-body system might effectively behave as a model 
with kinetic constraints. This type of plaquette models, 
and other similar sp i n models, were origina lly introduced 
( Sct hna et al. , 19911: iLipowskv et ai. 2000? ) ,o show how 
ultra-slow glassy dynamics can emerge because of grow- 
ing free energy barriers. 



3. Diffusing defects, excitation lines and space time bubbles 

On the quantitative side, research on the facilitation 
approach to glassy dynamics has mainly consisted in the 
analysis of the physical behavior of KCMs. Motivation 
for this work largely stems from the observation that 
KCMs behave, at the phenomcnological level, very 
much as real glass-fo rmers, as noted long ago both for 
thermal models dFredrickson and Andersen. 1984at 



Fredrickson and Andersen. 1984bl 
Fredrickson and Brawer. 1986t 



Fredrickson. 1986) 



and constrained lattice g a ses (I Jackie and Kronig. 1994 ; 
iKob and Andersen. 1991 ISellitto et ai, 19971 ). Actu- 
ally, theoretical research on dynamical heterogeneity 
was partly sparked by the numerical observations 



that KCMs display stron g spatial fluctuations of the 
local relaxation rate (iButler and Harrowell. 1991at 



iButler and Harrowell. 1991ri lHarrowell. 19931) . There- 
fore, KCMs provide a simplified context to understand 
glassy phenomena in detail, or at least study one of their 
possible explanations. 

In the following we shall not review all the different 
models an d their physical behavior, as this was done thor- 
oughly in ( Ritort and Sollich. 20031) . Neither shall we re- 
view in great detail the dynamic facilitation approach as 
a theory of the glass transition, as this is the object o f an- 
other recent review ( Chandler and Garrahan. 20101 ). In- 
stead, we shall present a simple physical picture, common 
to all KCMs, which is helpful to grasp the main physical 
behavior of KCMs and, therefore, the main predictions 
and limitations of the facilitation approach. 

A common feature of all KCMs is that their relax- 
ation can be accurately described in terms of the motion 
of sparse defects. In the simplest cases like the non- 
cooperative Fredrickson-Andersen model, these defects 
diffuse, but they might also undergo sub-diffusive mo- 
tion in some cases, like in the directional East model. In 
cooperative models, such as the Kob- Andersen lattice gas 
or the Fredrickson-Andersen model with k > 1, defects 
do not coincide with mobile sites, but can be formed by 
extended clusters moving in a cooperative manner. 

Since local relaxation at a given site occurs when it is 
visited by one such 'defect', explaining structural relax- 
ation is equivalent to explaining the defect dynamics. To 
make the discussion more concrete, let us explain how 
this defect description emerges for the simple case of the 
non-cooperative Fredrickson-Andersen model with k = 1 . 
In that case, a spin Hi can flip only if it has at least one 
spin rij = 1 among its nearest neighbors. At low temper- 
ature, defects diffuse in an activated manner. Typically 
an isolated defect facilitates the excitation of one of the 
neighboring site with rate exp(— J/T), and moves there 
with probability 1/2. Schematically, one gets: 



0100^ 0110^ 0010 



(67) 



Overall, defects perform random walks with a diffusion 
coefficient D ~ cxp( — J/T). Defects can also annihilate 
when they meet and be created from an existing defect. 
A representation of this diffusion, branching and annihi- 
lation process is shown in Fig. [1^] for the one-dimensional 
case. In particular, the central feature of the facilitation 
approach, namely that a defect cannot be created out of 
a completely immobile region, is obvious from this fig- 
ure, in particular at low temperature where defects are 
sparse. 

An approximate but quite good description of KCMs 
consists in focusing on sparse, possibly extended, de- 
fects of density p^, which diffuse anomalously with 
an exponent z and a generalized diffusion coeffi- 
cient D which is temp erature (or density) dependent 
( Toninelli et al.. 2005bt ). The relaxation time of the sys- 
tem can then easily be obtained by computing the time 
over which a finite fraction of the system, say 1/2, has 
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FIG. 16 Space-time representation of the dynamics of mo- 
bility defects in the non-cooperative Fredrickson- Andersen 
model with k = 1 in d = 1. Dynamic facilitation, im- 
posed by the kinetic constraints, implies the existence of ex- 
citation lines which can branch and coalesce. The physics 
changes from a homogeneous description at high tempera- 
tures to a strongly spatially and temporally heterogeneous 
dynamics at low temperatures w hen mobility is sparse. After 
(|Berthier and Garrahan. 2003bl ). 



been visited by at least one defect. The number of dis- 
tinct sites visited by a given defect after a time t is 
by definition equal to {Dt) dF / z , where dp is the frac- 
tal dimension of the walk (or the space dimensionality d 
when dp > d). For a random walk in d = 3, one has 
dp = z = 2, and therefore a number of sites visited per 
defect is simply equal to Dt. The relaxation time r a is 
then given by the equation Pd(DT a ) dF ^ z ~ 1/2. Both 
the defect density pd and the diffusion coefficient vanish 
when decreasing the temperature with expressions that 
are model-dependent. The relaxation times thus generi- 
cally scales as 



*D Pd 



z/dj. 



(68) 



In the case of the non-cooperative Fredrickson- Andersen 
model discussed before (k = 1), one gets r a cx e 2J / T for 
d > 1 and r Q oc e 3,/ / T in d — 1, i.e. an Arrhenius depen- 
dence in all dimensions, similar to the one describing the 
viscosity of strong glass-formers. 

As expected, stronger temperature dependencies are 



obtained when kinetic constraints become more and more 
severe. For the directional East model, for instance, the 
energy barrier to relax a domain of the form 10 ■ • • 01 
containing £ unexcited sites increases with I since the 
leftmost defect must in principle facilitate all sites to its 
right. The minimal energy barrier needed for this pro- 
cess to occur scales as E{£) ~ log(l), and corre sponds 
to a 'hierarchical' path ( Sollich and Evans. 19991) . Since 
the typical distance between defects at equilibrium is 
4 q ~ exp(J/T), one finds that logr Q ~ E(£ cq )/T ~ 
1/T 2 , which is reminiscent of the Bassler law in Eq. 
Even more fragile behaviour is obtained for cooperative 
Fredrickson- Andersen models and Kob- Andersen lattice 
gases. For instance, for the well-studied Fredrickson- 
Andersen model in d = 2 with k = 2 which is the orig i- 



nal model studied in (jFredrickson and Andersen. 1984al) , 
the relaxation time increases as T n ~ exp [exp(c/T)1 , 
with c a numerical c onstant ( Ritort and Sollich. 20031 : 
iToninelli et al, 2005al) . It is interesting to note that 
both dependencies found for the East and k = d = 2 
Fredrickson- Andersen fragile models seem to describe the 
visco sity data obtained in real molecular liquid s quite 
well (|Hecksher et al, 2008t lElmatad et al, 20091) . Actu- 
ally, it has been argued that in order to apply the East 
model results to real liquids, one should use a scaling 
£ cq ~ cxp (J/T — J/T on ). This leads to a modified ver- 
sion of t he Bassler law, which has been discussed and 
tested in ([Elmatad et al . 20091: lElmatad et al, 20~L0l) . 

The space-time representation of the defect diffusion 
dynamics in Fig. [16] can be understood as an alterna- 
tive, very illustrative w ay to represent the dynamical 
behaviour of the system (Garrahan and Chandler. 2002t 



IChandler and Garrahan. 20101) . To the extent that dy- 
namics is facilitated, mobility cannot disappear except if 
it is close to another mobile region. Furthermore, the de- 
fect dynamics implies that mobile regions organize along 
'excitation lines', when represented in a space-time phase 
diagram. Hence, in a system with facilitated dynamics 
space-time is structured with strings of excitations that 
are directed in time, can coalesce or branch but never 
die, nor appear spontaneously. In between these lines lie 
large inactive regions, also called 'bubbles'. The larger 
the spatial extension of the bubbles, the longer it takes 
an excitation line to relax the corresponding domain, and 
the slower is the structural relaxation. Thus, all dynam- 
ical properties can in fact be understood characterizing 
the space-time structure of bubbles and excitations lines. 
One should keep in mind, however, that this is fully 
equivalent to the defect dynamics laid out above. 

To summarize, KCMs can be thought of as 'defect' 
models, where defects can be non-trivial, extended ob- 
jects possibly with non-trivial sub-diffusive or coopera- 
tive dynamics. When temperature decreases both the 
density of defects become small and their dynamics 
slows down dramatically, implying that the overall re- 
laxation slows down in an activated, possibly supcr- 
Arrhenius manner. Additionally, space-time represen- 
tations, as shown in Fig. 1161 suggest the appearance 
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when temperature decreases of large immobility domains 
with broadly distributed spatial and temporal extensions. 
When temperature decreases, these fluctuations are am- 
plified and the corresponding timescales and lcngthscales 
increase rapidly and diverge (for most KCMs) in the limit 
of zero-temperature (or unit density for lattice gases), 
which can thus be cons idered as a dynamic critical point 
(j Whitelam et al, 20041) . as we now discuss. 



4. Main predictions and results 

Within the facilitation approach and, hence, KCMs, 
thermodynamic properties are trivial, and the interesting 
physics is in the dynamics. As a consequence, non-trivial 
predictions and results concern the dynamical behaviour, 
and, more precisely, dynamic heteroge neity, as was re- 
alized soon after their introduction (lHarrowell. 19931: 
iFranz et al, 2002t iGarrahan and Chandler. 20021) Re- 
markably, virtually all the aspects related to dynamic 
heterogeneity mentioned in Sec. IIIII can be investigated, 
and qualitatively or quantitatively rationalized in the lan- 
guage of KCMs. 

Detailed numerical and analytical studies have indeed 
shown that in these systems, non-exponential relaxation 
patterns do stem from a spatial, heterogeneous distribu- 
tion of timescales, direct ly connected to a distribution 
of dynamic lengthscales (jGarrahan and Chandler. 2002 ; 



Ritort and Sollich. 20031: 



Toninelli et al, 2005al 



Whitelam et al, 2005 



Jack et al, 20 06a') 



Pan et al, 2005 



A decoupling phenomenon between the structural re- 
laxation time and the self-diffusion coefficient naturally 
appears in KCMs and can be shown to be a very direct, 
qua ntifiable, consequ ence of the dynamic heterogene- 
ity (|Jung et al, 2004fh In fact, two distinct timescales 
emerge in the facilitation/KCMs description which both 
characterize the dynamics. One is the persistence time, 
t p , which corresponds to the typical time a tracer parti- 
cle has to wait to start moving, given an arbitrary start 
time for observation. The second timescale is the ex- 
change time, t Xl which is the typical time between ele- 
mentary moves. The average value of t p and a 2 /(t x ) are 
respectively estimates of the relaxation time r a and the 
self-diffusion coefficient D. with a the typical size of par- 
ticle jumps. As a consequence one finds Dr a ~ (t p }/(t x ). 
In some KCMs, this ratio increases substantially when 
decreasing the temperature, thus these models exhibit a 
substantial decoupling phenomenon. 

This is due to the fact that bubbles in space-time 
typically persist for the time scale of structural re- 
laxation, while particles can diffuse faster by 'surf- 
ing' on the excitation lines that surround these bub- 
bles. Another way to understand the decoupling 
is to express (t p ) and (t x ) as distinct moments of 
the distribution of waiting time between jumps per- 
formed by a tracer particle, which do not coincide 
when the distribution broadens with decreasing temper- 



atures us and Kivelson. 19951 : iHedees et al, 20071: 
iHeuer. 20081) 

The above description of self-diffusion in fluctuating 
mobility fields such as shown in Fig. [TH] has several 
non-trivial consequences such as for instance a strongly 
non-Fickian relation between timescales and lengthscales 
at distances smaller than a crossover lengthscale which 
grows when temperature dec reases and decoupling in- 
creases ( Berthier et al, 2005 rj ), as observed in num erical 
simulations of supercooled liquids ( Berthier. 2004] ) . 

Another useful aspect of KCMs is that multi-point 
susceptibilities, multi-point spatial correlation functions 
such as the ones defined in Eqs. (|2"2"|) and (|20[) can be 
studied in much greater detail than in molecular systems, 
to the point that scaling relations between timescales, 
lcngthscales . and dynamic su sc eptibilities can be es- 
tablished (iFranz et al, 2002t I Whitelam et al, 2004 : 
Toninelli et al, 2005ri iPan et al, 2005 ; 



Chandler et al, 2006t iBerthier et al. 2007bl ) 

The type of scaling behaviour with the temperature 
or the density depends on the details of the model at 
hand, in particular whether it has or does not have a 
transition. However, a unified qualitative description is 
still possible in terms of th e defect motion descri ption 
given in the previous section ( Toninelli et al. 2005bl ) . On 
timescales small compared to r a , all sites that have been 
visited by the same defect are dynamically correlated. 
As a consequence, the dynamical correlation length £i(t) 
increases as ^(t) ~ {Dt) 1 ^ and the corresponding four- 
point susceptibility, Xi(t) — Pd{Dt) 2dp / z , which is the 
square of the volume of the regions visited by the same 
defect times the defect density. When t is comparable to 
the relaxation time r Q , one finds Xi{ T a) — 1/ Pd which, 
as expected, increases when decreasing the tempera- 
ture and increasing the density. On longer timescales, 
Xi{^> T a ) decreases with t, because a site can now be 
visited by more than one defect, and dynamical corre- 
lations are progressively erased. As a direct example 
of this behaviour, we show in Fig. [T7] the time and 
temperature evolution of Xi{t) f° r f ne non-cooperative 
(k = 1) Fredrickson- Anders en model in three dimen- 
sions (I Whitelam et al, 20051). where the above scaling 
relations hold (jToninelli et al, 2005bl) . 

In summary, relaxation timescales and dynamic 
lcngthscales are found to diverge with well-define d 
critical laws (| Whitelam et al. 2004 ; I Jack et al, 2006af) . 
which however are model dependent. The discovery of 
such 'dynamic criticality' has proven useful, since it ac- 
tually triggered several works aimed at computing similar 
laws within competing theories. It also implies the possi- 
bility that some universal behaviour might truly emerge 
in the physics of supercooled liquids, precisely of the type 
observed numerically in Fig. 1 1 01 and experimentally in 

Fig.rnj 

As suggested by the space-time representations shown 
in Fig. 1161 dynamical heterogeneity studies are mostly 
concerned with distributions of dynamical quantities, 
most often through their variance, recall the definition 
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FIG. 17 Time evolution of x^it) f° r different temperatures for 
the non-cooperative (k — 1) Fredrickson- Andersen model in 
three dimensions. Temperature decreases from left to right. 
The peak height increases at low T as x-i{ T a) ~ 1/pd ~ 
exp(J/T), while the relaxation time scales as r a ~ exp(2J/T), 
such that x-i{ r a) ~ Tq^ 2 , as shown with a dotted line, while 
the approach to the peak is diffusive, X-i(i) ~ t 2 (dashed line). 



of X4 in Eq. (|22p . However, there is potentially useful 
information also in the overall shape of the distributions. 
It seems reasonable that regions where the correlation 
function is large possess rather stable structure at the 
molecular level, while regions where it is small corre- 
spond to relatively unstable local structure. To identify 
such trajectories over an observation time t bs, it is use- 
ful to define a measu re of dynamical activity, for example 
([Hedges et al. 20091 ) 



JV tobs/A* 



K(N,t ohs )=J2 E l*i(*i)-*i(*i-A*)| S 



(69) 



where rj(f) is the position of particle i at time t and 
the tj — jAt arc equally spaced times. For large N and 
t bs, the distribution of K becomes sharply peaked about 
its average, (K). In general, for large N and t obs , one 
expects the distribution of K to have the form 



P(K) ~ exp[-Nt obs f(K/Nt obs )}, 



(70) 



where the functio n f(k) resembles a fr ee energy density. 

In some KCM (jMerolle et al, 2005T ). the distribution 
P(K) has a characteristic shape, skewed towards small 
activity, with an apparently exponential tail. Further, 
on estimating f(k) from this plot, there is a range of 
K over which f(k) is non-convex (that is, f"(K) < 0). 
Within such a framework, non-convexity of f(K) has a 
direct interpret ation as a 'dynamical phase transition' 
in the system ( Garrahan et al.. 2009T) . The existence 
of these phase transitions ha s been proven in simple 
models ( Garrahan et al. 20091 ) and numerical results for 
Lennard-Jones model liquids are also consistent w ith the 
existence of such a transition ( Hedges et al.. 2009th which 
also exists in mea n-field models where RFO T theory of 
Sec. lIV.Bl applies ( Jack and Garrahan. 2010l ). This leads 



to the hypothesis that the nature of the dynamically het- 
erogeneous fluid state could be generically interpreted in 
terms of coexistence between an active liquid and inac- 
tive 'ideal glass' states, with seems to suggest intriguing 
connections to the mean-field picture based on the exis- 
tence of numerous metastable states. 



5. Current status of the facilitation approach 

We have conveyed the idea that despite the large num- 
ber of distinct KCMs capturing the idea of dynamic fa- 
cilitation, these models arc characterized by a common 
physics, based on the propagation of some rare 'defects'. 
Quantitatively, however, different models might behave 
quite differently. The variety of physical outcomes ob- 
tained from the variety of models and kinetic rules has 
proven useful, as it has widened the spectrum of all 
the possible behaviours that real material could display. 
Thus, in the recent effort aiming at devising tools to 
quantify and understand dynamic heterogeneity in glass- 
formers, as described in Sec. IIII1 KCMs have played a 
leading role. 

As far as building a unique theory for the glass tran- 
sition is concerned, this variety is also a source of worry 
since it is difficult to obtain a unified description: many 
models do not have a glass transition at finite tempera- 
ture, but some have; some models seem to describe strong 
glass-formers, other are better for fragile materials; quan- 
titative results about the behaviour of high order corre- 
lation functions or decoupling phenomena quite strongly 
depend on the chosen model, the value of an upper crit- 
ical dimension qualifying the strength of fluctuations is 
also different from one model to another. At the moment, 
it is not clear whether or why one of these models should 
be preferred to the others, since there exists no system- 
atic coarse-graining procedure to represent a molecular 
model with continuous degrees of freedom in terms of a 
lattice model with no interactions and kinetic constraint s 
(jVogel and Glotzer. 200l iDownton and Kennett. 20071 ). 
Thus, the choice of a 'reference' model to fit real data 
is an ambiguous issue, which would clearly deserve more 
work. 

It is interesting to remark that the theoretical sta- 
tus of the dynamic facilitation approach is almost oppo- 
site to that of the RFOT theory, since there are several 
well-defined finite dimensional theoretical models captur- 
ing the idea of dynamic facilitation which can be thor- 
oughly understood, but none of them can be derived 
from (even approximate) microscopic calculations. By 
contrast, RFOT theory is supported by several micro- 
scopic approaches, but theory is not confirmed within the 
framework of finite-dimensional models. This compari- 
son raises the following question: is it possible to obtain 
more microscopic theoretical insights from the dynamic 
facilitation approach beyond qualitative or scaling con- 
siderations? This absence of a microscopic derivation, or 
at least of some empirical procedure to back these ideas 
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makes the approach prone to criticisms. 

The idea that all thermodynamic aspects are unim- 
portant implies that the facilitation approach has little 
to say about the thermodynamic behavior of glass- 
formers close to T 3 . A possible coincidence between 
VFT and Kauzmann temperatures, To and Tk, is not 
expected (except if they are both zero), nor can the 
dynamics be deeply connected to thermodynamics, 
as in Ada m-Gibbs relations . Some view this as a 
deep flaw (jBiroli et al. 20051 ). while within the KCM 
approach, one is forced to disregard the significance 
of thermodynamics (IChandler and Garrahan. 20051: 
IChandler and Garrahan. 20101 ). This is certainly the 
point where KCMs and RFOT theory differ most 
evidently. Even though the dynamics of KCMs shares 
similarities with syste ms characterized with a com- 
plex energy landscape feerthier and Garrahan. 2003at 



IWhitelam and Garrahan. 2004 ) and KCMs even 
show a MCT-Iike tra nsition on Bethe lattices 
(jToninelli et al. 2005at ISellitto et al.. 2005h . their 
thcrmodynamical behaviours are widely different 
from RFOT theory, as was recently highlighted in 
( Jack and Garrahan. 2005h . It will be interesting to fol- 
low how the facilitation approach will handle the recent 
surge of interest in the definition and measurements of 
static 'point-to-set' correlations, which are nonexistent 
in KCMs in the defect formulation. 

A second central assumption made in KCMs is that 
mobility can neither be spontaneously created nor de- 
stroyed unless adjacent regions are mobile, or, at 
least, that events violating this constraint are rare and 
become rarer when approach in g the glass transitio n 
(|Garrahan and Chandler. 20031: lElmatad et al.. 20101) . 
Without these assumptions, KCMs become trivial mod- 
els where death and birth rates for mobility ultimately 
control the dynamics, and the features described in pre- 
vious sections become irrelevant. The fact that some 
degree of facilitation is present in the dynamics is very 
reasonable and partially p r oven by numerical simulations 
dVogel and Glotzer. 2004 iDownton and Kennett. 20071 
ICandelier et al.. 2009bt IChandler and Garrahan. 2010}) . 
However, the facilitation assumption goes much further 
and states that almost nothing else is possible and that 
rare events, violating the kinetic constraints, become 
rarer approaching T g , so that the mobility is effectively 
nearly conserved. There is at present no data support- 
ing this assumption. A recent experimental analysis of a 
granular system close to its glass transition actually sug- 
gests the opposite behaviour that mobility becomes less 
conserved and facilitation play s a decreasing role whe n 
the transition is approached ( Candelier et al. 2009al ). 
Certainly, this is an important issue to be addressed 
in the future in order to validate or disprove this the- 
oretical approach to the glass transition. It would 
be interesting in parti cular to repeat the analysis of 
(|Candelier et al.. 2009a[) in a model of supercooled liq- 
uid. 



D. Geometric frustration, avoided criticality, and Coulomb 
frustrated theories 

1. Physical picture and simple models 

In all of the above models, 'real space' was present 
in the sense that special attention was paid to different 
lcngthscales characterizing the physics of the models that 
were discussed. However, apart from the 'packing mod- 
els' with hard-core interactions, no or very little attention 
was paid to the geometric structure of local arrangements 
in molecular liquids close to a glass transition. This is 
generally justified using concepts such as 'universality' or 
'simplicity', meaning that one studies complex phenom- 
ena using simple models, a typical statistical mechan- 
ics perspective. However, important questions remain: 
what is the liquid structure within mosaic states? How 
do different states differ? What is the geometric origin 
of the defects invoked in KCMs? Are they similar to de- 
fects (disclinations, dislocations, vacancies, etc.) found 
in crystalline materials? 

There exists a line of research which attempts to pro- 
vide answers to these questions ( Tarius et al.. 200oT ). It 
makes heavy use of the concept of geometric frustration. 
Broadly speaking, frustration refers to the impossibility 
of simultaneously minimizing all the i nteraction term s 
in the energy function of the system ( Toulouse. 19771) . 
Frustration might arise from quenched disorder (as in 
the spin glass models described above), but liquids have 
no quenched randomness. In that case, frustration has a 
purely geometric origin. It is attributed to a competition 
between a short-range tendency to extend spatially a 'lo- 
cally preferred order', and global constraints that prevent 
the periodic tiling of space with this local structure. 

This can be illustrated by considering once more 
the packing problem of sp heres in three dimen- 
sions ( Sadoc and Mosseri. 19991 ). In that case, the locally 
preferred cluster of spheres is an icosahedron. However, 
the five-fold rotational symmetry characteristic of icosa- 
hedral order is not compatible with translational sym- 
metry, and formation of a periodic icosahedral crystal is 
impossible ( Frank. 19521) . By contrast, disks on a two- 
dimensional plane arrange locally as a regular hexagon, 
with one atom at the center and six neighbours at the 
vertices. If periodically repeated, this structure can then 
form a triangular lattice that can fill space with no influ- 
ence of geometric frustration. 

The geometric frustration that affects spheres in three 
dimen sional Euclide an space can be relieved in curved 
space ([Nelson. 20021) . This corresponds for instance to 
studying particles on a sphere, or on the hyperbolic 
plane (a surface of constant negative curvature). Indeed 
by changing the metrics and topology of the underlying 
space it may become possible for the local order to ex- 
tend over larger lengthscales. One can fruitfully exploit 
this idea in two ways. A first possibility is to start from 
a curved space carefully chosen such that geometric frus- 
tration is entirely absent for the considered system. The 
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structure of the system minimizing the energy can then 
be determined, and serves as a useful reference state. 
Changing the space curvature to go back to the physical 
Euclidean space then generates topological defects that 
disturb the initially perfect order. Detailed analysis along 
those lines showed in particular that a sphere packing 
possesses, in Euclidian space, topological defects (mainly 
disclination lines) , as the result of forcing the ideal icosa- 
hcdral ordering into a 'flat' space. Thus the relevant 
topological defects in that case are one-dimensional ob- 
jects forming a disordered network, presumably having 
a very complex dynamics. Nelson and coworkers have 
developed a solid theoretical framework based on this 
picture to suggest that the slowing down of supercooled 
liquids is d ue to the slow wandering of these topologi- 
cal defects (|Nelson. 20021 ). but their treatment remains 
so complex that few quantitative, explicit results have 
been obtained, in particular concer ning the dyna mical 
behaviour of the frustrated systems ( Nelson. 1983|) . 

The picture of sphere packing disrupted by geomet- 
ric frustration is appealing and provides handles to at- 
tack the problem of glass formation from an atom- 
istic perspective. Furthermore, the idea of uniform 
frustration can be incorporated into simple models al- 
lowing for a more abstract approach in terms of sim- 
ple statis tical mechanics ide as and scaling-type of ap- 
proaches ( Tarius et al.. 20051 ). To build such models, one 
must be able to identify, then capture, the physics of ge- 
ometric frustration. Using the concept of an ideal long- 
range ordering in a system of size L in a curved space, 
which is then strained back in the three-dimensiona l Eu- 
clidian space, Kivelson et al. ( Kivelson et al.. 19951 ) sug- 



gest that the corresponding free energy should scale as 



F(L, T) = <j{T)L 2 - (f>(T)L 3 + s{T)L 5 



(71) 



In this expression, the first two terms express the ten- 
dency of growing local preferred order and they repre- 
sent respectively the energy cost of having an interface 
between two phases and a bulk free energy gain inside 
the domain. It is assumed that without the last term 
long range order would sets in at T = T*. Geometric 
frustration is encoded in the third term which represents 
the strain free energy resulting from the frustration. This 
last term is responsible for the fact that the transition is 
avoided. The remarkable feature of Eq. (|TTj) is the super- 
extensive scaling of the energy cost due to frustration 
which opposes the growth of local order; in dimension d 
it would scale as L 2+d . As a consequence, when L is large, 
the last term becomes dominant and prevents thus the 
extension of order to the entire space. Thus, the system 
is broken up in a 'patchwork' of locally ordered domains 
separated by domains made of topological defects — hence 
the name of frustration limited domains. Furthermore, 
minimizing the free energy per unit volume, F/L 3 , one 
finds that the characteristic linear size L* of the patches 
scales as (er/s) 1 / 3 . Since a increases below T*, one finds 
that the characteristic size of the patchwork increases be- 
low T* too. 



Turning to the dynamics, Kivelson et al. further ar- 
gue that dynamics of the system involves restructuring 
of these domains in a thermally activated manner, us- 
ing arguments similar to the ones used within the mosaic 
picture of RFOT theory. The typical energy barrier is 
given by a high temperature constant plus a second term 
a(L*) 2 , increases below T*, which means that the as- 
sumption ip = 2 is made from the beginning, see Eq. (|61|) . 
Using more refined but still heuristic arguments, Kivel- 
son et al. argued that L* grows as (1 - T/T*f /K 1 ' 2 
where K is an adimensional parameter measuring the 
strength of frustration, and v is the exponent governing 
the growth of the correlation length of the unfrustrated 
transition (K = 0) . The corresponding prediction for the 
energy barrier is 



A(T < T*) = A ; 



AT* 



1 



T 

T* 
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where A> is the high temperature value of the bar- 
rier and A is a positive constant. Furthermore, it has 
also been argued that barrier fluctuations lead to typ- 
ical glassy features such as broad distribution of relax- 
ation times or spatially heterogeneous dynamics, which 
can thus be discussed in a phenomenological manner 
( Tarius and Kivelson. 199^ . We refer the reader to 
(IViot et ai. 20001) for more details on these aspects. An 
important prediction obtained from scaling considera- 
tions is the variation of the glass fragility with frustration: 
in this approach, larger frustration means smaller do- 
main sizes, and therefore less collective relaxation yield- 
ing smaller fragility. 

A complementary route consists in the analysis of rel- 
atively simple, finite dimensional statistical models that 
supposedly retain the basic physical elements of Eq. (|7T|) . 
We consider first an interacting spin model where the 
magnetization is meant to represent the 'preferred lo- 
cal order', nearest neighbour ferromagnetic interactions 
the tendency to local ordering, and longer-ranged anti- 
ferromagnetic interactions the opposite effect of the frus- 
tration. The following Hamiltonian possesses, in three 
dimensions, these minimal ingredients: 



(i.j) 
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where the spin Si occupies the site i at position Xj. It 
can be shown that the long-range Coulombic interac- 
tion plays a role analog to the super-extensive free en- 
ergy in Eq. (|TTj) . Of course, such Coulomb frustrated 
spin models can be for Ising or multi-component spins, 
for hard or soft spins, or gener alized to different spatial 
dimensions ( Tarius et al.. 2005T ). For theoretical studies 
it is convenient to study field-theoretical versions of the 
Hamiltonian ([731 : 
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where the ferromagnetic and Coulombic terms are easily 
recognized, and the magnetization </>(x) is now a field. 

It is interesting to notice that the competition of inter- 
actions acting at different lengthscales in fact describes 
a much larger body of problems, including the physics 
of diblock copolymers, magnetic multi-layer compounds, 
Raylcigh-Bcnard convection, doped Mott insulators, etc. 
A well-studied, related field theory is given by: 



H = d d 



1 7; 
-<Xx) (r + fc - 2 (V 2 + fc 2 ) 2 ) 0(x) + i 0(x) 4 

_ (75)" 

Although this field theory was mainly considered 
from the p oint of view of diblock copolymers 
Leibler. 1980T) and Ra yleigh-Benard convection 



Swift and Hohenberg. 1977^ . it was explicitly con- 



sidered in the conte xt of the glass problem in 
dSchmalian et ai. 20031: Geissler and Reichman. 20031: 



iGeissler and Reichman. 20041 ) . Dimensional analysis 
shows that, in that case, the strength of frustration is 
related to ko through ko ~ Q 1 ^ 4 . 

This family of statistical models yields a rich physical 
behaviour, with a complex phase diagram and dynami- 
cal behaviour. In the absence of frustration, Q = 0, they 
generally undergo a second order phase transition from 
a paramagnetic to a ferromagnetic phase, which should 
be viewed, in this context, as the analog of the ordering 
transition occurring in the curved space relieving the geo- 
metric frustration of the sphere packing. That transition 
occurs at some finite temperature, T C {Q = 0) = T*. The 
presence of frustration, Q > 0, generally has dramatic 
effects whose details depend, however, on the s tudied 
model (jChaves et ai. 19961 iNussinov et ai. 19991) . The 
transition can be entirely suppressed, T C (Q > 0) = 0, 
or severely depressed as soon as Q > yielding a gen- 
uine discontinuity when Q — > 0. By decreasing T at 
small but finite Q, the system gets close to, but nar- 
rowly 'avoids', the critical point at T*. This situation 
occurs for instance in the 0(N — > 00) and spherical ver- 
sions of the model. In the more canonical Ising case, the 
situation is different since the second order transition oc- 
curring for Q = becomes first order at finite Q, but 
there is no discontinuity of the transition temperature as 
Q —> ( Brazovskii. 19751) . At finite Q, the ordered phase 
is a spatially modulated phase (stripes in d = 2, lamellae 
in d = 3, etc.). Interestingly, the limit of stability of the 
paramagnetic phase (the spinodal) is depressed down to 
T = 0, and it is thus (in principle) possible to supercool 
the first order transition and study the disordered phase 
down to low temperatures. 

In all cases, therefore, there exists a temperature 
regime below the 'avoided critical point', T < T*, where 
the system is disordered and the dynamics is poten- 
tially affected by the presence of the frustration. De- 
tailed Monte Carlo studies show that, for a given frus- 
tration strength Q, the dynamics slows down consider- 
ably when T decreases below T* dGrousson et ai. 200 ll 
iGrousson et ai, 20021 iTarius et al.. 20051) . Grousson et 
al. quantified dynamics using spin-spin autocorrela- 



tion functions, C(t) = (Si(t) ■ Si(0)), in a variety of 
models (Ising, 5-state spins, and XY models) in three- 
dimensions. They show that the time decay of C(t) 
gets non-exponential at low T and extract a relax- 
ation timescale that shows activated, super-Arrhenius 
behaviour. Remarkably, they also find that the fragility 
decreases when frustrati on increases, in agree ment with 
the scaling approach of ()Kivelson et al.. 19951) . and sim- 
ulation s of Lennard- Jones models on the hyperbolic 
plane ( Sausset et al. 20081) . suggesting that Coulomb 



frustrated spin models do indeed capture an essential 
part of the physical behaviour of molecular glass-formers. 

This op timistic view was contradict ed in more re- 
cent work ( Geissler and Reichman. 20041 ). where numer- 
ical simulations of the Langevin dynamics of the field- 
theory of Eq. (|75[) were reported. The numerical analy- 
sis yields results for spin autocorrelation functions C(t) 
in agreement with those of Grousson et ai. However, 
Geissler and Reichman also measured the time decay 
of fluctuations of the Fourier components of the field, 
Ck(t) = (cf>(k, t)(j>{— k, 0)), and found exponential decay 
for all Fourier modes, Ck(i) ~ exp(— t/r(k)). There- 
fore they interpret the stretched exponential relaxation 
reported by Grousson et al. as being the result of a triv- 
ial superposition effect, C(t) cx J d d kC^(t). Even more 
striking was their finding that the relaxation times r(fc) 
were very accurately predicted by a simple dynamical 
Hartree approximation mirrorin g the classic static treat- 
ment of the field-theory ([75) in (|Brazovskii. 19751 ). This 



finding implies that the dynamical slowing down detected 
numerically results from the strong temperature depen- 
dence of static correlations which will eventually yield 
the system to undergo a first-order transition towards a 
modulated phase. This scenario bears little resemblance 
with the physics of supercooled liquids where dynamics 
is not obviously driven by simple static correlations, as 
we emphasized several times in this paper. 



2. Current status of the Frustration Limited Domains theory 

We believe that the generic idea of geometric frustra- 
tion is a very appealing one, because it directly addresses 
the physics in terms of the 'real space' at the molecular 
level. Moreover, it seems to yield quite naturally the 
idea that the system organizes itself, at low temperature 
and finite frustration level into some 'mosaic' of domains 
corresponding to some local order whose size increases, 
but does not diverge, when T decreases. Tarjus, Kivelson 
and co-workers clearly demonstrate that such an orga- 
nization of supercooled liquids into mesoscopic domains 
allows one to understand most of t he fundamental phe - 
nomena occurring in glass-formers (jTarius et al.. 20051) . 
Our use of the word 'mosaic', as in RFOT theory, is 
not accidental. It was demon str ated in several works 
dSchmalian and Wolvnes. 20001 IWestfahl et al.. 20011 
INussinov. 20oI ~ that, when treated as systems with 
one-step of replica symmetry breaking as discussed in 
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Sec. IIV.B.61 Coulomb frustrated systems such as in 
Eq. (|74| undergo a Kauzmann transition very similar 
to the one found in mean-field spin glasses. This is 
remarkable and suggests that the two theories may 
not be so different as they appear at first sight. The 
presence of structured domains also connects to ideas 
such as coopcrativity, dynamic heterogeneity and spatial 
fluctuations, that directly explains, at least qualitatively, 
non-exponential relaxation, decoupling phenomena or 
super-Arrhenius increase of the viscosity. The variation 
of fragility with the strength of the geometric frustration 
is perhaps the most striking and original prediction 
stemming from this approach. 

As for the RFOT mosaic picture, direct confirmations 
of this scenario from molecular simulations are difficult 
to obtain. Two independent routes have been followed 
recently. On the one hand there have been studies focus- 
ing on the role of locally preferred structures. The iden- 
tification of th e locally preferred struct ure is in general 
quite difficult ( Mossa and Tarius. 20061) . except when it 
can be linked to formation of local icosahedral order. 
Coslovich and Pastore have found by numerical simu- 
lations of several models of glass-forming liquids that 
the increase of icosahedral order upon super-cooling is 
more rapid and more pro nounced for liquids character - 
ized by a higher fragility ( Coslovich and Pastore. 20071) . 
This suc cess of the theory is tempered by recent sim- 
ulations ( Tanaka et al.. 2010t ) . that suggest instead us- 
ing different model systems, including three dimensional 
polydisperse hard spheres, that the increasing local order 
is the hexagonal bond orientational order (characteristic 
of the crystal state) a nd not the icosahedra l one. In yet 
another recent paper (jPedersen et al. 2010() . local order 
is again linked with slow structural relaxation but the in- 
terplay between local order and frustrated crystallization 
appears much more complicated than originally surmised 
in the frustration limited domain approach. Thus, we are 
still far from understanding whether the basic content of 
this approach are wrong, correct, or too simplistic to de- 
scribe supercooled liquids, and more work is needed to 
resolve this issue. 

The other route that has been followed consisted in an- 
alyzing a situation very similar to the one FLD advocates 
to happen for glass-formers, such as a two-dimensional 
monoat omic Lennard- Jones syst e m on the hyperbolic 
plane dSausset and Tarius. 20071 ISausset et al.. 20081 
ISausset and Tarius. 20081 ). The hyperbolic plane is a 
two-dimensional surface of constant negative curvature, 
and its relevanc e had been suggested earlier by Nelson 
and co-workers ( Rubistein and Nelson. 19831 ). but these 
numerical studies are te chnically not straightforward 
( Sausset and Tarius. 20071 ). In the hyperbolic plane, the 
hexagonal crystal order is frustrated. In this setting the 
crystal order represents mutatis mutandis the locally pre- 
ferred 'glass' order that is conjectured to be frustrated 
for three dimensional glass-formers. This provides a 
concrete example of a particle model characterized by 
an avoided phase transition. Furthermore, by chang- 



ing the curvature of the space one can test the predic- 
tions on the fragility dependence on frustration, which 
is related to the curvature of the plane. Sausset et al. 
confirmed numerically that by curving space, they can 
prevent crystal ordering, as expected, and that the low 
temperature dynamical behaviour of the liquid they ob- 
tain is very sluggish and has a behaviour very similar to 
glass-formers: time correlation functions decay in a two- 
step manner with a final decay which is not exponential, 
the d ynamics is spatially heter ogeneous, decoupling oc- 
curs (jSausset and Tarius. 2008T ). Interestingly, they also 
find that the fragility of the glass -formers they obtai n 
decreases when frustration grows (jSausset et al.. 20081 ). 
The physical explanation is that increasing the frustra- 
tion also increases the density of topological defects, so 
that the putative collective behaviour responsible for 
super-Arrhenius relaxation at small frustration becomes 
shorter-ranged when frustration increases. 

As a conclusion, this approach to the glass transi- 
tion is based on a very appealing real space descrip- 
tion. As others, it is compatible with several experi- 
mental results, and it has passed non-trivial benchmark 
tests. A possible criticism, which we also formulated in 
the case of the facilitation approach, is that no quan- 
titative or microscopic results beyond scaling have been 
obtained, and it is not even clear how they could actually 
be derived. In fact, approximate analytical approaches to 
FLD, that could have provided quantitative predictions, 
have found back results consistent with RFOT theory. 
In order to further test the FLD approach, it would be 
interesting to have more detailed, quantitative, testable 
predictions regarding the behaviour of thermodynamic 
(e.g., specific heat) and dynamic (e.g. four-point dy- 
na mic susceptibility) observab les from this perspective, 
see ( Sausset and Tarius. 2010l ) for recent work in this di- 
rection. 



V. OFF-EQUILIBRIUM DYNAMICS: AGING AND 
RHEOLOGY 

A. Why aging? Phenomenology and simple models 

We have dedicated most of the above discussion to 
physics taking place when approaching the glass tran- 
sition at thermal equilibrium. We discussed a rich phe- 
nomenology and serious challenges for both our numeri- 
cal and analytical capabilities to account for these phe- 
nomena. Most people, however, focus on the properties 
of glasses, i.e. below the glass transition, so deep in the 
glass phase that the material seems to be frozen forever in 
an arrested amorphous state, endowed with enough me- 
chanical stability for a glass to retain, say, the liquid it 
contains (preferentially a nice red wine either from France 
or Italy). One could naively think that in this state of 
matter the dynamical evolution is arrested. Is it true? 
The answer is clearly 'no'. There is still life (and physics) 
below the glass transition. We recall that for molecular 
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glasses, T g is defined as the temperature below which 
relaxation is too slow to occur within a given experimen- 
tal timescale. Much below T g , therefore, the equilibrium 
relaxation timescale is so astronomically large that ther- 
mal equilibrium is out of reach. One enters therefore the 
realm of off-equilibrium dynamics. A full physical un- 
derstanding of the non-equilibrium glassy state rema ins 
a central challenge (jYoung. 1998t iBarrat et at. 20031) . 

An immediate consequence of studying materials in 
a time window smaller than equilibrium relaxation 
timcscalcs is that the system can, in principle, remember 
its complete history, a most unwanted experimental sit- 
uation since all details of the experimental protocol may 
then matter. The simplest protocol to st udy aging phe - 
nomena in the glass phase is quite brutal ( Struik. 1978() : 
take a system equilibrated above the glass transition 
and suddenly quench it to a low temperature. The sys- 
tem then tries to slowly reach thermal equilibrium, even 
though it has no hope to ever get there. Aging means 
that the system never forgets the time £ w spent in the 
low temperature phase, its 'age', and that any measure- 
ment started at time t w might have an outcome that ex- 
plicitely depend on the value of t w , unlike the situation 
at thermal equilibrium. 

This implies in particular that any physical property 
of the glass becomes an age-dependent quantity in an 
aging protocol, and more generally dependent on how 
the glass was prepared. One can easily imagine using 
this property to tune mechanical or optical characteris- 
tics of a material by simply changing the way it is pre- 
pared, like how fast it is cooled to the glassy state. A 
very striking experimental realization of this idea was re- 
cently pr ovided for organic gla sses prepared in two differ- 
ent ways ( Swallen et ai, 2007t) . Ediger and collaborators 
have compared the properties of glasses prepared in a 
canonical way (slow cooling of bulk samples), to the ones 
of samples grown using slow vapor deposition at about 
50 K below the glass transition. They find that the latter 
samples are much more 'stable', in the sense that they 
behave as canonically prepared samples with a very large 
effective aging time (40 years for some sample). 

Coming back to the simplest situation of a sudden 
quench to low temperature, it is found that one-time 
physical observables such as density or energy evolve very 
slowly with the age of the sample. In polymer glasses 
for inst ance, the volu me of the sample slowly decreases 
with i w ( Struik. 19781 ). Power laws with small exponents, 
v{t w ) ~ Woo + (to/t w ) a , or even logarithmic relaxations, 
v(t w ) ~ vq — v\ log(twAo), are frequently reported. In 
some cases, the time evolution of these observables is so 
slow that aging behaviour is not obviously revealed by 
their study and the system might superficially appear to 
have reached equilibrium. 

In order to show that the system never equilibrates, 
two-time quantities, such as density-density or dipole- 
dipole correlation functions, are much more useful. A 
typical example is presented in Fig. [T5] where the self- 
part of the intermediate function in Eq. © is shown for a 




FIG. 18 Aging dynamics in a Monte Carlo study of a 
Lennard- Jones glass-forming liquid at low temperature. The 
system is quenched at time t w = to a low temperature which 
is then kept constant. Two-time self-intermediate scatter- 
ing functions are then measured for 20 logarithmically spaced 
waiting times i w from r, w = to i w = 10 5 (from left to right). 
The relaxation becomes slower when i w increases: the system 
'ages'. 



Lennard- Jones molecular liquid at low temperature. Im- 
mediately after the quench, the system exhibits a rela- 
tively fast relaxation: particles still move substantially. 
However, when the age of the system increases, dynam- 
ics slows down and relaxation becomes much slower. The 
relaxation then separates into two well-defined 'time sec- 
tors'. For short time differences, t — t w <C t w , which 
corresponds to short-time dynamics of the amorphous 
structure, correlation functions at different t w superpose. 
However, for large time differences, t — t w ^>t w , different 
curves relax at different rates, implying that structural 
relaxation becomes slower when the system gets older. 
Eventually, when t w becomes very large, the relaxation 
becomes too slow to be followed in the considered time 
window and the system seems frozen on that particular 
timescale. For practical purpose, it has now become a 
glass. 

A striking feature conveyed by these data is that an 
aging system not only remains out-of-equilibrium for- 
ever, but its typical decorrelation or relaxation time 3 is 
in fact directly set by its own age i w , which also sepa- 
rates short-time, equilibrated dynamics from long-time, 
aging relaxation. Therefore, although the system itself 
has no intrinsic characteristic relaxation timescales for 
relaxation, it remains able to relax during an aging ex- 
periment in a finite time which is set by the experimental 
protocol. More quantitatively, it is found in the simplest 
cases that two-time correlators C(t,t w ), or equivalently 
response functions, can be decomposed in the following 



3 This is sometimes called relaxation time. It does not mean that 
the system relaxes after this time but just that its structure 
changes substantially, i.e. dccorrclatcs, over this time scale. 
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way: 



C(t,t w ) ~ C eq (t t w ) + C ag i n g(i/i w ), 



(76) 



in which case the relaxation time T a (t w ) is directly 
proportional to t w . This 'time-aging time' superposi- 
tion principle is reminiscent of the more standard 'time- 
temperature' superposition principle often found at ther- 
mal equilibrium. It is found in many (but not all) poly- 
meric liquids that r a grows sub-linearly with i w , 



T~a (t-w ) 



(77) 



where fi is the so-called aging exponent ( Struik. 19781) . 
which usually takes value in the range /i £ 1] (p 
is equal to 0.8 in Fig. [T5]) . The simple descrip- 
tion (|76[) has been challenged by dielectric experiments 
of aging of molecular liquids dLehenv and Nagel. 19981 : 
iLunkenheimer et ai, 20051: iRichert. 20101) . These studies 
have revealed that the aging part of the correlation func- 
tion cannot be rescaled just by changing the relaxation 
time, i.e. the time-aging superposition is not always sat- 
isfied. One has to stress however that in experiments one 
can access only the behavior at frequency much higher 
than the inverse of the decorrclation time. This is due to 
the fact that the quench can be only performed at finite 
speed, of the order of lK/min, in order to avoid temper- 
ature gradients in the sample. As a consequence, when 
the final temperature is reached, the decorrelation time 
has already grown so much that it is out of the exper- 
imental window. Thus, it is not really known whether 
the violation of the time-aging time superposition princi- 
ple affects only the high frequency tails of the dielectric 
susceptibility or whether Eq. (|76p is incorrect even for 
t of the order of the decorrelation time for many glassy 
liquids. 

Since the complete history of a sample in the glass 
phase matters, there is no reason to restrain experimen- 
tal protocols to the simple aging experiment mentioned 
above. Indeed, experimentalists have investigated 
scores of more elaborated protocols that have revealed 
an incre dibly rich, a nd sometimes quite unexpected, 



physics ( Young. 19981 ). a very old tr a dition in th e 
field of polymer glasses ( Kovacs. 1963t IStruik. 19781) . 
In particular, striking 'memory' and 'rejuvenation' 
effects a re observed during te mperature cycling exper- 
iments ( Refregier et ai, 19871) (one can also imagine 
applying an electric field or a mechanical constraint, be 
they constant in time or sinusoidal, etc.). These two ef- 
fects were first observed in spin glasses, but the protocol 
was th en repeated in many d ifferent materials, from poly- 
mers (IBellon et al, 20001: iFukao and Sakamoto. 20051: 
iMontes et ai, 2006 ) and organic liquids 



Yardimci and Lehenv. 20031 ) 
Doussineau et al, 19991) . 



(jLehenv and Nagel. 19981: 
to disordered ferroelectrics 

On top of being elegant and quite intriguing, such pro- 
tocols are relevant because they probe more deeply the 
dynamics of aging materials, allowing one to ask more 
precise questions beyond the simplistic observation that 



'this material displays aging'. Moreover, the observa- 
tion of similar effects in many different glassy materi- 
als implies that these effects are probably quite generic 
to systems with slow dynamics. Interesting also are the 
subtle differences possibly observed from one material to 
another. A large body of experimental, numerical and 
theoretical pap ers have been devote d to this type of ex- 
periments, see ( Berthier et al., 20031 ) for a more extensive 
review on this topic. 

A popular interpretation of the aging phenomenon was 
obtained by considering trap models (jBouchaud. 1992t 
iMonthus and Bouchaud. 1996t h In this picture, remi- 
niscent of the Go ldstein view of the glass transition 
mentioned above ( Goldstein. 19691 ) . the system is de- 
scribed as a single particle evolving in a complex en- 
ergy landscape with a broad distribution of trap depths. 
Thus, this is a paradigmatic mean-field approach. In 
particular, it is not easy to be very specific about 
what the 'traps' of the 'rugged landscape' represent 
in physical terms, although much effort was recently 
dedicated to establish c onnections between landscape 
and real space pictures dBerthier and Garrahan. 2003at 
iBertin. 20051: iHeuer. 20081) . 

The simplest version of the trap model makes no refer- 
ence to a spatial organization of the traps, and therefore 
does n ot take into account possible in teractions between 
traps ( Monthus and Bouchaud. 19 96;). The dynamics of 
the model in written in terms of an evolution equation 
for P(E,t), the probability that the system is in a trap 
of depth E at time £, and it is assumed that dynamics is 
thermally activated: 

= - To e-f> E P(E, t) + T(t)p(E), (78) 

where T(t) = T J dE e- 0E P(E,t) is the average hop- 
ping rate at time t, Tq being an attempt frequency. The 
complexity of the glassy material is now hidden into 
the probability distribution of the trap depth, p(E). In 
Bouchaud's trap model, an exponential distribution of 
trap depth is assumed, 



P{E) 



exp 



E 



(79) 



but the Gaussian trap model was also often con- 
sidered in the context of the equilibr i um dynam- 
ics of supercooled liquids dBassler. 1987 : iDvre. 19951 
IMonthus and Bouchaud. 19961: IHeuer. 20081) . In the per- 
spective of aging studies, the exponential trap distribu- 
tion is more relevant because it induces a broader distri- 
bution of trapping times, r, 



(T T) 1 + r / r o ■ 



(80) 



Remarkably, the first moment of this distribution di- 
verges at low temperatures, T < T , so that the sys- 
tem undergoes a true glass transition with a relaxation 
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timcscalc which diverges at To- By contrast, the Gaus- 
sian distribution does not yield to ergodicity breaking 
(i.e. aging is interrupted at long times) and it is thus 
mostly used in equilibrium studies, where it seems to 
account well for properties of the potent ial energy lan d- 
scape accessed in computer simulations ( Heuer. 20081 ). 



When the system is suddenly quenched below To, ag- 
ing behaviour arises because the system visits traps that 
are increasingly deep when t w increases, corresponding 
to more and more stable states. It takes therefore more 
and more time for the system to escape, and the dy- 
namics slows down with time, in a manner reminiscent 
of Fig. [TS1 More quantitatively, one can compute two- 
time correlation functions, such as the persistence func- 
tion C(t,t w ) = J dEP(E,t w ) exp[-(r e-P E )} which rep- 
resents the probability the system has not changed trap 
between times t w and t + t w : 

C (t,t w ) = /' (l-xf^x-^dx 



(81) 



C(t/t w ), 



which represents therefore an explicit example where the 
aging exponent defined in Eq. ()77|) can be computed ex- 
actly, ji = 1, and correlation functions are scaling func- 
tions of the rescaled time t/t w , as in Eq. (J7BJ). 

In Eq. (|81[) . the correlation function is computed as an 
ensemble average. Important physical insights are also 
provided by trap models beyond averages at the level 
of the fluctuations, since the existence of waiting time 
distributions with fat tails indeed implies that time se- 
ries are typically highly intermittent, and that run-to- 
run fluctuations are also large, suggesting that it is in- 
teresting to consider also statistical distributions of the 
fluctuations. For instance, one can easily imagine that 
van-Hove functions in aging s ystems described by tra p 
models are far from Gaussian ( Barkai and Cheng. 20031) . 
and that non-ergodic effects can be quite strong in sys- 
tems describ ed by Eq. (1501). see (iHe et al, 2008|) for an 
example and ( Margolin and Barkai. 200*61) for a more for- 
mal approach. 



B. Mean-field aging and effective temperatures 

Theoretical studies of mean-field glassy mod- 
els have provided important insights into 
the aging dynamic s of both structural and 

spin glasses (ICugliandolo and KurchanT 19931: 

ICugliandolo and KurchanT 19941) . We have already 
encountered these models in Sec. IIV.B.1I They provide 
a simple setting to study glassy models with a rugged 
free energy landscape. As a consequence it is natural to 
analyze their aging dynamics and use the corresponding 
results as a mean- field guide line for real systems. 

In Sec. IIV.B.11 we described several alternative theo- 
retical paths leading to essentially similar results for the 
equilibrium properties of glasses, which could all be de- 



scribed as 'mean-field'. Because aging studies now di- 
rectly deal with states that are non-stationary, proto- 
col dependent, and far from equilibrium, not all of the 
mean-field equilibrium approaches are easily extended to 
low temperature aging studies. We expect for instance 
that mode-coupling theory should be able to treat such 
time-dependent phe nomena, bu t this seems to be techni- 
cally quite involved ( Latz. 20"ool) and a compete extension 
of MCT to the aging regime remains an open problem. 
For completeness, we mention a recent work where the 
phenomenological RFOT mosaic scaling arguments were 
used to describe also nonequilibrium relaxation b elow the 
glass transition ( Lubchenko and Wolvnes. 20041 ). This 
framework makes contact with the older t heorical ap- 
proach of Nayaranaswamy-Moynihan-Tool (Tool. 19461 : 
iNaravanaswamv. 19711 : iMovnihan et al, 1976 ) but it also 
predicts important deviations from their findings. 

Even by focusing on fully connected disordered mod- 
els, often described as 'simple' models, we note that it 
took several years to derive a proper asymptotic solution 
of the l ong-time dynamics for a series of mean-field spin 
glasses ( Cugliandolo. 20031) . These results have then trig - 
gered an enormous activity dCrisanti and Ritort. 20031 ) 
encompassing theoretical, numerical and also experimen- 
tal work trying to understand further these results, and 
to check in more realistic systems whether they have some 
reasonable range of applicability beyond mean-field. This 
large activity, by itself, easily demonstrates the broad in- 
terest of these results. 

Mathematically, understanding the aging dynamics of 
mean-field glass models means solving a closed set of dy- 
namical equations. For concreten ess, let us consider th e 
following spin glass Hamiltonian ( Nieuwenhuizen. 1995|) , 



* = -E E 



j 



(82) 



P=2 ji<---<j f 



where Si(i = ,N) arc spin variables interacting 

through coupling constants which are random Gaussian 
variables of zero mean and variance p\ Jp / (2N P ^ 1 ). This 
model is a straightforward generalization of the p-spin 
model of Eq. (|40|) . We consider soft spin variables, 
which are real variables satisfying the spherical con- 
straint, Y^ii s f = N. Due to the mean-field nature of the 
Hamiltonian ([521 , a closed set of dynamical equations in- 
volving two -time correlation a nd response functions can 
be derived ( Cugliandolo. 2003T ). Defining 



I N 



N 



R(t, t v 



1 d( 8i (t)) 
N ^ dh t {t w ) 

I— 1 x ' 



(83) 



hi=0 



where hi is a magnetic field that couples to spin Si, and 

_j oo 

.9(z) = ? E J p V > (84) 



o=2 



40 



one gets the time evolution of the two-time dynamic func- 
tions following a quench from a completely disordered 



J 



state at t v 



dC{t gt tw) = - n(t)C(t, f w ) + 2TR(t, t w ) + J J dt'D(C(t, t'))R(t w ,t') 

dR{t gt tw) = - n(t)R(t, t w ) + 6(t - t w ) + dt'£(t, t')R(t', 

n(t)=T+ f dt'\D{C{t,t'))R{t,t') + T,{t,t')C{t,t')], 
Jo 



dt'E(t,t")C(t",t w ), 



(85) 



r 



where the kernels are defined as 

D(t,t w ) =g'[C(t,U)}, S(t,t w ) = g"[C(t,t w )]R(t,U). 

(86) 

The unique feature that makes the dynamics of mean- 
field spin glass models soluble is that the dynamical equa- 
tions (|85p are closed and only involve two-point functions. 
This great simplification stems from the fully-connected 
nature of the Hamiltonian ([82} , and allows one to for- 
mulate an exact as ymptotic solution for the dynami cs of 
mean- field models (|Cugliandolo and Kurchan. 19931 ). 

A comparison with Eq. (|42[) immediately reveals why 
the aging regime is much harder to treat analytically, 
since one has to face two difficulties. First, two-time cor- 
relation functions now depend on both their arguments, 
C(t,t w ) C(t — t w ). Second, the equations of motion 
(|85p in the aging regime not only involve time correla- 
tions, but also time-dependent response functions. At 
thermal equilibrium response and correlations are not 
independent, since the fluctuation-dissipation theorem 
(FDT) relates both quantities, 



be sh own that two-time correlation functions take th e 
form (Cuglian dolo. 20031 lAndreanov and LefevreT 20061) 



R(t,t w ) — 



1 dC(t,U 
T dt w 



(87) 



Indeed, imposing both FDT and time translational in- 
variance in Eqs. ([55)1 yields the much simpler Eq. (|4"9"} 
for the p-spin model. 

The solutio n to Eqs. (1551 has been reviewed in great 
detail before ( Cugliandolo. 20031 ). so we only briefly de- 
scribe the most striking physical outcomes. As men- 
tioned in Sec. IIV.B.11 in these mean-field models, ther- 
mal equilibrium is never reached when the quench is per- 
formed below a critical temperature T c , below which the 
relaxation time is infinite. It can be shown that aging 
proceeds forever by downh ill motion in an increasingl y 
fiat free energy landscape ( Kurchan and Laloux. 199a ). 
with subtle differences between spin glass and structural 
glass models. In both cases, however, time translational 
invariance is broken, and two-time correlation and re- 
sponse functions explicitly depend on both their time ar- 
guments. When g(x) contains a single term with p > 2, 
as for instance in the spherical p-spin model, it can 



C(t,t w ) ~ C C q(t t w )+C 



( Ht) 

U(*w) 



h{t) = exp(t 1 



where n(T) is an exponent that cannot be computed an- 
alytically. Numerical solutions show that < /i(T) < 1 
and n{T — > 0) — >• 1. The scaling form ([88} is appealing 
since it is very similar to the empirical form in Eq. (|76l) . 
and provides an explicit example where the use of the 
aging exp onent ii defin ed in Eq. ([77} and introduced 
by Struik ( Struik. 19781 ) is analytically justified. From a 
broader perspective, it is interesting to note that the ex- 
act dynami c solution of the equa tions of motion for time 
correlators ( Kim and Latz. 200l[ ) displays a behaviour in 
strikingly good agreement with the numerical results re- 
ported for structural glasses, e.g., in Fig. [TBI 

Note that a more complex time dependence for C(t, t') 
can be obtained for some mean-field glassy systems. 
These generally display a phenomenology reminiscent of 
spin-glasses and not structural glasses (the transition is 
continuous, the spin glass susceptibility diverges, etc). 
This is for example the case of the 'p = 2 + 4 model' for 
which g(x) takes a more complex form. In these cases the 
scaling form with a single function h(t) in Eq. ([88} docs 
not hold anymore, but has to be generalized to include a 
continuous hierarchy of such functions. Physically, this 
implies that the relaxation of correlation functions in the 
aging regime is associated to an infinite number of aging 
timcscalcs, which all diverge with the age of the system. 

The temporal behaviour of time correlation functions 
already shows that mean-field spin glass models display 
a rich aging phenomenology In aging systems, however, 
there is no reason to expect the FDT in Eq. ([87} to hold 
and both correlation and response functions carry, at 
least in principle, distinct physical information. Again, 
the asymptotic solution obtained for mean-field models 
quantitatively establishes that the FDT does not ap- 
ply in the aging regime. The solution also shows that 
a gene ralized form of the FDT holds at large waiting 
times ( Cugliandolo and Kurchan. 19931) . This general- 
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izcd form of the FDT reads 



R(t, t w ) 



X(t,t w ) dC(t,t w ) 



T 



(89) 



which defines X(t,t w ), the so-called fluctuation- 
dissipation ratio (FDR). At equilibrium, correlation and 
response functions are time translational invariant and 
equilibrium FDT imposes that X(t,t w ) = 1 at all times. 
A parametric fluctuation-dissipation (FD) plot of the 
step response or susceptibility 



against 



X (t,t w ) = T J dt'R(t,t r ) 



AC(t,t w ) = C(t,t)-C(t,t w ), 



(90) 



(91) 



is then a straight line with unit slope. These 
simplifications do not necessarily occur in a non- 
equilibrium system. But the definition of an 
FDR through Eq. (|89|) becomes significant for mean- 
field aging systems ([Cugliandolo and Kurcham 19931: 
ICugliandolo and Kurchan. 19941) . In mean-field spin 
glass models the dependence of the FDR on both time 
arguments is only through the correlation function, 



(92) 



valid at large wait times, t w — > oo. 

For mean-field structural glass models (such as the p- 
spin model with p > 2), time correlation functions dis- 
play a two-step relaxation process as in Eq. . Corre- 
spondingly, the simplification (|92[) is even more spectac- 
ular since the FDR is shown to be characterized by only 
two numbers instead of a function, namely X ~ 1 at short 
time differences (large values of the correlator) corre- 
sponding to a quasi-equilibrium regime, with a crossover 
to a non-trivial number, X ~ X°° for large times (small 
values of the correlator). This implies that parametric 
FD plots are simply made of two straight lines with slope 
1 and X°°, instead of the single straight line of slope 1 ob- 
tained at equilibrium. Formally, the infinite-time FDR, 
X°°, is defined as 



A : 



lim lim X(t,t w ). 



(93) 



These theoretical predictions were tested with suc- 
cess in numerical simulations ( Barrat and Kob. 19991 



iDi Leonardo et ai. 20001) . In Fig. UM we present more 
recent numerical dat a obtained in an aging silica 
glass ( Berthier. 2007bl ). presented in the form of a para- 
metric response-correlation plot. The measured correla- 
tion functions are the self-part of the intermediate scat- 
tering functions defined in Eq. ©, while the conjugated 
response functions quantify the response of particle dis- 
placements to a spatially modulated field conjugated to 
the density. Plots for silicon and oxygen atoms at dif- 
ferent ages of the system are presented. They seem 



x 




FIG. 19 Parametric correlation-response plots measured in 
the a ging regime of a numerical model for a silica glass, 
Si0 2 (Berthier. 2007bt ). The plots for both species smoothly 
converges towards a two-straight line plot of slope 1 at short 
times (large C values), and of slope X°° « 0.51 at large times 
(small values of C) , yielding an effective temperature of about 
T cS = T/X°° « 4900 K. Note that T eS > T M ot « 3330K > 



to smoothly converge towards a two-straight line plot, 
as obtained in mean-field models for structural glasses 
(note, however, that this could be a pre-asymptotic, fi- 
nite "£ w " , effect). Note also that the value of T e g is much 
larger than other relevant temperature scales, in partic- 
ular than Tmcti while both quantities nearly coincide in 
mean-field models. 

The behaviour of spin glass mean-field models, such 
as the p = 2 + 4 model, is again more involved. It is 
found that to each timcscale of the continuous hierarchy 
corresponds to a given value of the FDR. This means that 
the FDR has become a continuous function, X ~ x(C) 
in the aging regime, instead of the single number found 
for the p-spin model. For the spherical p = 2 + 4 model, 
for instance, one finds 



x(C) 



T g'"(C) 
2 [g"{C)Y/ 2 ' 



(94) 



for C values corresponding to the aging regime; g{x) was 
defined in Eq. (|84|) . Thus, the parametric FD plot is now 
made of a straight line of slope 1 for large value of the 
correlation function (quasi-equilibrium regime), followed 
by a continuous curve with slope < 1 for smaller values 
of the correlation, that differs from the FDT even in the 
infinite time limit, i w — > oo. 

Since any kind of behaviour is in principle allowed 
in non-equilibrium situations, getting such a simple, 
equilibrium-like structure for the FD relations is a re- 
markable result. This immediately led to the idea that 
aging systems might be characterized by an effective ther- 
modynamic behaviour and th e idea of 'effective equilibra - 
tion' at different timescales ( Cugliandolo et ai, 1997al ). 
In particular, generalized FD relations in Eq. (|g9"]l sug- 
gest to define an effective temperature, as 



T 



A(Mw) : 



(95) 
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such that mean-field structural glasses are char- 
acterized by a unique effective temperature, 
T ff = T/X°°, in their aging regime. It can be in- 
terpreted as the temperatur e at which slow modes 
are quasi-equilibrated ( Cugliandolo et al.. 1997at 
iFranz and Virasoro. 2000l ). One finds in general that 
< X°° < 1, such that T cS > T, as if the system had 
kept some memory of its high temperature initial state. 
It was then proposed that FDRs, or cquivalcntly effective 
temperatures, can be m easured by focusing on an y type 
of physical observables ( Cugliandolo et al.. 1997al ) . The 
FDR is then defined in terms of the two-time connected 
correlation function for generic physical observables A(t) 
and B(t), 

C AB (t,t w ) = (A(t)B(t w )) - (A(t))(B(t w )), (96) 

with t > t w , and the corresponding two-time (impulse) 
response function 



RAB(t,t w ) = 



S(A(t)) 



Sh B (t w ) 



(97) 



Here h denotes the thcrmodynamically conjugate field to 
the observable A so that the perturbation to the Hamil- 
tonian is SH = —hs-B. A practical measurement of the 
FDT is then performed by building a parametric FD plot 
of the integrated response function, 



XAB = T dt'RAB{t,t'), 



(98) 



versus the correlation function [Cab(£, i) — CAB(t,t w )}. 

The name 'temperature' for the quantity defined in 
Eq. (|95p is not simply the result of a dimensional anal- 
ysis but has a deeper, physically appealing meaning 
that is revealed by asking the following questions. How 
does one measure temperatures in a many-body sys- 
tem whose relaxation involves well-separated timescales? 
What is a thermometer (and a temperature) in a far 
fro m equilibrium aging mater i al? Answers ar e provided 
in (jCugliandolo et al. 1997at iKurchan. 20051 ) both for 
mean-field models and for additional toy models with 
multiple relaxation timescales. The idea is to couple an 
additional degree of freedom, such as a harmonic oscilla- 
tor, x(t), which plays the role of the thermometer oper- 
ating at frequency u>, to an observable of interest A(t) via 
a linear coupling, — \x(t)A(t). Simple calculations show 
then that the thermometer 'reads' the following temper- 
ature, 



-k R T„ 



-uj 2 (x 2 ) 



2R AA (u;,t w ) 



(99) 



where C' aa (oj, t w ) is the real part of the Fourier trans- 
form of Eq. (|9l))) , and R AA (uj,t w ) the imaginary part of 
the Fourier transform of Eq. (|97|) . with h = Xx. The 
relation (|99[) indicates that the bath temperature is mea- 
sured, T moas = T, if frequency is high and FDT is satis- 



to be tuned to that of the slow relaxation in the aging 
material. More generally, relaxation in mean-field glassy 
syste ms may occur in several, well-s eparated time sec- 
tors ( Cugliandolo and Kurchan. 19941 ) . It is then easy to 



imagine that each sect or could be asso ciated with its own 
effective temperature ( Kurchan. 20051 ). 



The direct link between the FDR in Eq. (JHSJ) and the 
effective temperature measured in Eq. (|99[) was numeri- 
cally confirmed in computer simulations of gla ssy molec- 
ular liquids. In ( Berthier and Barrat. 2002a]) . a tracer 



particle was used as a thermometer, its frequency be- 
ing tuned by modifying its mass. It was verified that 
its kinetic energy was controlled by the temperature of 
the heat bath for light (high frequency) tracers, while it 
was related to the effective temperature of the slow de- 
grees of freedom for heavy (low frequency) tracers. In the 
same spirit, a t wo-level system with a djustable frequency 
was studied in (|llg and Barrat. 20071) . and the activation 



rate changed from being proportional to exp(—E/T) to 
cxp(— E/T e g) when decreasing the frequency. These two 
examples show that T c ff truly deserves the name of a 
temperature in a fundamental sense. We note that these 
modern concepts are related to, but make much more pre- 



~- 

tures in aging glasses ( 


Tool. 19461: iNaravanaswamv. 197ll 


iMovnihan et al.. 1976 


) . Attempts were even made to de- 



velop a thermodynamic of the glass state, heavily relying 
on the idea of effecti ve temperatures. This has becom e 
the subject of a book ( Leuzzi and Nieuwenhuizen. 20071) . 



The final piece of information extracted from the 
behaviour of mean-field spin glass models is the exis- 
tence of a connection between a non-equilibrium dynamic 
quantity, namely the asymptotic FDR x(C) defined in 
Eq. (|92j) . and the thermodynamic behaviour of the sys- 
tem in the low temperature phase. It turns out that 
the thermodynamics of mean-field models for structural 
glasses is characterized by one-step replica symmetry 
breaking, while a full-step solution is need ed to solve th e 
thermodynamics of models for spin glasses ( Parisi. 2003[) . 



Remarkably, the structure of the spin glass order param- 
eter, the Parisi function P cq (q) describing the probability 
distribution of overlaps between equilibrium states, is di- 
rectly related to the structure of the function x(C). For 
the p = 2+4 model, for instance, s tatic calculations yield s 



an explicit expression for P cq ( Nieuwenhuizen. 19951) 



Comparing with Eq. (fM)) . it turns out that the following 
equality holds: 



x(C) 



dqP cq {q). 



(100) 



ficd, while Tj T 



T e ff > T if frequency is slow enough 



The situation for mean-field structural glasses is more 
complicated, since Eq. (jlOOp does not hold, but the in- 
tegrated Parisi function (the right hand side in (|100|) ) 
has t he same structure as the FDR (the left hand 
side) (jCugliandolo and Kurchan. 19931 ). Therefore, the 
full-step or one-step replica symmetry breaking schemes 
needed to solve the static problem in these models 
have a direct dynamical counterpart, the FDR be- 
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ing a function or a number, respectively, in the ag- 
ing regime. It was further argued that Eq. (|100[) 
might hold for finite d imensional glassy systems as 
well (jFranz et al. 1998bl) . raising the exciting possibil- 
ity that the Parisi function might become experimentally 
accessible through agi ng experiments, which trig gered a 
large research activity ( Crisanti and Ritort. 20031) . How- 
ever, it is not easy to deter mine the cond itions under 
which Eq. ([TOO)) might hold (jParisi. 20031 ). Moreover, 
aging studies are very often performed very far from any 
asymptotic regime, and very little is know about how 
Eq. (|10JJ|) is modified when the £ w — » oo limit cannot b e 
taken ( Barrat and Berthier. 2001al ICastillo et al. 20021) . 
making the thermodynamic interpretation of the out- 
come of aging measurements delicate. 

Taken together, these results make the mean-field de- 
scription of aging very appealing, and they nicely com- 
plement the mode-coupling/RFOT description of the 
equilibrium glass transition described above. More- 
over, they have set the agenda for a very large 
bo dy of numerical and exper imental work, as reviewed 
in ( Crisanti and Ritort. 2003T ). It should be clear, how- 
ever, that these results are strictly valid in the mean-field 
limit in the sense discussed in Sec. IIV.B.1I Nonequi- 
librium aging dynamics in mean-field spin glasses turns 
out to describe the slow descent of the system in an 
energy landscape wh ich becomes more and mo re flat 
as the age increases ( Kurchan and Laloux. 19961) . with 
no access to the deeper metastable states that are sup- 
posed to play an important role in glasses near the ex- 
perimental glass transition. This view is thus in strik- 
ing contrast with the purely activated description given 
for instance by trap models. Additionally, it is impor- 
tant to understand the role of spatial fluctuations which 
are not naturally included in the mean-field d escription, 
see (|Castirlo et al. 2002t [Ch anion et al. 20041) for recent 
work in this direction. Thus it is important to test further 
the mean-field concepts to understand how they apply to 
the three-dimensional world. 



C. Beyond mean-field: Experiments, critical points and 
kinetically constrained models 

Despite successes such as shown in Fig. [121 
the broader applicability of the mean-field sce- 
nario of aging dynamics remains unclear, however. 
While some experiments and simulations indeed 
seem to support the existence of well-behaved ef- 
fective 



found non-m onotonic or ev en ne gative response 
functions (IViot et al. 20031: iNicodemi. 1999t 



temperature s 



Abou and Gallet. 2 004: 



(jGrigera and Israelo ff. 1999: 
~ IWang et al. 2006 ; 



Oukris and Israeloff. 2010|) . other studies also reveal the 
limits of the mean-field scenario. Experiments have for 
instance reported anomalously large FDT violations as- 
sociated with intermittent dynami cs dBellon et al. 2001 ; 
Bellon and Ciliberto. 2001 iBuisson et al. 2003a: 



Buisson et al. 2003bl : iGreinert et al. 20061 ) . while 
theoretical studies of model systems have also 



iKrzakala. 20051 : Depken and Stinchcombe. 2005[) . 

and ill- defined or observable-dependent 

FDRs ([Fielding and Sollich. 20"02l) . Some ex- 

periments even reported an absence of vi- 
olation s to the fluctuation -di ssipation theo - 
rem (jJabbari-Farouii et al. 20081 : IJop et al. 20091) . 
In principle, these discrepancies with mean-field predic- 
tions are to be expected, since there are many systems 
of physical interest in which the dynamics are not 
of mean-field type, in particular displaying activated 
processes. However, it is not possible to draw any 
consistent picture from the experiments at this stage. As 
a consequence it is not quite clear yet how and to what 
extent mean field results are violated in real systems. 

It is thus an important task to understand from the 
theoretical point of view when the mean-field concept 
of an FDR-related effective temperature remains viable. 
However, studying theoretically the interplay between 
relevant dynamic lengthscales and thermally activated 
dynamics in the non-equilibrium regime of disordered 
materials is clearly a challenging task. As a consequence, 
much work has been devoted to analyze simple effective 
models. A lot of attention has been focused on spin mod- 
els. 

A first class of system that we discuss are coarsen- 
ing systems. Although not directly related to the glass 
problem, they provide a simple, yet non-trivial, theoreti- 
cal framework to study situations where both aging and 
spatial heterogeneity are present, and where time corre- 
lation and response functions display interesting scaling 
behaviour. The paradigmatic situation is that of an Ising 
ferromagnetic model (with a transition at T c ) suddenly 
quenched in the ferromagnetic phase at time t w = 0. For 
i w > domains of positive and negative magnetizations 
appear and slowly coarsen with time. The appearance of 
domains that grow with time proves the presence of both 
aging and heterogeneity in this situation. 

The case where the quench is performed down to 
T < T c is well understood. The only relevant lengthscale 
is the growing domain size, £(t w ) and the physical be - 
havior can be understood by scaling theory ( Brav. 19941) . 
Correlation functions display aging, and scale invariance 
implies that C(t,t w ) ~ f (£(t) / £(t w )) . Resp onse functions 
can be decomposed i nto two contributions ( Barrat. 19981 : 
iBerthier et al. 19991 ): one part stems from the bulk of 
the domains and behaves as the equilibrium response, 
and a second one from the domain walls and becomes 
vanishingly small in the long time limit where £(t w ) — > oo 
and the density of domain walls vanishes. This im- 
plies that for coarsening systems in d > 2, one has 
X°° = 0, or equivalently an infinite effective temper- 
ature, T c ft = oo. The case d = 1 is special be- 
cause T c — and the response function remains dom- 
inated by the domain wal ls, which yields the non- 
trivial value X°° = 1/2 dGodreche and Luck. 2000at 



iLippiello and Zannetti. 2000f ) 
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Another special case has retained attention. When the 
quench is performed at T = T c , there is no more dis- 
tinction between walls and domains and the above argu- 
ment yielding X°° = does not hold. Instead one stud- 
ies the growth with time of critical fluctuations, with 

1 1 z 

£(iw) ~ t-w ° the correlation length at time t w , where 
z c is the dynamic critical exponent. Both correlation 
and re sponse functions become non -trivial at the critical 
point ( Godreche and Luck. 2000bh . It proves useful in 
that case to consider the dynamics of the Fourier com- 
ponents of the magnetization fluctuations, 



C q (t,t w ) = (m q (t)m_ q (t w )), 
and the conjugated response 



5{m q {t)) 



(101) 



(102) 



h_ o =0 



where h q (t) is the Fourier component of the mag- 
netic field at time t. From Eq. (|59"|) a wavevector 
dependent FDR follows, X q (t,t w ), which has interest- 
ing properties that can be computed by a number 
of means, i ncluding dynamical renormaliz ation tech- 
niques, see ([Calabrese and Gambassi. 2005|) for a re- 
view. One of the main outcomes of these studies is 
that the effective temperature for quenches at the crit- 
ical point might in some cases dep end on the observ- 
able ( Calabrese and Gambassi, 20041) and on the initial 
condition ( Calabrese et al, 20061) . thus weakening the in- 
terpretation of X q in terms of effective temperature. 

In dimension d = 1, it is possible to compute X q (t, t w ) 
exactly in the aging regime at T = T c = 0. An interest- 
ing scaling form is found, and numer ical simulations per - 
formed for d > 1 confirm its validity (jMaver et al., 20031 ) : 



X q (t,t w )^ *[gf(t w )], 



(103) 



where the scaling function X(x) is X{x — > oo) — > 1 at 
small lengthscale, q£ > 1, and X(x -> 0) -> AT 00 = 1/2 
(in d = 1) at large distance, q£ <C 1; recall that z c = 2 in 
that case. 

Contrary to mean-field systems where geometry played 
no role, here the presence of a growing correlation length- 
scale plays a crucial role in the off-equilibrium regime 
since £(i w ) allows one to discriminate between fluctua- 
tions that satisfy the FDT at small lengthscale, X q ~ 1, 
and those at large lengthscale which are still far from 
equilibrium, < X q ~ A°° < 1. These studies sug- 
gest therefore that generalized fluctuation-dissipation re- 
lations in fact have a strong lengthscale dependence — a 
result which is not predicted using mean-field approaches. 

Another interesting result is that the FDT violation 
for global observables (i.e. those at q = 0) takes a par- 
ticularly simple form, since the introduction of a sin- 
gle number is sufficient, the FDR at zero wavevector, 



X q= o(t, t w ) 



X'- 



1/2 (in d = 1). This univer- 



sal q uantity takes non-trivial va lues in higher dimen- 
sion (jGodreche and Luck. 2000bl ). e.g. X°° w 0.34 is 



J? -2 - 
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l-C q (t,t w )/C q (t,t) 



FIG. 20 Parametric response-correlation plots for the Fourier 
components of the mobility field in the d = 3 Fredrickson- 
Andersen model. Symbols are from simulations, lines from 
analytic calculations, and wavevectors decrease from top to 
bottom. The FDT is close to being satisfied at large q corre- 
sponding to local equilibrium. At larger distance deviations 
from the FDT are seen, with an asymptotic FDR which be- 
comes negative. Finally, for energy fluctuations at q = 
(bottom curve), the plot becomes a pure straight line of (neg- 
ative!) slope —3, as a result of thermally activated dynamics. 



measured in d = 2 ( Maver et al., 20031 ). This shows that 
the study of global rather than local quantities makes 
the measurement of X 00 much easier. Finally, having a 
non-trivial value of X°° for global observables suggests 
that the possibility to define an effective temperature re- 
mains valid, but it has become a more complicated ob- 
ject, related to global fluctuations on large lengthscale. 
The first experimental determination of the value of X°° 
near a critical point was reported only very recently in a 
system of liquid crystals, wh ere the director is th e rele- 
vant fluctuating observable ( Joubaud et al., 2009t ). An 
intriguing value X°° ~ 0.31 was reported, which has 
received, to our knowledge, no theoretical justification 



( Calabrese and Gambassi. 20051) . 

Kinetically constrained spin models represent a second 
class of non-mean-field systems wh ose off-equilibrium ha s 
been thoroughly studied recently ( Leonard et al., 2007j) . 
This is quite a natural thing to do since these systems 
have local, finite ranged interactions, and they combine 
the interesting features of being defined in terms of (effec- 
tive) microscopic degrees of freedom, having local dynam- 
ical rules, and displaying thermally activated and hetero- 
geneous dynamics. 

The case of the Fredrickson- Andersen model, de- 
scribed in Sec. IIV1 ha s been studied in great de- 



tail ( Leonard et al., 20071 ). and we summarize the main 
results. Here, the relevant dynamic variables are the 
Fourier components of the mobility field, which also cor- 
respond in that case to the fluctuations of the energy 
density. Surprisingly, the structure of the generalized 
fluctuation-dissipation relation remains once more very 
simple. In particular, in dimension d > 2, one finds a 
scaling form similar to Eq. (|f 03|) . X q (t,t w ) = X[q^(t w )], 
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with a well-defined limit at large distances, X g —o(t, t w ) = 
X°° . The deep analogy with critical Ising models stems 
from the fact that mobility defects in KCMs diffuse in a 
way similar to domain walls in coarsening Ising models. 
It is in fact by exploiting this analogy that analytic re- 
sults are obtaine d in the aging regime of th e Fredrickson- 
Andcrscn model ( Maver and Sollich. 20071 ). 

There is however a major qualitative difference be- 
tween the two families of model. The (big!) surprise 
lies in the s ign of the asymptot ic FDR, since calculations 
show that ([Mayer et al, 20061) 



-3, d>2. 



(104) 



q=o(t,t w ) — f{t/t w ) 

w -3.307. Numeri- 



In dimension d = 1, one finds X 
with X q=Q (t oo,t w ) = 
cal simulations in various spatial dimensions nicely con- 
firm these calculations. In Fig. [20l we show such a 
comparison between simulations (symbols) and theory 
(lines) in the case of the d = 3 Fredrickson-Andersen 
model ([Maver et al., 20061) . Fourier components of the 
mobility field yield parametric FD plots that follow scal- 
ing with the variable q£(t w ), as a direct result of the 
presence of a growing lengthscale for dynamic hetero- 
geneity, with a simple diffusive behaviour in that case, 
£,(t w ) ~ v^w- Again, generalized fluctuation-dissipation 
relations explicitly depend on the spatial lengthscale con- 
sidered, unlike in mean-field studies. In Fig. [201 the limit 
q = corresponding to global obscrvables is also very 
interesting since the plot is a pure straight line, as in 
equilibrium. Unlike equilibrium, however, the slope is 
not 1 but —3. A negative slope in this plot means a neg- 
ative FDR, and therefore suggests a negative effective 
temperature, a very non- intuitive result at first sight. 

Negative response functions in fact directly follow from 
the ther mally activated nat ure of the dynamics of these 
models ( Maver et al., 20061 ). First, one should note that 
the global observable shown in Fig. [20] corresponds to 
fluctuations of the energy, e(t w ), whose conjugated field 
is temperature. In the aging regime the system slowly 
drifts towards equilibrium. Microscopic moves result 
from thermally activated processes, corresponding to the 
local crossing of energy barriers. An infinitesimal change 
in temperature, T —> T + ST with ST > 0, accelerates 
these barrier crossings and makes the relaxation dynam- 
ics faster. The energy response to a positive tempera- 
ture pulse is therefore negative, Se < 0, which directly 
yields Se/ST < 0, and explains the negative sign of the 
FDR. This line of reasoning does not apply to mean- 
field glasses, where thermal activation plays no role. To 
our knowledge, negative effective temperatures in realis- 
tic aging materials has not yet been observed. 

Finally, another scenario holds for local observables in 
some KCMs when kinetic constraints are stronger, such 
as the East model ( Leonard et al., 20071) or a bidimen - 
sional triangular plaquette model ( Jack et al., 2006bT) . 
Here, relaxation is governed by a hierarchy of energy 
barriers that endow the systems with specific dynamic 
properties. In the aging regime following a quench, in 



particular, the hierarchy yields an energy relaxation that 
arises in discrete steps which take place on very different 
timcscales, reminiscent of the 'time sectors' encountered 
in mean-field spin glasses (but not in mean-field struc- 
tural glass models!). Surprisingly, it is found that to 
each of these discrete relaxations one can associate a well- 
defined (positive) value of the fluctuation-dissipation ra- 
tio, again reminiscent of the dynamics of mean-field spin 
glass models. Therefore, even in models that are very 
far from the mean-field limit the physical picture of a 
slow relaxation taking place on multiple timcscales with 
each timescale characterized by an effective temperature 
seems to have some validity. 

In conclusion, we have described multiple non mean- 
field situations where collective behaviour (critical point 
in standard ferromagnets, dynamic criticality emerging 
in kinetically constrained models) produce non-trivial ag- 
ing dynamics characterized by the emergence of univer- 
sal fluctuation-dissipation properties. These results are 
instructive for understanding how and to what extent 
mean-field results can change in finite dimensional sys- 
tems. However, it is not clear that aging in structural 
glasses may be understood in terms of domain coarsen- 
ing (domains between what?), nor that kinetically con- 
strained models are faithful effective models of glasses. 
As a consequence, it is difficult to reach for theoretical 
models a consistent picture of effective temperature and 
aging valid for finite dimensional structural glasses. How- 
ever these results suggest that, most certainly, noncqui- 
librium effective temperatures have a validity much be- 
yond the realm of mean-field disordered spin models 
where they were first discovered ( Kurchan. 2005f ). 



D. Driven dynamics of glassy materials 

We have introduced aging phenomena with the argu- 
ment that in a glass phase, the timescale to equilibrate 
becomes so long that the system always remembers its 
complete history. This is true in general, but one can 
wonder whether it is possible to invent a protocol where 
the material history could be erased, and the system 
'rejuvenated' (jMcKenna and Kovacs. 19841) . This con- 
cept has been known for decades in the field of polymer 
glasses, where complex thermo-mechanical histories are 
often considered, for obvious practical reasons. 

Let us consider an aging protocol where a supercooled 
liquid is quenched to a low temperature at time t w = 0, 
but is simultaneously forced by an external mechanical 
constraint. Experimentally, one often finds that a sta- 
tionary state can be reached, which explicitly depends on 
the strength of the forcing: a system which is forced more 
strongly relaxes faster than a material which is less per- 
turbed, a phenomenon called 'shear-thinning'. The ma- 
terial has therefore entered a driven steady state, where 
memory of its age is no longer presen t and dynamics has 
become stationary: aging is stopped ( Kurchan. 200 1[ ). 

Many studies of these driven glassy states have 
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been performed in recent years. Note that the drive 
is in general mechanical and, hence, these are rel- 
evant for the rheology of supercooled liquids and 
glasses, and the T -C T g limit corresponds to stud- 



a broad field in itself, see ( 


. — „ „ , ... — - 

Falk and Langer. 1998 


Malonev and Lemaitre, 20061: 


iTanguv et al, 2006 


Schall et al, 20071 iBocauet et al, 20091) for a few mod- 



ern perspectives in this direction. In the colloidal world, 
such studies are also relevant for the newly-defined 
field of the rheology of 'soft glassy materials'. These 
materials are (somewhat tautologically) defined as 
those for which the non-linear rheological behaviour 
is believed to result precisely from the competition 
between intrinsically slow rela xation processes of glassy 
origin, and an external forcing ( Sollich et al., 1 997). It is 
believed that the rheology of dense colloidal suspensions, 
foams, emulsions, pastes, or even biophysical systems 
are ruled by such a competition , quite a broad field of 
application indeed ( Cates. 20031 ). 



From the point of view of statistical mechanics 
modeling, the rheology of soft glassy materials can 
be naturally studied from the very same angles as 
the equilibrium glass transition a nd aging phenom- 
ena itself. A s such, trap models ( Sollich et al., 19971 
ISollich. 1991 iFielding et al, 2000l). mean-field spin 



glasses ( Cugliandolo et al, 1997b : Berthier et al, 2000T ) 
and the related mode-coupling theo 



proaches ( 


Mivazaki and Rcichman. 2002; 


Mivazaki et al, 2004; 


Mivazaki et al, 2006 


Fuchs and Cates. 2002; 


Fuchs and Cates. 2005 


Brader et al, 20071 Fuchs and Cates. 20091) have been 



explicitly extended to include an external mechanical 
forcing. In all these cases, one generically finds that a 
driven steady state can be reached and aging is indeed 
expected to stop at a level that depends on the strength 
of the forcing. A most interesting aspect is that the 
broad relaxation spectra predicted to occur in glassy 
materials close to a glass transition directly translate 
into 'anomalous' laws both for the linear rheological 
behaviour (seen experimentally in the broad spectrum 
of elastic, G'(oj), and loss, G"(cj), moduli) 4 , and the 
non-linear rheological behaviour (such as a strong 
dependence of the steady state viscosity r\ upon an 
imposed shear rate 7), while aging and rheology would 
compete deep in the glass phase to produce interesting 
phenomena such as slow creep behaviour. However, 
all these different theoretical approaches have their 
strengths and weaknesses, as we now briefly describe. 

The trap model described in Sec. IV. Al was 
extended into the 'soft glassy rhe ology' (SCR) 



model (|Sollich et al, 19971 : ISollich. 199$ ). Its simplic- 



ity makes it a nice tool to investigate complex thermo- 



mcchanic al histories where agi ng and mechanical forcing 
compete ( Fielding et al, 2000h . leading to several inter- 
esting predicted behaviours. Thus, the model is often 
used by experimentalists to rationalize non-trivial rhe- 
ological behaviour s commonly encoun tered in complex 
materials, see e.g. ( Rouver et al, 20081 ). 

The SGR model is a direct extension of the trap model, 
where each trap is characterized not only by its depth, 
E, but additionally by a 'strain' variable, I. The evo- 
lution equation now involves P(E,£,t), which general- 
izes P(E, t) to include the statistical fluctuations of the 
strain variable. Using notations similar to the ones used 
in Eq. (|78p . the dynamics of the SGR model is defined 
as: 



T(t)p(E)S(£), 



(105) 



where T(t) = T Q J d£ j dE e -^ E -^ ke ^ P{E, £,t). The 
first term represents the effect of the global shear rate, 7, 
in the absence of hopping between traps, namely affinc 
deformation of the traps. The second term describes the 
probability to leave the trap occupied at time t, and takes 
into account in a linear manner the fact that shearing 
promotes hopping between traps by lowering the barrier 
heights, thereby defining an elastic constant, k. Note 
that activated dynamics with an 'effective temperature 
factor', T e ff = 1//3, is assumed, even though its mean- 
ing in the context of a driven colloidal system or emul- 
sion is not clear (jCates. 20031 ). The last term now in- 
cludes the new factor 6(£), implying that after a hop, the 
newly found trap is unstrained. Of course, for such a 
mean-field description (a single particle hopping between 
traps), 'strain', 'shear rates', or 'shear stress' are just 
names and are not intended to carry physical information 
about a real three-dimensional flow. As in the original 
trap model, an exponential for m is adopted for th e dis- 
tribution of trap depths, p{E) (jSollich et al, 19971 ). The 
global shear stress is defined by 

<j{t) = k(£) P = k d£ dEP(E,£,t)£, (106) 



4 G(ui) relates the stre ss to an imposed sin usoidal strain oscillating 
at frequency ui, see l lBarnes et al., 199ll ). 



so that knowledge of P{E, £, t) from solving Eq. (|105[) 
allo ws one to pr edict any needed rheological quan- 
tity ( Sollich. 19981 ). The success of the SGR model stems 
partly from the fact that depending on the value of 
the effective temperature, a broad variety of non-trivial, 
but experimentally realistic, rheological behaviour can 
be predicted both for linear and non-linear response. 
The model can therefore easily be used to fit a set of 
experimental data by adjusting a few quantities, such 
as T c ff or r , see ( Fabry et al, 200il ) for a biophysi- 
cal example. Moreover, it is simple enough that ex- 
tremely complex thermo-mechanical histories can be 
easily implemented, and com pared to experiments, see 
( Viasnoff and Leaueux. 20021) for an elegant illustration. 
The weaknesses of the approach are the same as for the 
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original trap model, as far as the interpretation of the 
traps in real space is concerned. Moreover, the lack 
of a spatial representation of the physics implies that 
the different traps are completely independent from one 
another, and the model cannot describe shear and ki- 
netic heterogeneities. It would be interesting to develop 
spatial variations of the original SGR trap model to in- 
clude mechanically realistic interactions between traps 
( Hebraud and Leaueux. 19981) . as is also done in theo- 



retical modelling of elasto - plastic response of am orphous 
solids (jPicard et al., 20051 iBocauet et al. 20091) . 



Mean-field glass models can also be modified to in- 
clude the physical effec t of an extern a l mechanical forc- 
ing in the dynamics ([Horner. 1996 : Thalmann. 19981 : 



ICugliandolo et al. 1997bl iBerthier et al. 2000h . Since 
these models contain fully-connected Hamiltonians de- 
fined with no reference to geometry, the modeling of an 
external flow is necessarily crude. In the case of a sheared 
glassy material, it is argued that the main effect of the 
imposed shear flow in the equations of m otion is to inject 
energy into the system ( Kurchan. 200ll ). Taking again 
the example of the p-spin model in Eq. (|4T))) . one now 
considers the driven dynamics 

dJ sr = -*«)»<»- £® +r-w+*w. do?) 

where /f rive (t) stands for an external driving force. A 



natural choice in this context is to consider a driving 
force which has a functional form similar to the p-spin 
interaction but involves coupling constants which contain 
an asymmetric part, so that the resulting force cannot 
be derive d from an energy fun ction. The specific choice 
made in ( Berthier et al. 20001 ) is 



/f™(i) = a{t) 



£ ■ 

ji<— <jk~i 



(108) 



with coupling constants that are random Gaussian vari- 
ables of variance fc!/(2iV fc_1 ), which are symmetrical 
about permutations of (ji,--- ,jk-i), but arc uncor- 
rected about permutations of i with any element of 
(ji, • • • , Jfe-i). With this particular choice of asymmet- 
rical couplings and a constant amplitude of the driv- 
ing force, a(t) = a, a numerical solution of the two- 
time dynamical equations of the form similar to Eq. ([55)1 
shows that the dynamics becom es stationary for any 
a > (jCugliandolo et al. 1997bh and any temperature 
T > (even for T < T c ). Therefore, in the stationary 
state following a quench at time t w = — oo in the pres- 
ence of a constant driving force, the dynamical equations 
become: 



dC{t) 

dt 
dR{t) 

dt 



fiR(t)+ / dt"E(t-t')R(t'), 
Jo 



nC{t)+ / dt'E(t-t')C(t') + / dt'[T,(t + t')C{t') + D(t + t')R(t')] 
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dt' [D(t')R(t') + £(t')C(t')] 



r 



with kernels given by: 



D(t) 



°\C{t) 



k-l 



E(f) = E^l_Jlc(t)P' 2 R(t). 



(110) 

(111) 



Mathematically, the asymmetry in the coupling con- 
stants of the driving force shows up in the expressions of 
the kernels, since only D(t) contains the driving term pro- 
portional to <7, with no counter part in the expr ession for 
S(t). It can be shown formally ( Kurchan. 20031 ) that de- 
tailed balance imposes specific symmetries of the kernels 
D and S, which is indeed explicitly broken in Eqs. (j 109[) 
by the term proportional to a. As in the aging regime, 
Eq. (|85)1 , these equations involve both correlation and re- 
sponse functions, and are thus more difficult to solve than 
the equilibrium, high temperature dynamics in Eq. (1491) . 
Howeve r, the solution procee ds very much as in the aging 
regime ( Berthier et al. 200 ll ). 



The connection to rheological quantities is done us- 
ing energetic considerations. Using Grecn-Kubo type 
of arguments, the viscosity r\ is related to the relax- 
ation timescale obtained from the time decay of C(t) 
from Eq. (|109[) , so that rj(T,a) ~ r a . The power dis- 
sipated by the driving force is then shown to be cx 
<7 2 /tq, ~ a 2 /i], while it is oc aj = <r 2 /r] for a sheared 
fluid. Thus one identifies a as the 'shear stress' for 
the present system, and steady state constitutive rhe- 
ological relations can readily be studied. As for the 
SGR model, the results show that the interplay be- 
tween the glass transition physics and external forcing 
induces strong shear-thinning behaviour which depends 
quantitatively on the temperature, with experimentally 
releva nt scaling laws for th e viscosity across the (a, T) 
plane ( Berthier et al. 2000h . For T < T c for instance, a 
shear-thinning behaviour is predicted at low shear rates, 



-u{T) 



(112) 
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where v(T) is a temperature dependent exponent, with 
v(T c ) = v c = 2/3 and v — > 1 as T — > 0. This behaviour 
is ty pical of a 'pow er-law fluid' in the rheological litera- 
ture (|Larson. 19991) . The non-trivial shear-thinning ex- 
ponent v c = —2/3 at T c reveals a complex interplay be- 
tween thermal processes and mechanical forcing, while in 
the low-T limit the 'natural' exponent v = —1 is recov- 
ered, as expected on dimensional grounds. Above T c , the 
following scaling form is obtained: 



r/(7,T)~?7o(T) [1 + 7/70]^ 



(113) 



where ?]o{T) is the equilibrium value of the 
viscosity and 70 is a constant. This form 
of scaling is well-known an d commonly found 
in the rheology literature ( Larson. 19991 ). and 



these predictions comp are rather well with com- 



Barrat and Berthier, 2001b 


), and experi- 


Berthier and Barrat. 2002b 


ments (Crassous et al. 20061: iBesseling et al. 20071). 



As compared to the SGR model, the equations of mo- 
tion (|109[) are mathematically more involved, even in the 
steady state, so that it is technically more difficult to 
solve for the response of the system to more complicated 
histories, although this is a technical rather than a fun- 
damental limitation. A more fundamental difference, as 
already discussed in equilibrium and aging situations, is 
the absence, at the mean-field level, of thermally acti- 
vated processes which would allow the system to visit the 
numerous deep metastable states that contribute to its 
configurational entropy between the mode-coupling and 
Kauzmann transitions. 

Just as these activated processes explain why mode- 
coupling dynamic singularities are avoided in finite di- 
mensional systems, activated processes could in principle 
affect the rheological behaviour of glassy materials. As 
discussed above, including these processes in the con- 
text of mean-field glass theory is a challenging task. A 
phenomenological description in the framework of R FOT 
theory was recently proposed ( Lubchenko. 20091) . ex- 
tending in particular the validity of Eq. (|113|) to lower 
temperatures. Not e that distinct predictions were re- 
cently obtained in ( Biroli and Bouchaud. 20091 ). Stay- 
ing a t the mean-fiel d level, however, it can be ar- 
gued ( Berthier. 2003al) that thermal activation should al- 



low the system to visit metastable states with lower free 
energy. In such a state, it can be shown that the system 
must be forced above a finite level of stress in order to 
flow. In rheological terms, this means that the presence 
of metastable states leads to materials with a static finite 
yield stress (although dynamically there no finite yield 
stress). Very recent calculations in the replica frame- 
work exposed in Sec. IIV.B~6l have confirmed the existence 
of a finite shear mod ulus originating from the ex istence 
of metastable states ( Yoshino and Mezard. 20101) . In a 



realistic description, yet to be developed, one should de- 
scribe within a single framework the effect of shear and 
thermal activation within a complex energy landscape. 



Another difference with the SGR approach concerns 
the study of effective temperatures, which again mirrors 
the mean-field results obtained in aging situations (see 
Sec. IV.B[) . Indeed, many of the results obtained in ag- 
ing systems about deviations from fluctuation-dissipation 
relations can be shown to apply to the driven case as 
well. In particular, a fluctuation-dissipation ratio can 
still be defined from Eq. (IM1) leading to the notion of 
effective temperatures for driven systems. The solution 
of Eqs. (|109p yields the time dependence of both correla- 
tion and response functions, from which a parametric FD 
plot can be built. As for aging materials, it is predicted 
that these FD plots conserve an equilibrium shape, be- 
ing piecewise linear, with each relaxation timescale being 
associated to its own value of the fluctuation-dissipation 
ratio and being interpreted in terms of an 'effective equi- 
librium' for slow degrees of freedom. An interesting fea- 
ture is that effective temperatures are predicted to oc- 
cur even above the glass transition, provided the driving 
force is large enough that non-linear rheological effects 
are observed. Thus, a deep relation between anomalous 
response functions and deviations from t hermal equilib- 
rium is established ( Berthier et al, 2000l ). 



These 
ber of 



predictions 
numerical 



were confirmed in a num- 
studies of sheared super- 



cooled liquids ( 


Berthier and Barrat, 2002a: 


Ik and Barrat, 20071: 


Barrat and Berthier. 2001b 


Berthier and Barrat. 2002bt Ono et al, 20021). and 



the main physical result, the existence of an effec- 
tive temperature describing the relaxation of sheared 
glassy systems, now forms the basis of several recent 
phenomenological description of the plastic def or- 



mation of amorphous so lids (IHaxton and Liu. 2007 



Shi et al.. 20071 



Langer and Manning. 20071 
Bouchbinder. 20081 ). We are aware of no experimental 
attempt to quantify violations of the fluctuation- 
dissipation theorem under stationary conditions created 
by a shear flow, although this should in principle be 
much easier than in non-stationary, aging situations 
where several experiments have already been performed. 



We mentioned in Sec. IV.Bl that mode-coupling theory 
had not been fully extended to deal with nonequilibrium 
aging situations. By contrast, in recent years, a large 
research activity has been dedicated to the derivation 
of mode-coupling approximations to deal with the rhe- 
ology of glassy liquids and colloidal suspensions. A first 
derivation is obtained starting from generalized fluctuat- 
ing hydrodynamic equations, as in field-theoret ic deriva- 
tions of the equilibrium MCT ( Das et al. 19851 ). Among 
the several approximations involved, the standard mode- 
coupling decomposition of four-point static correlations 
as product of two-point functions is performed, yielding 
closed dynamical equations for the time evolution of the 
intermediate scattering function under shear. For the 
case of a stationary simple shear flow with an imposed 
strain in the x-dire ction 7(0, one gets the followin g dy- 
namical equations ( Mivazaki and Reichman. 20021) : 
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M(q,f) = ^^ydq'y(q J q')F(q(t),q'(t))F(q(t)-q'(t),t)F(q'(t),t) J 



where q(t) = (q x ,qy + 'y(t)q x - ) Qz) is a time depen- 
dent wavevector, the vertex function has the standard 
mode-coupling expression as in Eq. (|47p . and F(q,t) — 
(p q (_t)(t)p-q(0)) is the intermediate scattering function 
modified to take into account the global advection by the 
shear flow. Formally, these equations are very similar to 
the ones derived at equilibrium, see Eqs. (|461 1471) . This 
implies that the physics captured by this approximation 
stems from the advection, and thus the distortion, of den- 
sity fluctuations along the x-direction by the shear flow. 
Due to the mode-coupling mechanism in Eq. 1|114|) . re- 
laxation of density fluctuations in the x-direction trig- 
gers t he relaxation of all th e modes. In a paralle l 
effort (jFuchs and Cat es. 2002t iFuchs and Cates. 20051) . 
very similar MCT dynamic equations for glasses un- 
der shear flow were derived employing the technique 
of projection operators also used at thermal equilib- 
rium to derive mode-coupling equations. Here also, a 
similar decoupling of four-point correlations into prod- 
ucts of two-p oint quantities is perfo rmed, and in its 
latest version ( Fuchs and Cates. 20091 ) the final dynam- 
ical equations ar e very similar to the ab ove expressions, 
Eqs. (fTT4"l) , see (|Fuchs and Cates. 2009]) for a detailed 
discussion of the technical differences between the two 
approaches. 

The rheological behaviour in steady shear flow above 
the mode-coupling transition resembles the description in 
Eq. (|113[) , although the 'simple' shear-thinning exponen t 
v = — 1 is found. In the glass ( Fuchs and Cates. 2002t) . 
the system develops a finite yield stress, which jumps 
discontinuously from zero when the mode-coupling tran- 
sition is crossed. Below T c , the r heological beh aviour re- 
sembles that of a Bingham fluid (jLarson. 19991) . We note 
that although fully connected glass models and schematic 
mode-coupling models are fully equivalent at thermal 
equilibrium, both approaches seem to predict different 
behaviours out of equilibrium. It is not yet clear whether 
these differences are due to the set of approximations in- 
volved in the derivation of Eqs. (|114[) . or to the very pe- 
culiar form of the non equilibrium drive (|108|) chosen for 
mean field glass models, which may be unrealistic. Both 
approaches suffer from the same fundamental limitation 
that emerging 'critical' or 'universal' properties near the 
mode-coupling transition will be drastically modified in 
realistic numerical simulations or real experiments, as the 
mode-coupling singularity is not present in real materi- 
als, even in colloidal ha rd sphere systems wh ere the the- 
ory was often applied ( Crassous et ai. 20061 ). This im- 
plies that all the problems and ambiguities due to the 



'crossover' nature of the dynamic transition encountered 
by MCT at thermal equilibrium (see Sec. IIV.B.4|) will 
be again present under shear, and the theory should fare 
with experimental data no better or no worse than at 
equilibrium. 

Although initially developed to study stationary 
shear flows, 7(t) = jt, MCT-based rheological equa- 
tions have bee n now derived for arbitrary flows and 
shear histories (|Brader et ai. 20071 : 1 rader et ai, 20081: 
iBrader et ai, 20091 ). The resulting equations are more 
complicated than Eq. (| 1 14[) . involving in particular 
three-time memory functions in the general case. The 
case of oscillatory simpl e shear flows was also stud- 
ied ( Mivazaki et al., 2006h . The possibility to study time 
dependent shear flows makes the MCT approach quite 
appealing, as it can thus compete with the flexibility of- 
fered by the SGR model described above, with the ad- 
vantage that one is dealing with a microscopically real- 
istic description of the liquid. However, it is possible to 
study t he interplay of shear flow and aging in the SGR 
model ( Fielding et ai, 2000f ). which is not yet feasible 
within MCT approaches. 

As concrete applications of the theoretical frame- 
work, step strain protocols ( Brader et ai. 20"07l ) or si- 
nus oidal shear flows of ar bitrary amplitudes were stud- 
ied (|Mivazaki et ai. 20061 ). In Fig. [2U we show the be- 
haviour obtained for the storage and loss moduli of a 
dense suspension of hard spheres close to the colloidal 
glass transition as a function of frequency (as in linear vis- 
coelasticity) and strain amplitude (as in nonlinear steady 
shear flows). These results illustrate both the broad vis- 
coelastic spectra emerging in fluids near the glass transi- 
tion, and the competition between intrinsic slow dynam- 
ics and the external perturbation imposed by the shear 
flow which accelerates the dynamics. 

The issue of a nonequilibrium extension to the 
fluctuation-dissipation theorem in shear flows was 
also discussed within MCT. A formal deriva- 
tion of the Einstein relation between self-diffusion 
and mobility under shear flow is presented in 



(ISzameL 20041). while FDT violatio ns are acknowl- 
edged ( Mivazaki and Reichman. 20021 ). but do not play 
a role during the derivation of Eqs. (| 1 14[) to arrive at 
dynamical equations involving correlation functions 
only, as opposed to coupled equations for response and 
correlations in Eqs. p09p . Very recently, an approxi- 
mate expr ession for time dependen t susceptibilities was 
obtained ( Kriiger and Fuchs. 20101 ) and compared to 
the correlation obtained from the solution of Eq. f|l 14|) . 
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FIG. 21 Frequency and strain dependence of (a) the storage 
modulus and (b) the loss modulus for a hard sphere suspen- 
sion near the mode-coupling glass transition. The broad vis- 
coelastic spectra emerging in fluids near the glass transition is 
affected in a non-linear manner by a la rge external shear flow 
which accelerates the dynamics. From (|Mivazaki et al, 20061 ) 



In standard equilibrium MCT, the response function is 
obtained a posteriori using the FDT. Complementary 
approaches, like field theory are plagued by difficulties 
related to maintaining FDT in a self-consistent manner. 
As a consequence, the s tatus of the approximation in 
( Kriiger and Fuchs. 2010() is not entirely clear at this 
stage. Moreover, an important difference is found with 
mean-field results since a nea rly constant value of the 
FDR, X = 1/2, is predicted (jKrfiger and Fuchs. 2010h 
in contrast to mean-field models where the FDR is 
generically dependent of the state point, w ith X ~ T 
at low temperatures ( Cugliandolo. 20031). as also 



found in numerica l work (Di Leonardo et al, 20001 : 
IZamponi et al, 20051: IBerthier. 2007bl ) This suggests 
that still more work is needed to clarify the status 
(derivation and physical behaviour) of non-equilibrium 
effective temperatures in glassy systems under shear. 



E. Current status of nonequilibrium studies 

Aging studies have a long history in the field of the 
glass transition ( Struik. 19781) . since glasses are all by 
definition nonequilibrium aging materials. It is only since 
the 90's, however, that a large body of work has been per- 
formed to describe nonequilibrium glasses at a more fun- 
damental level, using tools and concepts from statistical 
mechanics, as briefly reviewed above. Despite much con- 
ceptual progress, it is not clear whether these theoretical 
advances have led to a better description of the phys- 
ical properties of, say, polymer glasses under complex 
thermo-mechanical histories, and thus have had concrete 
experimental consequences (|Angell et al, 20001) . It cer- 
tainly remains true that a large number of experiments 
have been devised to analyze two-time response and cor- 
relations functions in a number of materials, in order to 
test theoretical predictions, although it appears that less 
studies of aging materials were performed in more recent 



years. 

We wish to suggest that recent theoretical and experi- 
mental work to characterize and understand memory and 
rejuvenation effects in spin glasses could be fruitfully re- 
visited in the field of structural glasses, see the experi- 
mental chapters in (lYoung. 19981) . In particular, a clear 
link between memory effects and typical lengthscales over 
which the slow dynamics takes place was established 
in the context of spin glasses (iBouchaud et aZ.. 20011: 
IBerthier and Bouchaud. 20021: |j5nsson et al, 2004) . We 
have discussed in previous sections that dynamic length- 
scales also depend very sensitively on temperature in 
structural glasses, as illustrated for instance in Fig. [TT] 
At present, it is not really known how these ideas ap- 
ply to aging of structural glasses and whether they 
would yield memory effects similar to the ones ob- 
served in spin glasses. In the same vein, very lit- 
tle work was performe d to understand how the mo- 
saic picture of RFOT (iLubchenko and Wolvnes. 2004t 
iBiroli and Bouchaud. 20091) or the frustration-limited 
domains discussed in Sec. IIV.DI would respond to com- 
plex thermal histories, nor has mode-coupling theory 
be fully developed to handle non-stationary dynamics 
far from equilib rium except in the context of m ean-field 
glassy systems (jCugliandolo and Kurchan. 19991 ). 

More work should also certainly be performed to 
characterize dynamic heterogeneities in nonequilibrium 
systems. Indeed, when experiments were devised to de- 
tect for instance two-time correlation functions in aging 
materials, it became immediately clear that slow dy- 
namics in these far from thermal equilibrium conditions 
was also ch aracterized by intriguingly large dynamic 
fluctuations ( Cipelletti and Ramos. 20051) showing up in 
the form of intermitte nt time s ignals dBellon et al, 200 1[ 
Buisson et al, 2003al iDuri and Cipelletti. 200~ ; 



Mamane et al, 20091) . very large spatial correla- 
tions (IDuri et al, 2009), peculia r forms of relaxation 
functions (ICipelletti et al, 20001: iBellour et al, 20031: 



iBandvopadhvav et al, 2004), or superd iffusive processes 
(|Guo et al, 20091: iMazover et al, 20091) . Although well- 
documented experimentally and displaying some form of 
'universality', these non-trivial dynamics have received 
too little attention from the theoreti cal community, see 
however ( Bouchaud and Pitard. 200"ll) . It is in particular 
intriguing that these effects are apparently not observed 
in computer simulations where dynamic heterogene- 
ity seems t o proceed ro ughly as under equilibrium 
conditions dParisi. 19991: Castillo and Parsaeian. 20071 ; 



iParsaeian and Castillo. 20081: |E1 Masri et al, 2010j) . 
In this context, an interesting development is the 
extension of some of the concepts derived from mean- 
field spin glass models, in particular the notion of 
time reparametrization invariance of the dynamical 
equations of motion, to finit e d imensional aging ma - 
terials (|Castillo et al, 2002t IChamon et al, 2004 . 
In particular, this approach naturally explains the 
appearance of nontrivial mcsoscopic dynamic fluctu- 
ations in aging materials, and provides specific new 
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tools to analyze a number of physical quantities, 
such as distributions of time correlation functions 
( Chamon and Cugliandolo. 20071 ). 

In striking contrast with aging studies, the rheology of 
soft glassy materials has been increasingly actively stud- 
ied over the last decade, and has developed as a research 
field on its own that will certainly continue to expand 
in the coming years. Thus, we close this prospective 
section with just a few selected issues on the rheology 
of glasses that are currently the object of intense 
research. First, since all the above mentioned theoretical 
modeling of glassy rheology are somehow 'mean-field' 
in nature, it is not clear how spatial heterogeneities 
or correlation s can be described th eoretically. Yet , 



experiments dCoussot et ai. 20021 : Becu et ai. 2006 



Ianni et al. 20081 ) and simulations (Varnik et al. 2003 ; 
Zink et al.. 20061 : IShi and Falk. 20061 ) clearly reveal for 



instance that soft glassy materials commonly display 
the phenomenon of shear-banding. Namely, when 
submitted to a macroscopic shear force, the system 
spontaneously 'phase separates' between a flowing state 
supporting the shear, and an immobile state with no 
flow. This observation means that at least two dynami- 
cal states exist for a given level of external forcing, so 
that the flow curve a(j) is multi- valued, but there is 
presently no ag reement on the micr osco pic origins of this 
o bservation (IVarnik et ai. 20031: _ iBerthier. 2003at 



Manning et al. 2007 



Besseling et al.. 201 



iFielding et ai. 200S : 
A second relevant question 
related to spatial aspects is the subject of dynamic 
heterogeneity. We have underlined in Sec.|llT]the impor- 
tance of spatial fluctuations of the dynamic at thermal 
equilibrium, and have described in Sec. IIVI how different 
theories describe these fluctuations. Much less is known 
under shear, although numerical simulations have re- 
vealed the_e2ostejice_cJJa£ge scale heterogeneities under 
flow ( Furukawa et al.. 2 009) that arc in fact strongly 
reminisce nt of the plastic flow o f v ery low temperature 
glasses dFalk and Langer 1998t iTanguv et ai. 20061: 
iLemaitre and Caroli. 2009^ Further studies of these 
issues could also shed some light on the flow of soft glassy 
materials under confinement, a situation of great interest 
due to the rapid development of microfluidic techniq ues 
and applications ( Govon et al.. 20081 : llsa et ai. 20091 ). 



VI. SOME GENERAL AND CONCLUDING REMARKS 

We conclude this long review with few remarks on some 
key topics that arc often debated in the literature and on 
which we would like to present our point of view. 



rapidly approaching T g , whereas the latter remained elu- 
sive for a long time. Recently, this state of affairs has 
changed. First, the whole topic of dynamical heterogene- 
ity made it clear that a complete theory of the glass tran- 
sition has to be able to explain growing dynamical length- 
scales. Furthermore, there are general theoretical, and 
not just phcnomenological reasons to believe that grow- 
ing lcngthscalcs should play an important role. In fact, 
a system with a finite number, say N, of degrees of free- 
dom is expected to have a relaxation time no larger than 
exp(KN), where K is independent of N and possibly de- 
pendent on temperature, but in a non singular way ex- 
cept at T = 0. A system, whose largest correlation length 
is £, can be viewed as a collection of independent sub- 
systems of linear size £. The relaxation time r is there- 
fore equal to the one of a given sub-system and there- 
fore cannot be larger than exp(A'^ <i ). This implies that a 
large relaxation times should be related to large numbers 
of spatially correlated degrees of freedom. This intuitio n 
was made rigorous in (|Montanari and Semeriian. 2006aD . 
where it is shown that indeed the point-to-set length 
defined in Sec. IIV.B.6I has to grow at least as fast as 
c(logr Q ) 1 / ti , where d is the spatial dimension and c a 
proportionality constant. Since c actually depends on 
temperature (as K in the previous expression) this in- 
equality makes a growing length a necessity only if the 
relaxation timescale diverges at a finite temperature. If 
not, one should evaluate the constant and check whether 
the inequality indeed implies a lengthscale that indeed 
becomes 'large' at low temperature. In any case, all 
these results, and many others, have been so influential 
in stressing the importance of possibly growing length- 
scales that by now many - possibly too many - length- 
scales have been defi ned and studied in the literature 
( Berthier et al.. 201 lh . The more pressing open ques- 



tions in this field concern the relation between the well- 
studied dynamical lcngthscalcs (e.g. £4), and the more re- 
cently devised static point-to-set lengthscales. Are these 
two types of spatial correlations related? Do structural 
correlations drive dynamical ones? What is their pre- 
cise relation with the overall increase of the viscosity? Is 
there a unique way of defining static correlation length- 
scales, and are point-to-set correlations equivalent to al- 
tcrnati ve, more geometric ways of defining static length- 
scales (|Tarius et al.. 20051 ICoslovich and Pastore. 20071: 
iKurchan and Levine. 20091 : iTanaka et al. 201of )? We 
believe that the intensive study of dynamical lcngthscalcs 
in the past decade will be continued by a similar intensive 
search of static correlations in future work to disentangle 
all these issues. 



A. Growing lengthscale(s) 

For a long time, the research on the glass transition 
has been focused on timescales more than lengthscales. 
The reason is simple: the former are clearly increasing 



B. Glass and jamming transitions 

We introduced in Sec. lII.BTTl thc idea that many differ- 
ent materials undergo a fluid to amorphous solid transi- 
tion reminiscent of the glass transition of molecular liq- 
uids, and indeed we included experimental or numerical 
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results on colloids or granular materials without further 
discussion in the rest of the article. It is time to discuss 
more critically the assumption that all these materials 
'jam' in a similar way. 

From a practical point of view, one should distinguish 
two distinct 'solidity' transitions that can both be ob- 
served in the system of hard spheres, which we men- 
tioned several times in this review. At thermal equi- 
librium at temperature fc^T 5 , hard spheres undergo 
in three dimensions a glass transition in the regime 
(fig « 0.57 — 0.59. Above this transition, the system 
appears as a solid, at least on experimental timescales. 
However, this system is compressible (it is a hard sphere 
glass), its pressure is finite and its equation of state, 
Z = P((p)/ pksT = Z(tp) is a smooth function of <p across 
tpg. At present, there no indication that the physics of 
this first fluid - solid transition in the hard sphere sys- 
tem at finite temperature is any different from the glass 
transition observed, say, in a Lennard- Jones liquid. This 
means that all the concepts and theories reviewed in this 
paper are actually relevant to this situation. For in- 
stance, experimental evidence suggests that hard spheres 
undergo a change from MCT-li ke behaviour to an act i- 
vated regime when increasing ip (|Brambilla et ai, 20091) . 
they display a simi lar spatially heterogeneous dynamics 
( Weeks et ai. 20001 ) . an d display in analyt ical calcula- 
tions a mode-cou pling dBarrat et ai. 19891) and Kauz- 
mann transitions dParisi and Zamponi. 20101 ) completely 
analogous to model liquids. Interestingly, the same phe- 
nomenology is found in dens e granular materials driven 
for instance by cyclic shear jMartv and Dauchot. 20051) 
or air flow (jKevs et ai, 20071 ). In these cases, a non- 
equilibrium steady state is reached thanks to the me- 
chanical driving, which plays a role similar to Brownian 
forces in colloids. When increasing the density or de- 
creasing the strength of the driving, these systems ap- 
pear to display a 'granular glass transition' with proper- 
ties which are again very similar to the ones observed for 
supercooled liquids, eve n at the most microscopic level 
(jCandelier et ai, 2009bl ). 

A second, distinct 'solidity' transition occurs in hard 
spheres out-of-equilibrium. As discussed above, it is 
not possible to compress hard spheres above a certain 
density and keep the system in equilibrium. This does 
not mean the system cannot be compressed anymore, 
it can, but in a non-equilibrium and protocol depen- 
dent manner. This second 'solidity' transition takes 
place when the system cannot be compressed anymore. 
This is mainly a geometrical problem and thermal en- 
ergy plays no role in this transition. For a three di- 
mensional system, this occurs near 'random close pack- 



0.63 - 0.65 (jBernal. 19591 ). At this den- 



5 For hard spheres, the temperature simply sets the microscopic 
timcscale since it only affects the microscopic timescale. As long 
as the system is in equilibrium, the thermodynamic and dynam- 
ical properties at different temperatures are identical up to a 
trivial rescaling. 



mg ip Icp 

sity, the compressibility vanishes, the reduced pres- 
sure Z(ip) diverges, and the number of contacts per 
particle is exactly equal to the minimal number re- 
quired for the sys tem to behave mechanically as a solid 
( Alexander. 19981 ) . This second transition is thus directly 
relevant to understand the static properties of granu- 
lar materials. It is interesting also for systems made 
of large (athcrmal) particles such as foams and emul- 
sions: because the particles are soft, these system s can 
be compressed above <y9 rcp ([Liu and Nagel. 200ll ). In 
these soft systems, the (osmotic) pressure in the solid 
phase is now proportional to the particle surface ten- 
sion, and thus solidity is driven by the elasticity of the 
particles rather than temperature. From a dynamical 
point of view, much less is known about this second 
transition. Connections with the glass problem are still 
rather speculative but they are the focus on an intense 
res earch activity, trigge red by the seminal contribution 
of ( Liu and Nagel. 19981 ). For instance, dynamic hetero- 
geneity at T = near rando m close packing has only 
been studied very recently dLechenault et ai, 2008bt 
ISessoms et al., 20091 iHeussinger et ai. 20101 ). and im- 
portant differences with the dynamics of viscous liquids 
were noted. Therefore, although glass and jamming tran- 
sitions might be observed in the same system (say, hard 
spheres), they likely correspond to two distinct ways for 
the system to become a solid. 



C. Metastability and the role of the crystal 

In this review we have discussed very little the role of 
the crystalline state. As many (but not all) researchers, 
we have assumed that the crystal does not play an im- 
portant role for the glass transition phenomenon, apart 
from the fact that crystal nucleation has to be avoided by 
supercooling. (For all known glass-formers, the melting 
temperature is larger than the experimental glass tran- 
sition.) This may be questionable for several reasons. A 
first objection is that equilibrium thermodynamic theo- 
ries of the glass transition are problematic because the 
true thermodynamic phase is indeed the crystal. This 
is not a real concern: supercooled liquids are in a long 
living metastable state. As long as the nucleation time 
is much larger than the structural relaxation time, r a , 
they can be considered as bona fide equilibrium states, at 
least using Feynman's definition: "When all fast things 
have happened and all slo w things have n ot, then the 
system is in equilibrium" ( Fevnman. 19721 ). After all, 
the people who value diamonds do not worry that they 
will turn into graphite anytime soon. A more serious 
concern is that the relaxation time of the supercooled 
liquid cannot really diverge at a finite temperature if 
the thermodynamically stable state is a crystal, because 
the nucleation time is necessarily finite at finite temper- 
ature if we assume a finite Gibbs free energy difference 
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between the crystal and the supercooled liquid 6 . Thus, 
t q will inevitably hit the (finite) nucleation time before 
diverging. Below this temperature the supercooled liq- 
uid is no more a metastable state because it nucleates 
the crystal before actually being able to relax its struc- 
ture. This, however, is also not necessarily a serious 
problem. First, not all theories are based on a diver- 
gence at finite temperature. Second, several theories ex- 
plain the slowing down of the dynamics by the proxim- 
ity to a phase transition but none of them needs that 
the transition actually takes place to be proven correct. 
Thus, although important conceptually, the existence of 
the crystal does not imply that theories with no singu- 
larity should be preferred to describe the physically rel- 
evant temperature regime where nucleation is unimpor- 
tant. A final, more physical, reason to take into account 
the crystal in explaining the glass transition would be 
that the slow dynamics in the supercooled state is due to 
the existence of some local order reminiscent of the crys- 
tal structure, as in fr ustration-limited domain approaches 
( Tarius et al., 20051 ). and reemphasized in several re- 
cent numerical papers (ICoslovich and Pastore. 2007t 
iTanaka et al. 20ld IPedersen et al. 2010h . 



D. The 'ideal' glass transition 

This is certainly a recurrent, probably Utopian, but 
nevertheless fascinating, topic for discussion. Several 
theories explain the slowing down of the dynamics by 
the proximity to a phase transition. None of them re- 
quires that the transition actually takes place. This 
holds either because the transition is avoided by con- 
struction as in the frustration-limited domain theory, or 
because is not accessible experimentally as in dynami- 
cal facilitation theory or random first order transition 
theory. What happens at lower temperatures, where 
no real system can be equilibrated, is after all just a 
matter of curiosity. It is however interesting, just for 
a short while, to dwell on what an ideal glass transi- 
tion could possibly be. It should correspond to a true 
divergence of the relaxation time and the viscosity at 
a finite temperature Tidcai- Because of the proof by 
Montanari and Semerjian for lattice models, we now 
strongly suspect that static correlations, actually point- 
to-set c orrelations, would generally diverg e at the tran- 
sition 7 (jMontanari and Semerjian. 2006b ). Thus, the 



6 This is due to the fact that the maximum barrier to form the criti- 
cal nucleus and also for its subsequent growth is necessarily finite, 
thus nuclei will form and expand even though these processes can 
be extremely slow. In practice, crystal nu cleation can be much 
slower than what is naively expected, see fCavagiia, 200JJ). 

7 If the relaxation timescale only diverges at zero temperature but 
faster than Arrhenius then the ideal glass transition would take 
place at T = 0. Again, on the basis of the bound of Montanari 
and Semerjian, we would expect a diverging length in this case 
too. 



transition would correspond to the development of some 
long-range order (likely amorphous order), since suitable 
boundary conditions fix the density field in an amorphous 
configuration in the entire (infinite) sample. 

It is also important to discuss what properties an ideal 
glass transition would not display. Actually, it is some- 
times assumed that at an ideal glass transition, as dis- 
cussed in the context of theories in Sec. IIV1 the dynam- 
ics would be completely arrested. This is not necessary 
and, actually, impossible if Tidcai > 0. For example, a 
probe particle will still be able to move at Xideai and be- 
low. The situation would be similar to the crystalline 
state, where particles diffuse, although very slowly, leav- 
ing the crystalline order intact. What would diverge at 
Ildeai is only the time to destroy the amorphous correla- 
tion in the density field. For example, at the ideal glass 
transition advocated in RFOT theory, the density field 
orders in one of the possible amorphous low temperature 
configurations, so that the time to relax the structure 
becomes infinite whereas, instead, the self-diffusion coef- 
ficient stays finite. Also, it is sometimes believed that the 
ideal glass transition is related to a fragmentation of the 
configuration space for finite size systems, i.e. that below 
a certain temperature or above a certain density it is no 
more possible to go from any typical equilibrium config- 
uration to another one. This is clearly not true, as it can 
be seen easily for soft particle models, and actually even 
for hard spheres. This can be harder to prove for some 
effective models ( Eckmann and Procaccia. 2008| ). In any 
case, this is not what the ideal glass transition would 
be: such kind of dynamical arrests can only take place 
in effective theories whose domain of applicability must 
break down at a some point (maybe not accessible in 
experiments). 

Thus, the theoretical possibility of a finite temperature 
'ideal' glass transition towards a genuine 'glass' state ex- 
ist, and will no doubt continue to obsess many physicists 
in coming years. 



E. Final words 

The problem of the glass transition, already very ex- 
citing in itself, has ramifications well beyond the physics 
of supercooled liquids. Glassy systems figure among the 
even larger class of 'complex systems'. These are formed 
by a set of interacting degrees of freedom that show an 
emergent behaviour: as a whole they exhibit properties 
not obvious from the properties of the individual parts. 
As a consequence the study of glass-formers as statisti- 
cal mechanics models characterized by frustrated inter- 
actions is a fertile ground to develop new concepts and 
techniques that will likely be applied to other physical, 
and more generally, scientific situations. The glass tran- 
sition in supercooled liquids can in fact be considered as 
one of the simplest situations where frustration, geome- 
try, ergodicity, and disorder compete to produce a glassy 
state. Hence, the concepts reviewed in this article should 
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be directly applicable to scores of more complex systems, 
such as for instance physical gels, liquids in confined ge- 
ometries, diffusion in crowded biological environments, 
dense granular media, self-assembly, or microfluidic flows 
of dense emulsions and colloidal suspensions. Thus, we 
certainly expect more progress to emerge in the future 
along these interdisciplinary routes. 
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